[CF338D]GCD Table（扩展中国剩余定理）

题目：

我是超链接

题意：

一个数表，其中G(i,j)=gcd(i,j)，给出一个序列a1…ak，判断这个序列是否在数表中出现过

题解：

首先可以确定的是行的范围至少是lcm(a1..ak)，因为作为行必须让(a1..ak)都成为ta的因数嘛，然后n小于这个值直接退出
然后假设我们从第x列开始，不难列出
x=a1*b1
x+1=a2*b2
x+2=a3*b3
…
x+k-1=ak*bk
把ta化成同余方程组，即x=a[i]-i+1(mod a[i])
然后扩展中国剩余定理，然后再把x带回去验证一下
因为乘出来太大了，我们使用快速乘

代码：

#include <cstdio>
#define LL long long
using namespace std;
const int NN=10005;
LL N,M,a[NN],m[NN],c[NN];int k;
LL gcd(LL a,LL b){if (!b) return a;else return gcd(b,a%b);}
void exgcd(LL a,LL b,LL &x,LL &y){if (!b) x=1,y=0;else exgcd(b,a%b,y,x),y-=a/b*x;}
LL ksj(LL a,LL k,LL mod)
{
    LL ans=0;
    if (k<0) k=-k,a=-a;
    for (;k;k>>=1,a=(a+a)%mod)
      if (k&1) ans=(ans+a)%mod;
    return ans;
}
LL inv(LL a,LL b)
{
    LL x,y;
    exgcd(a,b,x,y);
    x=(x%b+b)%b;
    if (!x) x+=b;
    return x;
}
int main()
{
    scanf("%I64d%I64d%d",&N,&M,&k);
    LL lcm=1;
    for (int i=1;i<=k;i++) 
    {
        scanf("%I64d",&a[i]);
        lcm=lcm/gcd(a[i],lcm)*a[i];
        if (lcm>N) {printf("NO\n");return 0;}
    }
    for (int i=1;i<=k;i++) m[i]=a[i],c[i]=a[i]-i+1;
    for (int i=2;i<=k;i++)
    {
        LL d=gcd(m[i-1],m[i]),m1=m[i-1],m2=m[i],c2=c[i],c1=c[i-1];
        if ((c2-c1)%d) {printf("NO\n");return 0;}
        m[i]=m1/d*m2;
        c[i]=inv(m1/d,m2/d);
        c[i]=ksj(c[i],c2-c1,m2);
        c[i]=ksj(c[i],m1/d,m[i]);
        c[i]=(c[i]+c[i-1])%m[i];
        c[i]=(c[i]%m[i]+m[i])%m[i];
    }
    if (!c[k]) c[k]+=m[k];
    if (c[k]+k-1>M) {printf("NO\n");return 0;}
    for (int i=0;i<k;i++)
      if (gcd(c[k]+i,lcm)!=a[i+1]) {printf("NO\n");return 0;}
    printf("YES\n");
}