luogu 2312 解方程 乱搞+取模

思路非常好想,但是你很难想到去用这个算法,因为这个几乎就是个乱搞~

我们发现多项式中每一个系数都很大,但是 $m$ 却很小,即最多只用 $10^6$ 个整数需要验证.

我们知道,如果一个数等于 $0$,那么这个数模任何一个数也都应该该等于 $0$

所以可以直接取 $3$ 个左右的质数当模数,分别带值,取模,然后判一下等不等于 $0$.

当然,带值的部分可以用秦九昭算法,但是我感觉这只算是常数上的优化吧~

只能在 luogu 上过,bz 上过不去~

复杂度 $O(n\times m)$

#include <bits/stdc++.h>        
#define N 1000005    
#define LL long long 
#define setIO(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)      
using namespace std;         
// 0 ~ 15                 
const LL mod[]={998244353,19,1e9+7, 1e9+9, 233,233233,23,17,19,11,613317,119,911,2332,2323,1415,1717};    
int n,m;       
char str[103][N];       
vector<int>v;     
LL a[20][N];    
int check(int tmp) 
{      
	// 15 个模数    
	int i,j;        
	for(i=0;i<=0;++i) 
	{             
		LL temp=0ll;   
		a[i][n+1]=0ll;  
		for(j=n+1;j>=0;--j)           
		    temp=(temp*1ll*tmp%mod[i]+a[i][j])%mod[i];                 
		if(temp!=0) return 0; 
	}
	return 1;   
}
inline void Init() 
{ 
	int i,j; 
	for(i=0;i<=0;++i) 
	{    
		for(j=0;j<=n;++j)        
		{ 
			LL tmp=0,base=1ll;   
			int len=strlen(str[j]);    
			for(int k=len-1;k>=0;--k) 
			{    
				if(k==0&&str[j][k]=='-') 
				{ 
					tmp=(mod[i]-tmp%mod[i])%mod[i]; 
				}
				else 
				{
					tmp=(tmp+(str[j][k]-'0')*base)%mod[i], base=base*10%mod[i]; 
				}
			}
			a[i][j]=tmp;    
			// printf("%d %lld\n",j,a[j]); 
		}   
	}
}
int main() 
{             
	// setIO("input"); 
	int i,j;
	scanf("%d%d",&n,&m);    
	for(i=0;i<=n;++i)  scanf("%s",str[i]);      
	Init();               
	for(i=1;i<=m;++i)  if(check(i))   v.push_back(i); 
	printf("%d\n",v.size()); 
	for(i=0;i<v.size();++i)   printf("%d\n",v[i]);     
	return 0;
}

  

posted @ 2019-11-05 20:32  EM-LGH  阅读(134)  评论(0编辑  收藏  举报