NOIP2014 解方程

30’

$0<n \le 2,|ai| \le 100,an \not=0,m<100$
显然$n\in\{1,2\}$
直接求根+判断：

$n=1$：

$a_0+a_1x=0$
$x=- \dfrac{a_0}{a_1}$

$n=2$ :

$a_0+a_1x+a_2x^2=0$
$x=\dfrac{-a_1\pm\sqrt{a_1^2-4a_0a_2}}{2a_2}$
复杂度：$O(1)$

50’

直接枚举$i\in[1,m]$
高精判断是否为$0$。
复杂度：$O(nm\ log_{10}a_{max})$

70’

假设$f(x)=a_0+a_1x+a_2x^2+...a_nx^n$
原题就是求$f(x)=0$的解
若有$i \in [1,m]$使得$f(i)=0$
那么：$f(i)\ mod\ k=0$（$k$ 为任意数）
$f(i)\ mod\ k=(a_0+a_1i+a_2i^2+...+a_ni^n)\ mod\ k$
$=a_0\ mod\ k+a_1\ mod\ k \times i \ mod\ k+...+a_n\ mod\ k\times i^n\ mod\ k=0$
我们在读入时取模，可以不用高精，但是考虑到：
$f(i)\ mod\ k=0$
那么只能说明：$k|f(i)$
我们可以取$t$个$k$($num>1$)，使这些$k$互质来保证答案的正确性（及只有$[k_1,k_2,..k_t]|f(i)$且$f(i)\not=0$时，才会发生误判）
复杂度：$O(nmt)$

100’

续70’算法。
我们发现
当$f(i)\ mod\ t=0$时，$f(i+pt)\ mod\ t=0$（$p \in N^+$）
当$f(i)\ mod\ t\not=0$时，$f(i+pt)\ mod\ t\not=0$（$p \in N^+$）
所以用类似筛选法筛掉就行了

#define s 10
int prime[s]={4621,8053,9227,3001,2333,5011,7121,8009,9973,1117};
int g[100010][s],p[1000010],n,m;
int f(int x,int y)
{
    int r=0,t=1;
    fr(i,0,n)
    {
        r=(r+t*g[i][y])%prime[y];
        t=t*x%prime[y];
    }
    rt r;
}
int main(){
    n=read();
    m=read();
    fr(i,0,n)
    {
        char c=gc;
        while(c<'0'&&c>'9'&&c!='-')
            c=gc;
        int k=c=='-'?-1:1;
        c=c=='-'?gc:c;
        while(c>='0'&&c<='9')
        {
            fr(j,0,s-1)
                g[i][j]=(g[i][j]*10+c-48)%prime[j];
            c=gc;
        }
        fr(j,0,s-1)
            g[i][j]*=k;
    }
    fr(i,0,s-1)
        fr(j,0,prime[i]-1)
            if(f(j,i))
                for(int k=j;k<=m;k+=prime[i])
                    p[k]=1;
    int ans=0;
    fr(i,1,m)
        if(!p[i])
            ans++;
    printf("%d\n",ans);
    fr(i,1,m)
        if(!p[i])
            printf("%d\n",i);
    rt 0;
}

复杂度：$O(ntk_{max})$