洛谷 2312 / bzoj 3751 解方程——取模

题目：https://www.luogu.org/problemnew/show/P2312

　　　https://www.lydsy.com/JudgeOnline/problem.php?id=3751

惭愧……先听了同学讲它，又看了题解，才A。

取一些模数就能限制时间复杂度并有一些概率正确。取模后常数也能取模，就不用高精度了！而且m以内的数只用算 0~模数-1 的值，表示取模后是该值的数在模该值意义下答案是否为0。

bzoj上的数据需要精心选取模数。据说2e4左右的效果最好？试了几次终于A了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,M=1e4+,A=1e6+,K=,C=,D=4e4+;
int n,m,a[N][M],nm[N],c[N],ans[A],cnt;
bool fx[N],flag;
int zs[K]={,,,,,,};
bool vis[K][D];
void solve(int p)
{
    for(int i=;i<=n;i++)
    {
        c[i]=;
        for(int j=;j<=nm[i];j++) c[i]=(c[i]*+a[i][j])%p;
        if(fx[i]) c[i]=(-c[i]+p)%p;
    }
}
void calc(int bh,int x,int p)
{
    int cr=;
    for(int i=n;i>=;i--)
        cr=(cr*x%p+c[i])%p;
    vis[bh][x]=(!cr);
}
int main()
{
    scanf("%d%d",&n,&m); char ch=getchar();
    for(int i=;i<=n;i++)
    {
        while(ch>''||ch<''){ if(ch=='-') fx[i]=; ch=getchar();}
        while(ch>=''&&ch<='') a[i][++nm[i]]=ch-'',ch=getchar();
    }
    for(int i=;i<C;i++)
    {
        solve(zs[i]);
        for(int j=;j<zs[i]&&j<=m;j++)
            calc(i,j,zs[i]);
    }
    for(int i=;i<=m;i++)
    {
        flag=;
        for(int j=;j<C;j++)
            if(!vis[j][i%zs[j]]) {flag=;break;}
        if(!flag) ans[++cnt]=i;
    }
    printf("%d\n",cnt);
    for(int i=;i<=cnt;i++) printf("%d\n",ans[i]);
    return ;
}