Noip 提高组 2014 Day2 T3 解方程 同余系
2017-10-27 13:24
246 查看
Description
已知多项式方程:
求这个方程在[1,m]内的整数解(n和m均为正整数)。
Input
第一行包含2个整数n、m,每两个整数之间用一个空格隔开。
接下来的n+1行每行包含一个整数,依次为a0,a1,a2,…,an。
Output
第一行输出方程在[1,m]内的整数解的个数。
接下来每行一个整数,按照从小到大的顺序依次输出方程在[1,m]内的一个整数解。
Sample Input
2 10
2
-3
1
Sample Output
2
1
2
HINT
对于100%的数据,0 < n<=100,|ai|≤1010000,an≠0,m≤1000000。
传送门
……好可恶早上考试写这题结果写渣了= =
考虑一下如何来解题。
首先当然想到暴力枚举一个1~m的x了,然后代入方程检验。
然而……m很大,而且ai还是高精度+负数所以= =
那么可以想到把ai化简,使得不用高精度去做,
很容易想到哈希了,取几个模数pi,分别求出x代入方程%pi的值,
如果都是0,那就近似认为x是一个解。
如果取得好、冲突少,就可以A了,所以往这里去想想。
取模数T个,那么这种方法的时间复杂度是O(nmT)
有些人竟然就这么卡常过掉了= =
那不卡常怎么办呢。。
可以发现,假如x代入方程%pi=0,那么(x+k*pi)代入方程%pi肯定也是0
由此发现,假如x和y同余,且y%pi=0,那么x%pi=0,
那就可以简化了:我们只用代入0~(pi-1)的数,并且记录值是否为0,
然后枚举1~m的x时,只要判断x%pi是否为0,就知道x代入方程的结果了。
所谓“同余系”。
因此我们取的质数不能太大(不然就会T了QAQ)
而且为了避免long long(long long慢啊。。),
取的质数在4W左右就好了。。
接过TLE了一年= =气死网上随便找了组模数(看上去差不多啊喂。。)
uoj,luogu……都AC啊(ccccccc)
后面卡常到怀疑人生终于8sA掉bzoj的加强数据……(吐血)
已知多项式方程:
a0+a1*x+a2*x^2+...+an*x^n=0
求这个方程在[1,m]内的整数解(n和m均为正整数)。
Input
第一行包含2个整数n、m,每两个整数之间用一个空格隔开。
接下来的n+1行每行包含一个整数,依次为a0,a1,a2,…,an。
Output
第一行输出方程在[1,m]内的整数解的个数。
接下来每行一个整数,按照从小到大的顺序依次输出方程在[1,m]内的一个整数解。
Sample Input
2 10
2
-3
1
Sample Output
2
1
2
HINT
对于100%的数据,0 < n<=100,|ai|≤1010000,an≠0,m≤1000000。
传送门
……好可恶早上考试写这题结果写渣了= =
考虑一下如何来解题。
首先当然想到暴力枚举一个1~m的x了,然后代入方程检验。
然而……m很大,而且ai还是高精度+负数所以= =
那么可以想到把ai化简,使得不用高精度去做,
很容易想到哈希了,取几个模数pi,分别求出x代入方程%pi的值,
如果都是0,那就近似认为x是一个解。
如果取得好、冲突少,就可以A了,所以往这里去想想。
取模数T个,那么这种方法的时间复杂度是O(nmT)
有些人竟然就这么卡常过掉了= =
那不卡常怎么办呢。。
可以发现,假如x代入方程%pi=0,那么(x+k*pi)代入方程%pi肯定也是0
由此发现,假如x和y同余,且y%pi=0,那么x%pi=0,
那就可以简化了:我们只用代入0~(pi-1)的数,并且记录值是否为0,
然后枚举1~m的x时,只要判断x%pi是否为0,就知道x代入方程的结果了。
所谓“同余系”。
因此我们取的质数不能太大(不然就会T了QAQ)
而且为了避免long long(long long慢啊。。),
取的质数在4W左右就好了。。
接过TLE了一年= =气死网上随便找了组模数(看上去差不多啊喂。。)
uoj,luogu……都AC啊(ccccccc)
后面卡常到怀疑人生终于8sA掉bzoj的加强数据……(吐血)
#include<bits/stdc++.h> #define ll long long using namespace std; const int prime[6]={30011,11261,14843,19997,10007,21893}; int n,m,ANS[1000005],a[105][6]; bool f[30015][6]; char s[10005]; void get(int o){ scanf("%s",s+1); int f=1,sta=1; if (s[1]=='-') f=-1,sta=2; for (int i=sta;s[i];i++) for (int j=0;j<6;j++) a[o][j]=(a[o][j]*10+s[i]-'0')%prime[j]; for (int j=0;j<6;j++) a[o][j]*=f; } bool ok(int x,int mod){ int p=0,t=1; for (int i=0;i<=n;i++){ p=(p+a[i][mod]*t)%prime[mod]; t=t*x%prime[mod]; } return p==0; } int main(){ scanf("%d%d",&n,&m); for (int i=0;i<=n;i++) get(i); for (int i=0;i<6;i++) for (int j=0;j<prime[i];j++) f[j][i]=ok(j,i); int cnt=0; for (int i=1;i<=m;i++){ bool ff=1; for (int j=0;j<6;j++) if (!f[i%prime[j]][j]){ff=0;break;} if (ff) ANS[++cnt]=i; } printf("%d\n",cnt); for (int i=1;i<=cnt;i++) printf("%d\n",ANS[i]); return 0; }
相关文章推荐
- NOIP 2014 提高组 Day2 T1 无线网路发射器选址
- 解方程 NOIP2014 提高组 Day2 T3
- 2014 Noip提高组 Day2
- 【暴力】洛谷 P2038 NOIP2014提高组 day2 T1 无线网络发射器选址
- [解题报告] NOIP 2014 提高组Day2试题
- 【DFS】【最短路】【spfa】【BFS】洛谷P2296 NOIP2014提高组 day2 T2 寻找道路
- 【noip2014提高组】【Day2】【解题报告】
- NOIP2014复赛提高组day2(A:无线网络发射器选址 B:寻找道路 C:解方程)
- NOIP2014提高组Day2
- 总结(【中山市选2014】投票 【中山市选2014】dwarf tower 【中山市选2014】图 【NOIP2013提高组day2】华容道)
- NOIP 2014 提高组 Day2
- NOIP 2014 提高组 Day2 T2 寻找道路
- 无线网络发射选址 NOIP2014 提高组 Day2 T1
- Noip2014 Day2 T3 解方程(数论+取模运算)
- [NOIP2014] 提高组 洛谷P2038 无线网络发射器选址
- 【NOIP】提高组2014
- 【NOIP2016提高组day2】愤怒的小鸟
- NOIP2014提高组 酱油记
- NOIP2014提高组初赛总结(题目+易错点+解析)
- NOIP2014提高组D.发射器选址(解题报告)