记录编号 |
139880 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[NOIP 2014]解方程 |
最终得分 |
100 |
用户昵称 |
cstdio |
是否通过 |
通过 |
代码语言 |
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;
}