C語言求所有四階幻方？

#include <stdio.h>#include <stdlib.h>#include <windows.h>#define EVEN_DOUBLE_4 4 //雙偶的最基本類型，4階雙偶#define SCREEN_SIZE 19 //屏幕顯示不變形的最大尺寸（主要是因為窗口大小限制）#define MIN_SIZE 3 //最小階數為3#define MAX_SIZE 30/*原則上是任意階，算法是相同的，這里就以30為上限吧，當然如果你愿意，可以修改的更大一些*/#define PRINT printf("Esc退出，Enter繼續"); //打印的宏#define CLEAR row = 0; column = 0; //清零int Magic[MAX_SIZE][MAX_SIZE]={0}; //全局，幻方數組int row = 0; column = 0; //全局，幻方的行列數int main(void){int read(); //讀取函數void odd(int size, int ini_value ); //奇數階幻方生成void mean_double_4(int size); //4階雙偶生成void mean_double(int size); //雙偶生成void mean_single(int size); //單偶生成void print_magic(int size); //打印幻方void sum_print(int data[], int size); //行、列、對之和打印void clear_sum(int data[]);void check_magic(int data[], int size ); //檢查所得矩陣是否為幻方陣int size; //幻方階數int sum[2*MAX_SIZE+2] = {0};//行、列、對之和do{CLEARclear_sum(sum);size=read();system("cls");if(size%2 != 0 ) odd(size, 0);else if(size %4 == 0) mean_double(size);else mean_single(size);print_magic(size);sum_print(sum, size);check_magic(sum,size);PRINT}while (getch() != 27);return 0;}/*讀入數據*/int read(){int min_size = MIN_SIZE;int max_size = MAX_SIZE;int size;do{printf("請輸入幻方階數n，n∈[%d,%d]\n", min_size, max_size);scanf("%d", &size);getchar();if(size<3 || size > MAX_SIZE) printf("非法輸入，請重新輸入[%d,%d]的正整數\n", min_size, max_size);else if(size > SCREEN_SIZE){printf("大于屏顯最大階數，輸出將變形\n");Sleep(2000);}}while(size < MIN_SIZE || size > MAX_SIZE);return size;}/*奇數階幻方，采用house法*/void odd(int size, int ini_value){int num; //填充數字int min_num = 1+ini_value;int max_num = 1+size*size+ini_value; //填充范圍int x = size/2;int y = 0; //初始填充坐標for(num = min_num; num < max_num; num++){Magic[y+row][x+column] = num;if(num % size == 0) y++; //跳步else x++, y += 2; //馬步，其實Horse法和Siamese是完全類似的x = x % size ;y = y % size ;/*越界反彈，即觸碰到邊界，從另一邊返回*/}}/*雙偶數階幻方，采用對調法*//*對調法的基礎，4階雙偶。注意要求是將非主副對角線上的元素對調，其實換個角度，你也可以說就是把祝福對角線中心對調。只不過兩種思路得到的矩陣正好是反著來的*//*本函數實現主副對角線互調*/void mean_double_4(int size){int i;int total = size * size +1 ;for(i = 0; i < EVEN_DOUBLE_4; i++){Magic[row+i][column+i] = total - Magic[row+i][column+i];Magic[row+i][ EVEN_DOUBLE_4+column-i-1] =total - Magic[row+i][ EVEN_DOUBLE_4+column-i-1];}}/*任意階雙偶*/void mean_double(int size){int num; //填充數字int min_num = 1;int max_num = 1+size*size; //填充范圍int i = 0; //循環變量