jzoj3935【NOIP2014day2官方数据】解方程（高精度取模）

1. 【NOIP2014day2官方数据】解方程 (Standard IO)

Description

Input

Output

Sample Input
输入1：
2 10
1
-2
1
输入2：
2 10
2
-3
1
输入3：
2 10
1
3
2

Sample Output

输出1：
1
1
输出2：
2
1
2
输出3：
0

Data Constraint

分析：显然，把方程两边同时对一个质数取模，方程解不变，然后高精度对每个a取模即可，考虑到有偶然性，所以最好用两个模数，但是直接枚举答案会超时，根据各种定理，只需枚举到最小那个模数，假设x能使方程成立，那么只要判断x+p，x+2p，x+3p……成立即可。

代码

#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#define N 200
#define mo1 47895201
#define mo2 100007
#define ll long long
#pragma GCC optimize(3)
using namespace std;

char s[20000];
ll a[N][3];
int b[N],n,m,ans;

void get(int x)
{
    int i=0;
    int len=strlen(s);
    bool fl=false;
    if (s[0]=='-')
    {
        fl=true;
        i=1;
    }
    for (int j=i;j<len;j++)
    {
        a[x][1]=(a[x][1]*10LL%mo1+(s[j]-'0'))%mo1;
        a[x][2]=(a[x][2]*10LL%mo2+(s[j]-'0'))%mo2;
    }
    if (fl)
        {
            a[x][1]=mo1-a[x][1];
            a[x][2]=mo2-a[x][2];
        }
}

bool check(int x,int j,int mo)
{
    ll sum=0;
    ll bin=1;
    for (int i=0;i<=n;i++)
    {
        sum=(sum+1LL*a[i][j]*bin%mo)%mo;
        bin=1LL*bin*x%mo;
    }
    if (sum==0) return true;
    return false;
} 

int main()
{
//  freopen("equation8.in","r",stdin);
    scanf("%d%d",&n,&m);
    for (int i=0;i<=n;i++)
    {
        scanf("%s",s);
        get(i);
    }
    for (int i=1;i<=mo2;i++)
        if (check(i,2,mo2)) for (int j=i;j<=m;j+=mo2) if (check(j,1,mo1)) b[++ans]=j;
    printf("%d\n",ans);
    sort(b+1,b+ans+1);
    for (int i=1;i<=ans;i++)
    {
        printf("%d\n",b[i]);
    }
}