codeforces 338D - GCD Table 扩展中国剩余定理

题意：给一个n*m的表G，和一个含k个整数的M序列。 (1 ≤ n, m ≤ 1012; 1 ≤ k ≤ 10000).
表中每个元素值为G[i][j]=GCD(i,j).判断表G中是否存在一行从某列开始的元素和所给的序列重合
链接：http://codeforces.com/problemset/problem/338/D
首先如果存在，那么一定出现在I=lcm(M1,M2,…,Mk) ，cf上的证明没怎么看懂。如果I>n 可以直接判NO了。
对j：
j=0 mod m1
j+1=0 mod m2
.
.
j+k-1=0 mod mk

=>j=1-i mod mi
中国剩余定理解出最小的j代入原始条件判断，令LCM=lcm(m1,m2,…,mk）理论上通解为J=j+LCM*x(x>=0),GCD(i,j)=GCD(i,j+LCM*x)=GCD(i,j+i*x)=GCD(i,j). 所以判断j就足够了。

#include<cstdio> // j=0 mod m1   =>  J=m1*x0+a1 mod lcm(m1,m2)
#include<queue>  // j=-1mod m2
#include<iostream>//j=-2mod m3
#include<vector>
#include<map>
#include<cstring>
#include<string>
#include<set>
#include<stack>
#include<algorithm>
#define cle(a) memset(a,0,sizeof(a))
#define inf(a) memset(a,ox3f,sizeof(a))
#define ll long long
#define Rep(i,a,n) for(int i=a;i<=n;i++)
using namespace std;
const int INF = ( 2e9 ) + 2;
const int maxn = 1e4+10;
ll A[maxn],M[maxn];
ll n,m,k;
ll exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(!b){x=1;y=0;d=a;}
    else 
    {
        exgcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
ll Chinese_remainder(ll *A,ll *M) // 扩展中国剩余定理：通解为:LA+k*LM (k>=0) 
{
    ll d,x,y,LM=M[1],LA=A[1];
    for(int i=2;i<=k;i++)
    {
        ll a=LM,b=M[i],c=A[i]-LA;
        exgcd(a,b,d,x,y);
        if(c%d)return -1;
        ll mod=b/d;
        x=x*c/d;
        x=(x%mod+mod)%mod;
        LA = LM*x+LA;
        LM = LM*b/d;
    }
    return LA;  // 最小解 
}
bool check()
{
    ll lcm=1;
    for(int i=1;i<=k;i++)lcm=lcm/__gcd(lcm,M[i])*M[i];
    if(lcm>n)return false;      //  至少得在lcm行上 
    for(int i=1;i<=k;i++)A[i]=(1-i);
    ll J=Chinese_remainder(A,M);        //中国剩余定理解最小的j 
    if(J==0)J=lcm;  // 解为0？ 
    if(J>m-k+1||J<=0)
    {
        return false;
    }
    for(int i=0;i<k;i++)
        if(__gcd(lcm,(ll)J+i)!=M[i+1])return false;
    return true;
}
int main()
{
//  freopen("in.txt","r",stdin);//
    scanf("%I64d%I64d%I64d",&n,&m,&k);
    for(int i=1;i<=k;i++)scanf("%I64d",&M[i]);
    if(check())printf("YES\n");
    else printf("NO\n");
}