【模拟】Codeforces 710C Magic Odd Square

题目链接：

　　http://codeforces.com/problemset/problem/710/C

题目大意：

　　构造一个N*N的幻方。任意可行解。

　　幻方就是每一行，每一列，两条对角线的和都相等。

题目思路：

　　【模拟】

　　分为奇幻方、单偶幻方和双偶幻方三种构造。

　　具体分类可以查看百度。幻方的N种构造方法

//
//by coolxxx
//#include<bits/stdc++.h>
#include<iostream>
#include<algorithm>
#include<string>
#include<iomanip>
#include<map>
#include<memory.h>
#include<time.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
//#include<stdbool.h>
#include<math.h>
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
#define abs(a) ((a)>0?(a):(-(a)))
#define lowbit(a) (a&(-a))
#define sqr(a) ((a)*(a))
#define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define eps (1e-8)
#define J 10
#define mod 1000000007
#define MAX 0x7f7f7f7f
#define PI 3.14159265358979323
#define N 54
using namespace std;
typedef long long LL;
int cas,cass;
int n,m,lll,ans;
int a[N][N];
void work1(int xx,int yy,int nn,int s,int t)
{
	int i,x,y;
	x=xx+0;y=yy+nn/2;
	a[x][y]=s;
	for(i=s+1;i<=t;i++)
	{
		x=(x-1+nn)%nn+xx;
		y=(y+1)%nn+yy;
		if(a[x][y])x=(x+2)%nn+xx,y=(y-1+nn)%nn+yy;
		a[x][y]=i;
	}
}
void work2()
{
	int i,j,k;
	m=n/2;k=m/2;
	work1(0,0,m,1,sqr(m));
	work1(m,m,m,sqr(m)+1,sqr(m)*2);
	work1(0,m,m,sqr(m)*2+1,sqr(m)*3);
	work1(m,0,m,sqr(m)*3+1,sqr(m)*4);
	for(i=0,j=m+m/2;i<m;i++)
		swap(a[i][j],a[i+m][j]);
	for(i=0;i<m;i++)
	{
		if(i==k)continue;
		for(j=0;j<k;j++)
			swap(a[i][j],a[i+m][j]);
	}
	for(i=0;i<k;i++)
		swap(a[k][i+k],a[k+m][i+k]);
}
void work3()
{
	int i,j,k;
	m=n*n+1;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			a[i][j]=i*n+j+1;
	for(i=0;i<n;i+=4)
	{
		for(j=0;j<n;j+=4)
		{
			for(k=0;k<4;k++)
			{
				a[i+k][j+k]=m-a[i+k][j+k];
				a[i+k][j+3-k]=m-a[i+k][j+3-k];
			}
		}
	}
}
int main()
{
	#ifndef ONLINE_JUDGE
//	freopen("1.txt","r",stdin);
//	freopen("2.txt","w",stdout);
	#endif
	int i,j,k;
//	for(scanf("%d",&cas);cas;cas--)
//	for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
//	while(~scanf("%s",s+1))
	while(~scanf("%d",&n))
	{
		if(n&1)
			work1(0,0,n,1,sqr(n));
		else if(n%4==0)
			work3();
		else
			work2();
		for(i=0;i<n;i++)
		{
			for(j=0;j<n;j++)
				printf("%d ",a[i][j]);
			puts("");
		}
	}
	return 0;
}
/*
//

//
*/