洛谷 P2312 解方程
JDOJ 2890: [NOIP2014]解方程 D2 T3
Description
已知多项式方程:
a0+a1x+a2x2+⋯+anxn=0
求这个方程在[1, m]内的整数解( n和 m均为正整数)。
Input
输入共 n+2行。
第一行包含 2个整数 n、m,每两个整数之间用一个空格隔开。
接下来的 n+1行每行包含一个整数,依次为 a0,a1,a2,…,an。
Output
第一行输出方程在 [1, m]内的整数解的个数。
接下来每行一个整数,按照从小到大的顺序依次输出方程在 [1, m]内的一个整数解。
Sample Input
Input I: 2 10 1 -2 1 Input II: 2 10 2 -3 1 Input III: 2 10 1 3 2
Sample Output
Output I: 1 1 Output II: 2 1 2 Output III: 0
HINT
【数据说明】
对于 30%的数据,0< n ≤2, |ai|≤100,an≠0, m≤100;
对于 50%的数据,0< n ≤100,|ai| ≤10100,an≠0,m≤100;
对于 70%的数据,0< n ≤100,|ai| ≤1010000,an≠0,m≤10000;
对于 100%的数据,0< n ≤100,|ai| ≤1010000,an≠0,m≤1000000。
Source
题解:
高一重点实验班同学(比如\(JDFZ\)高一6班)和高二同学应该很熟悉这个式子(滑稽.jpg)
这就是秦九韶算法的裸题嘛!
(为了照顾高一非重点实验班同学和高一以下同学,友情提供秦九韶算法的讲解:)
但是真的这么简单么?
注意数据范围,这一看就是高精度的题。但是可以不用高精做(高精要狗命)可以在读入和计算的时候不断取模。
当\(f(x)=0\),那么肯定会有\(f(x)mod \,\,\,p=0\)。\(p\)最好是质数。
这种方法很玄学,正确性也无法保证(多试试就好了)。
代码:
#include<cstdio> #define int long long using namespace std; const int maxn=110; const int mod=1e9+7; int a[maxn],ans,n,m,b[1000010],tot; bool flag; char *p1,*p2,buf[100000]; #define nc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++) inline int read() { int x=0,f=1; char ch=nc(); while(ch<48){if(ch=='-')f=-1;ch=nc();} while(ch>47) x=((((x<<2)+x)<<1)+ch-48)%mod,ch=nc(); return x*f; } inline bool check(int x) { int ret=a[n]; for(register int i=n-1;i>=0;i--) ret=(ret*x+a[i])%mod; return ret?0:1; } inline void print(int x) { if(x<0) putchar('-'),x=-x; if(x>9) print(x/10); putchar(x%10+48); } signed main() { n=read(),m=read(); for(register int i=0;i<=n;i++) a[i]=read(); for(register int i=1;i<=m;i++) if(check(i)) { flag=1; ans++; b[++tot]=i; } if(!flag) { print(0); return 0; } else { print(ans); puts(""); for(register int i=1;i<=tot;i++) print(b[i]),puts(""); } return 0; }