2018.11.02【NOIP2014】【洛谷P2312】解方程（乱搞）（哈希）

传送门

---

解析：

真的是一道纯粹的乱搞题。

首先，高次方程式没有求根公式的，所以我们唯一的方法就剩下了暴力枚举。

高达$1e1000$的系数和只有$1e6$的取值范围让我明白了这道题的乱搞本质。

计算直接用秦九韶公式来搞一波。
随便取几个稍微大一点的质数哈希（我选的是10000以下最大的四个质数）。选太小会被卡，选少了也会被卡。

然后处理出所有数在$mod$意义下的最终结果，就可以做到$O(1)$验证取值范围内的每个数了。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

cs int prime[4]={9973,9967,9949,9941};
int num[105][4];

inline void getbignum(int pos){
	re char c;
	re bool f=0;
	while(!isdigit(c=gc()))if(c=='-')f=1;
	for(int re i=0;i<4;++i)num[pos][i]=c^48;
	while(isdigit(c=gc()))for(int re i=0;i<4;++i)
	num[pos][i]=((num[pos][i]<<1)+(num[pos][i]<<3)+(c^48))%prime[i];
	if(f)for(int re i=0;i<4;++i)num[pos][i]=prime[i]-num[pos][i];
}

int n,m;
int res[10000][4];
inline void init(){
	for(int re p=0;p<4;++p){
		for(int re k=0;k<prime[p];++k){
			int sum=0;
			for(int re i=n;~i;--i)sum=(sum*k+num[i][p])%prime[p];
			res[k][p]=sum;
		}
	}
}

int a[1000006],cnt;
signed main(){
	scanf("%d%d",&n,&m);
	for(int re i=0;i<=n;++i)getbignum(i);
	init();
	for(int re i=1;i<=m;++i){
		bool flag=true;
		for(int re j=0;j<4;++j)if(res[i%prime[j]][j])flag=false;
		if(flag)a[++cnt]=i;
	}
	printf("%d\n",cnt);
	for(int re i=1;i<=cnt;++i){
		printf("%d\n",a[i]);
	}
	return 0;
}

转载于:https://www.cnblogs.com/zxyoi/p/10047119.html