记录编号 143758 评测结果 AAAAAAAAAA
题目名称 [NOIP 2014]解方程 最终得分 100
用户昵称 Gravatar天一阁 是否通过 通过
代码语言 C++ 运行时间 1.462 s
提交时间 2014-12-17 17:56:59 内存使用 0.35 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#define LL long long
using namespace std;
LL M[11]={9,42227,100000007,99999931,999999913,937713711,1313131313,1717171717,999999893,133333119,214748337};
LL P[110][12];
int n,m;
char str[30010];
LL calc(LL Mod){
	LL res=0,si=0;
	for(int i=0;i<strlen(str);i++){
		if(str[i]=='-') si=1;
		else (res=10*res+str[i]-'0')%=Mod;
	}
	if(si) res=(Mod-res)%Mod;
	return res%Mod;
}
LL work(LL x,int k){
	x%=M[k];
	LL res=0,Pow=1,tmp;
	for(int i=0;i<=n;i++){
		tmp=(Pow*P[i][k])%M[k];
		res=(res+tmp)%M[k];
		Pow=(Pow*x)%M[k];
	}
	return res;
}
bool pass(LL x){
	for(int i=2;i<=M[0];i++) if(work(x,i)) return false;
	return true;
}
vector<int> S1,ans;
int main(){
	freopen("equationa.in","r",stdin);
	freopen("equationa.out","w",stdout);
	scanf("%d %d",&n,&m);
	for(int i=0;i<=n;i++){
		scanf("%s",str);
		for(int j=1;j<=M[0];j++){
			P[i][j]=calc(M[j]);
			//printf("P = %d ",P[i][j]);
		}
	}
	for(int i=1;i<=M[1];i++)
		if(work(i,1)==0) S1.push_back(i);
	for(int i=0;i<S1.size();i++)
		for(int j=S1[i];j<=m;j+=M[1])
			if(pass(j)){
				ans.push_back(j);
				if(ans.size()==n) goto DONE;
			}
	DONE:;
	sort(ans.begin(),ans.end());
	printf("%d\n",ans.size());
	for(int i=0;i<ans.size();i++) printf("%d\n",ans[i]);
	return 0;
}