[noip2014]解方程 题解

看看数据范围，woc，什么鬼，让我在初中第一眼看见的时候惊呆了，然后留下了阴影。现在才补完。

其实也不算太难。

首先，我最先先到的就是在一个模大质数的意义下解，复杂度就成O(mn)的了，再用两个比较大的质数同时做，如果这样都wa的话就没办法了（几率太小）。不过，这样是过不了最后两个点的，怕是这样只能O(wys)了233。

但是，我们可以改进一下思路，先选取一个比较小的质数p，算它的答案，假设为x，那么（x+p）也是一个解。显然，这样的解不超过n个，我们再用一个比较大的数在这些解里面判断是否为真正的解就好了。

#include<bits/stdc++.h>#define LL long longusing namespace std;const LL p[]={67891,1000000207};char s[10005];LL a[10][105];int tot;LL n,m,st,ed,x;bool flg,E[1000005];inline void write(int x){    if(x==0)putchar('0');    else    {        char buf[15];        int len=0;        if(x<0)putchar('-'),x=-x;        while(x)buf[++len]=x%10+'0',x/=10;        for(int i=len;i>=1;i--)putchar(buf[i]);    }    putchar('\n');}bool check(LL x,int k){    //for(int k=0;k<5;k++)    //{        LL tmp=1,sum=0;        for(int i=0;i<=n;i++)        {            sum=(sum+tmp*a[k][i])%p[k];            tmp=(tmp*x)%p[k];        }        if(sum!=0)return false;    //}    return true;}int main(){    freopen("in.txt","r",stdin);    cin>>n>>m;    for(int i=0;i<=n;i++)    {        flg=false;        scanf("%s\n",s);st=0,ed=strlen(s);        if(s[0]=='-')st++,flg=true;        for(int k=0;k<2;k++)        {             x=0;            for(int j=st;j<ed;j++)            {                x=x*10+s[j]-'0';                x=x%p[k];            }            if(flg)a[k][i]=(p[k]-x);            else a[k][i]=x;           }    }    for(int i=1;i<=p[0];i++)        if(check(i,0))            for(int j=i;j<=m;j+=p[0])                if(check(j,1))E[j]=true,tot++;    write(tot);    for(int i=1;i<=m;i++)if(E[i])write(i);    return 0;}