Atcoder agc027D

我们令$m=1$，对网格黑白染色后，要求白格的权值比相邻的黑格大，并且模相邻的黑格余$1$。
不考虑$a_{i,j}\le 10^{15}$的限制的话，有一个显然的构造，令黑格取不同的奇质数，白格取相邻黑格的$lcm+1$，这样显然互不相同且合法，问题是权值太大了。
考虑调整一下这个做法，我们给每条主对角线和反对角线分别对应一个奇质数，每个黑格是对应两条对角线的奇质数的乘积。那么这样白格的权值至多是四个质数的乘积$+1$，稍微注意一下质数的顺序即可通过。
时间复杂度$\mathcal{O}(n^2\log V)$。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

inline ll lcm(ll x,ll y) {
  return x/__gcd(x,y)*y;
}

int prime[10005];
bool check[10005];

void pre(int n) {
  int tot=0;
  for(int i=2;i<=n;i++) {
  	if (!check[i]) prime[++tot]=i;
  	for(int j=1;j<=tot&&i*prime[j]<=n;j++) {
  		check[i*prime[j]]=1;
  		if (i%prime[j]==0) break;
	  }
  }
}

ll ans[505][505],val1[1005],val2[1005];

int main() {
  pre(1e4);
  int n;
  scanf("%d",&n);
  if (n==2) {
  	puts("4 7");
  	puts("23 10");
  	return 0;
  }
  int cnt=1;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      if ((i+j)&1) {
      	if (!val1[i-j+n]) val1[i-j+n]=prime[++cnt];
      	ans[i][j]=val1[i-j+n];
	  }
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      if ((i+j)&1) {
      	if (!val2[i+j]) val2[i+j]=prime[++cnt];
      	ans[i][j]*=val2[i+j];
	  }
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      if (!((i+j)&1)) {
      	ans[i][j]=1;
      	if (i>1) ans[i][j]=lcm(ans[i][j],ans[i-1][j]);
      	if (i<n) ans[i][j]=lcm(ans[i][j],ans[i+1][j]);
      	if (j>1) ans[i][j]=lcm(ans[i][j],ans[i][j-1]);
      	if (j<n) ans[i][j]=lcm(ans[i][j],ans[i][j+1]);
      	ans[i][j]++;
	  }
  for(int i=1;i<=n;i++) {
    for(int j=1;j<=n;j++) printf("%lld ",ans[i][j]);
    printf("\n");
  }
  return 0;
}