奇数阶幻方

奇数阶幻方有公式的，
第一行正中填1，
接下来斜右上45度方向填下一个数，
如下情况解决（3阶举例）：
1.行出界了，在所在列的最下方（最大行）填，如2就是在第3列第3行
2.列出界了，在所在行的最左方（最小列）填，如3就是在第1列第2行
3.碰到右上方已有数字，那就填正下方，如3右上方就有1，那就再3的下边填4，
4.行列都出界了，填在上一个数的正下方，如6是最右上角，7就是行列都出界，填在6的下方，第3列第2行

5.8，9正好按如上规则填在空白处。

#include <iostream> using namespace std; void range(int **a, int N) { if(N % 2 == 0) cerr << "Please input an odd number!/n"; int x = 0; int y = N /2; a[x][y] = 1; for(int i = 2; i <= N * N; i++) { x--; y++; if(x < 0) { if(y >= N) // row and col over { x += 2; y--; } else // row overflow { x = N - 1; } } else { if(y >= N) // col overflow { y = 0; } else { if(a[x][y] != 0) // the position is used { x += 2; y--; } } } a[x][y] = i; } } bool check(int **a, int N) { long M = N * (N * N + 1) / 2; for(int i = 0; i < N; i++) { long sumofrow = 0; // check each row long sumofcol = 0; // check each col for(int j = 0; j < N; j++) { sumofrow += a[i][j]; sumofcol += a[j][i]; } if(sumofrow != M || sumofcol != M) return false; } int sumofdiag = 0; // check the diagonal int sumofbackdiag = 0; // check the back-diagnoal for(int i = 0; i < N; i++) { sumofdiag += a[i][i]; sumofbackdiag += a[N - i - 1][i]; } if(sumofdiag != M || sumofbackdiag != M) return false; // OK return true; } void show(int **a, int N) { for(int i = 0; i < N; i++) { for(int j = 0; j < N; j++) cout << "--------"; cout << "-/n|"; for(int j = 0; j < N; j++) cout << " |"; cout << "/n|"; for(int j = 0; j < N; j++) { cout.width(5); cout << a[i][j] << " |"; } cout << "/n|"; for(int j = 0; j < N; j++) cout << " |"; cout << "/n"; } for(int j = 0; j < N; j++) cout << "--------"; cout << "-/n"; } int main() { int N = 0; cin >> N; //------------new int **p = new int *[N]; for(int i = 0; i < N; i++) { p[i] = new int[N]; for(int j = 0; j < N; j++) { p[i][j] = 0; } } //------------range range(p, N); //------------check bool is = check(p, N); if(is) { cout << "the range of " << N << " magic square:/n"; //------------show show(p, N); } else { cout << "Sorry, can't get the range of " << N << " magic square:/n"; } //------------detete for(int i = 0; i < N; i++) { delete [] p[i]; } delete [] p; return 0; }