【BZOJ】3751: [NOIP2014]解方程

题意

求(sum_{i=0}^{n} a_i x^i = 0)在([1, m])内的整数解。（(0 < n le 100, |a_i| le 10^{10000}, a_n eq 0, m le 1000000)）

分析

神题。

题解

我们可以取几个质数然后对应取模来计算即可。可是在经过变态的加强数据后，不是tle就是wa。
于是我们可以用一个正确率很低的优化。
令$f(x, p) = sum_{i=0}^{n} a_i x^i pmod{p} (，容易知道)f(x, p) = f(x mod p, p)(，所以我们取一些较小的)p(，暴力预处理一下然后)O(1)$就可以判断了。
最后为了提高正确率，我们再将验证通过的数再用大质数验证一次。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int p[4]={11261, 14843, 21893, 9851197}, f[25005][3], n, m, a[105][4], ans[1000015];
int main() {
	scanf("%d%d", &n, &m);
	for(int i=0; i<=n; ++i) {
		char c=getchar();
		int f=1;
		for(; c<'0'||c>'9'; c=getchar()) {
			if(c=='-') {
				f=-1;
			}
		}
		for(; c>='0'&&c<='9'; c=getchar()) {
			for(int k=0; k<4; ++k) {
				a[i][k]*=10;
				a[i][k]+=c-'0';
				a[i][k]%=p[k];
			}
		}
		if(f==-1) {
			for(int k=0; k<4; ++k) {
				a[i][k]=p[k]-a[i][k];
			}
		}
	}
	for(int k=0; k<3; ++k) {
		for(int x=1; x<p[k]; ++x) {
			for(int i=0, now=1; i<=n; ++i, now=now*x%p[k]) {
				f[x][k]+=now*a[i][k]%p[k];
				if(f[x][k]>=p[k]) {
					f[x][k]-=p[k];
				}
			}
		}
	}
	int cnt=0;
	for(int i=1; i<=m; ++i) {
		bool flag=1;
		for(int k=0; k<3; ++k) {
			if(f[i%p[k]][k]) {
				flag=0;
				break;
			}
		}
		if(flag) {
			int ret=0;
			for(int j=0, now=1; j<=n; ++j, now=(ll)now*i%p[3]) {
				ret+=(ll)now*a[j][3]%p[3];
				if(ret>=p[3]) {
					ret-=p[3];
				}
			}
			if(ret==0) {
				ans[cnt++]=i;
			}
		}
	}
	printf("%d
", cnt);
	for(int i=0; i<cnt; ++i) {
		printf("%d
", ans[i]);
	}
	return 0;
}