码迷,mamicode.com
首页 > 其他好文 > 详细

[NOIP2014提高组]解方程

时间:2017-09-17 19:05:50      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:正整数   style   bsp   高精度   解题思路   const   src   std   为什么   

题目:BZOJ3751、洛谷P2312、UOJ#20、Vijos P1910、codevs3732。

题目大意:已知多项式方程:

技术分享

求这个方程在[1, m]内的整数解(n 和 m 均为正整数)。

解题思路:因为$0=0$(废话),能得出$0+x·p\equiv 0(mod\ p)$。

也就是当方程右边为0时,方程左边mod p为0。

但方程左边mod p等于0时,方程右边不一定等于0。

但是也不一定不等于0。

所以我们如果多引入几个p(最好是素数),对其进行测试,发现都为0的话,那我们就可以认为它就是一个方程的解。

在读入a时,每次读进来一位就对p取模,即可避免高精度。

然后设f[a][b]表示在模第a个模数意义下,x为b时方程的解,求出f即可(由于x和x+p在方程中答案一样,所以x不需要枚举到m,只需枚举到模数即可)。然后对于$[1,m]$的每个整数x,当f[0][x]=f[1][x]=...=f[k-1][x]=0(k为你选的模数个数)时,记录答案。

选模数时选4~5个10000~30000的素数基本就可以了。BZOJ的数据比较强,模数选的不好会错(搞不懂为什么是Output  Limit  Exceed),一般的话会发现耗时特别长(其他1200+ms,BZOj9400+ms)。

时间复杂度$O(nkp)$。

C++ Code:

#include<cstdio>
#include<cctype>
using namespace std;
const int P[]={11261,19433,21221,25537,14843};
int n,m,a[5][155],mi[33333],f[5][33333],ans[1000005],cnt;
bool check(int x){
	for(int t=0;t<5;++t){
		if(f[t][x%P[t]])return 0;
	}
	return 1;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=0;i<=n;++i){
		bool flag=false;
		char c=getchar();
		for(;!isdigit(c);c=getchar())flag=c==‘-‘;
		for(;isdigit(c);c=getchar()){
			for(int t=0;t<5;++t){
				a[t][i]=(a[t][i]*10+(c^‘0‘))%P[t];
			}
		}
		if(flag)
		for(int t=0;t<5;++t)a[t][i]=-a[t][i];
	}
	for(int t=0;t<5;++t){
		for(int x=1;x<=P[t];++x){
			mi[0]=1;
			for(int i=1;i<=n;++i)mi[i]=mi[i-1]*x%P[t];
			int s=0;
			for(int i=0;i<=n;++i)
			s=(s+a[t][i]*mi[i])%P[t];
			if(s<0)s+=P[t];
			f[t][x]=s;
		}
	}
	cnt=0;
	for(int i=1;i<=m;++i)
	if(check(i))ans[++cnt]=i;
	printf("%d\n",cnt);
	for(int i=1;i<=cnt;++i)printf("%d\n",ans[i]);
	return 0;
}

 

  

 

[NOIP2014提高组]解方程

标签:正整数   style   bsp   高精度   解题思路   const   src   std   为什么   

原文地址:http://www.cnblogs.com/Mrsrz/p/7536314.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!