任意阶幻方..

/*@coder Gxjun*/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define maxn 100
int map[maxn][maxn] ;
void creat_magic(int n,int x,int y ,int sn)   //奇阶幻方构造
{
    int i,j,k;
    i=0;
    j=n/2;
    for(k=n;k<=n*n;k++)
    {
        map[i+x][j+y]=k+sn;
        if(k%n!=0)
        {
            if(i!=0)  i--;
            else      i=n-1;
            if(n==j+1) j=0;
            else      j++;
        }
        else
        {
            if(i==n-1)i=0;
            else    i++;
        }
    }
}

void  magic_4(int n)
{
    int i,j;
    for(i=0 ; i<n ;i++)
    {
        for(j=0 ; j<n ;j++)
        {
            if((i%4==0||i%4==3)&&(j%4==0||j%4==3)||(i%4==1||i%4==2)&&(j%4==1||j%4==2))
            {
                map[i][j]=n*n-(i*n+j);
            }
            else
            {
                map[i][j]=i*n+j+1;   //i*n+j+n
            }
        }
    }
}

void magic_other(int n)
{
    int i,j,t;
    creat_magic(n/2,0,0,0);
    creat_magic(n/2,n/2,n/2,n*n/4);
    creat_magic(n/2,0,n/2,n*n/2);
    creat_magic(n/2,n/2,0,n*n/4);
    for(i=0;i<n/2;i++)
    {
        for(j=0;j<n/4;j++)
        {
            if(i!=n/4||j!=0)
            {
                /*<swap>*/
                t=map[i][j];
                map[i][j]=map[i+n/2][j];
                map[i+n/2][j]=t;
            }
        }
    }
            t=map[n/4][n/4];
            map[n/4][n/4]=map[n/4+n/2][n/4];
            map[n/4+n/2][n/4]=t;
        for(i=0;i<n/2;i++)
        {
            for(j=n-n/4+1 ;j<n;j++)
            {
                t=map[i][j];
                map[i][j]=map[i+n/2][j];
                map[i+n/2][j]=t;
            }
        }
}
int main()
{
    int n,i,j;
    while(scanf("%d",&n)!=EOF)
    {
        memset(map,0,sizeof(map));
        if(n&1) creat_magic(n,0,0,0);
        else if(n%4==0)    magic_4(n);
        else magic_other(n);
        for(i=0;i<n;i++)
        {
            for(j=0;j<n;j++)
            {
                printf("%d ",map[i][j]);
            }
            putchar(10);
        }
    }
    return 0;
}