AGC 027D.Modulo Matrix(构造 黑白染色)

题目链接

---

\(Description\)

给定\(n\)，要求构造一个\(n\times n\)的矩阵，矩阵内的元素两两不同，且任意相邻的两个元素\(x,y\)，满足\(\max(x,y)\ \mathbb{mod}\ \min(x,y)\)等于一个非零常数。
\(n\leq500,\ 1\leq 矩阵中的元素\leq10^{15}\)。

\(Solution\)

https://blog.csdn.net/Tiw_Air_Op1721/article/details/82719507 orz!

考虑没有元素互不相同的限制，那么就直接填\(x,x+1,x,x+1...\)，即对矩阵黑白染色，黑格子填\(x\)，白格子填\(x+1\)。
如果有限制，依旧可以先黑白染色，然后令白格子比周围四个黑格子都大。就令\(\max(x,y)\ \mathbb{mod}\ \min(x,y)=1\)好了，这样白格子的权值等于周围四个黑格子的\(LCM+1\)。
但是这样直接随便给黑格子赋权值，元素是会超过\(10^{15}\)的。
因为是求\(LCM\)，不妨令黑格子根据对角线设成两个数的乘积。即，对于左上-右下和左下-右上这\(2n\)条对角线（只考虑黑格子），分别给它一个数。那一个黑格子的它所在的两条对角线上的数的乘积。
要保证任意元素互不相同，所以考虑把这\(2n\)个数设成\(2n\)个质数。
那么对于最大的白格子，权值是\(第499个质数\times第500个质数\times第999个质数\times第1000个质数<8\times10^{14}\)。那么就OK啦。

\(n=2\)时不对，要特判。

---

//54ms	4992KB
#include <cstdio>
#include <algorithm>
typedef long long LL;
const int N=505,M=8100;

int P[1005];
LL A[N][N];

void Init()
{
	static bool notP[M];
	for(int cnt=0,i=2; cnt<1000; ++i)
	{
		if(!notP[i]) P[++cnt]=i;
		for(int j=1,v; j<=cnt&&(v=i*P[j])<M; ++j)
		{
			notP[v]=1;
			if(!(i%P[j])) break;
		}
	}
}
inline LL LCM(LL x,LL y)//longlong
{
	return x/std::__gcd(x,y)*y;
}

int main()
{
	Init();
	int n; scanf("%d",&n);
	if(n==2) return printf("4 7\n23 10\n"),0;
	for(int i=0; i<=n+1; ++i) A[0][i]=1, A[n+1][i]=1, A[i][n+1]=A[i][0]=1;
	for(int i=1; i<=n; ++i)
		for(int j=i&1?1:2; j<=n; j+=2)
			/*if(!(i+j&1))*/ A[i][j]=P[i+j>>1]*P[(i-j+n+1)/2+n];
	for(int i=1; i<=n; ++i)
		for(int j=i&1?2:1; j<=n; j+=2)
			/*if(i+j&1)*/ A[i][j]=LCM(LCM(A[i-1][j],A[i][j-1]),LCM(A[i+1][j],A[i][j+1]))+1;
	for(int i=1; i<=n; ++i,putchar('\n'))
		for(int j=1; j<=n; ++j) printf("%lld ",A[i][j]);

	return 0;
}