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;
}

原文地址：https://www.cnblogs.com/SovietPower/p/10455399.html