记录编号 139880 评测结果 AAAAAAAAAA
题目名称 [NOIP 2014]解方程 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 1.503 s
提交时间 2014-11-17 14:35:17 内存使用 0.34 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<iomanip>
#include<cstdlib>
#include<vector>
#include<queue>
using namespace std;
typedef long long LL;
const int SIZEN=110,SIZEL=10010;
LL P[20]={34261,999999937,999999929,999999893,999999883,999999797,999999761,999999757,999999751,999999739,999999733,999999677};
int N,M,K=10;
LL A[SIZEN][20]={0};
LL calc(LL x,int k){//f(x) mod P[k]
	x%=P[k];
	LL nowx=1,ans=0,tmp;
	for(int i=0;i<=N;i++){
		tmp=(nowx*A[i][k])%P[k];
		ans=(ans+tmp)%P[k];
		nowx=(nowx*x)%P[k];
	}
	return ans;
}
bool test(LL x){
	for(int i=1;i<=K;i++){
		if(calc(x,i)!=0) return false;
	}
	return true;
}
vector<int> P0_ans;
vector<int> ans;
void work(void){
	for(int i=1;i<=P[0];i++){
		if(calc(i,0)==0) P0_ans.push_back(i);
	}
	for(int i=0;i<P0_ans.size();i++){
		for(int x=P0_ans[i];x<=M;x+=P[0]){
			if(1<=x&&x<=M&&test(x)){
				ans.push_back(x);
				if(ans.size()==N) goto OUT;
			}
		}
	}
	OUT:;
	sort(ans.begin(),ans.end());
	printf("%d\n",ans.size());
	for(int i=0;i<ans.size();i++) printf("%d\n",ans[i]);
}
LL calc_str(char s[],int n,LL P){
	LL sgn=1,now=0;
	for(int i=0;i<n;i++){
		if(s[i]=='-') sgn=-1;
		else{
			now=now*10+s[i]-'0';
			now%=P;
		}
	}
	if(sgn==-1) now=(P-now)%P;
	return now;
}
char str[SIZEL];
void read(void){
	scanf("%d %d\n",&N,&M);
	for(int i=0;i<=N;i++){
		scanf("%s",str);
		int n=strlen(str);
		for(int j=0;j<=K;j++){
			A[i][j]=calc_str(str,n,P[j]);
		}
	}
}
int main(){
	freopen("equationa.in","r",stdin);
	freopen("equationa.out","w",stdout);
	read();
	work();
	return 0;
}