[atAGC027D]Modulo Matrix

对网格图黑白染色，在黑色格中填不同的质数，白色格中填相邻黑色格的lcm+1，但这样会超过1e15的上限
将网格图划分为两类对角线，每一条对角线选一个质数，然后每一个点就是两条对角线的质数相乘，而白格的值就仅为四个较小质数的乘积+1（注意不能让两个大质数配到一起）

#include<bits/stdc++.h>
using namespace std;
#define ll long long
map<int,int>mat;
int n,p[10005],vis[10005];
ll a[505][505];
ll gcd(ll x,ll y){
    if (!y)return x;
    return gcd(y,x%y);
}
ll lcm(ll x,ll y){
    if ((!x)||(!y))return x+y;
    return x/gcd(x,y)*y;
}
int nex(int x){
    if (!mat[x])mat[x]=p[p[0]--];
    return mat[x];
}
int main(){
    scanf("%d",&n);
    if (n==2){
        printf("4 7
23 10
");
        return 0;
    }
    for(int i=2;i<=10000;i++){
        if (!vis[i])p[++p[0]]=i;
        for(int j=1;(j<=p[0])&&(i*p[j]<=10000);j++){
            vis[i*p[j]]=1;
            if (i%p[j]==0)break;
        }
    }
    for(int i=1;i<=p[0]/2;i++)swap(p[i],p[p[0]-i+1]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if ((i+j)%2==0)a[i][j]=nex(i+j);
    mat.clear();
    for(int i=n;i;i--)
        for(int j=n;j;j--)
            if ((i+j)%2==0)a[i][j]*=nex(i-j);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if ((i+j)%2)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++){
        for(int j=1;j<=n;j++)printf("%lld ",a[i][j]);
        printf("
");
    }
}