- 【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]);
}
}