【NOIP2014day2官方数据】解方程

传送门

https://jzoj.net/senior/#main/show/3935

题目大意

给你a0..an，求x属于[1..m]且a0+a1*x+a2*x^2+a3*x^3+……+an*x^n的所有数。
注意：

思路

因为数据太大，so要用特殊读入。
对于这个大数据，我们可以给每个值取模（多mod几个），读入取mod即可
然后我们吧公式稍微转换一下，不断提取公因数 得:
a1+x(a2+x(a3+x(a4+………..+x(an)))))))))…..
所以我们可以递推：ans=ans*i+a[i]。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,m,a[111],b[11111];
long long ans=0,mo=100000007;
//__attribute__((optimize("-O3")))
long long read()
{
    long long zf=1,sum=0;
    char c=getchar();
    while (c<'0'||c>'9') 
        if (c=='-') 
        {
            zf=-1;
            c=getchar();
        }
    while (c>='0'&&c<='9') sum=sum*10+c-'0',sum%=mo,c=getchar();
    return sum*zf;
}
long long check(long long i)
{
    long long sum=a[n];
    for(long long j=n-1; j>=0; j--)
    {
        sum=(sum*i+a[j])%mo;
    }
    return sum;
}
int main()
{
//  freopen("data.in","r",stdin);
//  freopen("data.out","w",stdout);
//  printf("done\n");
    n=read();m=read();
    for(long long i=0; i<=n; i++) a[i]=read();
    for(long long i=1; i<=m; i++) if(!check(i)) b[++ans]=i;
    printf("%lld\n",ans);
    for(long long i=1; i<=ans; i++) printf("%lld\n",b[i]);
    fclose(stdin); fclose(stdout);
}