[FROM LUOGU]解方程

传送门

SOL
首先有很大的数据规模,所以只能取模判断是否为0
其次计算这个多项式的值可能比较复杂,如果用pow函数或者快速幂复杂度就肯定炸了,常数小的话可以过70分
所以可以迭代:
计 算 第 一 次 : a n ∗ x + a n − 1 计算第一次:an*x+an-1 anx+an1
计 算 第 二 次 : ( a n ∗ x + a n − 1 ) ∗ x + a n − 2 计算第二次:(an*x+an-1)*x+an-2 (anx+an1)x+an2
… … ……
如此以来,求一个多项式就只需计算 n n n
除此之外,可以进行一些剪枝,例如对于一个小质数 p p p,先求 1 < = x < p 1<=x<p 1<=x<p的答案,然后对于 x < = p < = m x<=p<=m x<=p<=m,看一下 x   m o d   p x~mod~p x mod p即可
为了保险起见,最好多用几个质数(就像Hash一样),如果取模出来的答案都为 0 0 0,才确定为解,这样同时也可以减去一些无用状态

代码:

#include<bits/stdc++.h>
using namespace std;
#define re register
#define ll long long
inline char nc(){
    static char buf[10000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,10000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd(){
	int re data=0;static char ch=0;ch=nc();
	while(!isdigit(ch))ch=nc();
	while(isdigit(ch))data=(data<<1)+(data<<3)+(ch^48),ch=nc();
	return data;
}
const int N=1e6+5,pa=1e5+7,pb=1e9+7;
int n,m,a[105],b[N],ans[N],cnt;char ch;bool w,tag[pa];ll aa,bb;
inline int f(ll x,int mod,int *t,ll sum=0){sum=t[n];for(int re i=n-1;~i;--i)sum=(sum*x+t[i])%mod;return sum;}
signed main(){
	n=rd(),m=rd();
	for(int re i=0;i<=n;++i){
		aa=bb=0,w=0,ch=nc();
		while(!isdigit(ch)&&ch!='-')ch=nc();
		if(ch=='-')w=1,ch=nc();
		while(isdigit(ch))aa=((aa<<1)+(aa<<3)+(ch^48))%pa,bb=((bb<<1)+(bb<<3)+(ch^48))%pb,ch=nc();
		a[i]=w?pa-aa:aa,b[i]=w?pb-bb:bb;
	}
	for(int re i=0;i^pa;++i)if(!f(i,pa,a))tag[i]=1;
	for(int re i=1;i<=m;++i){
		if(!tag[i%pa])continue;
		if(!f(i,pb,b))ans[++cnt]=i;
	}printf("%d\n",cnt);
	for(int re i=1;i<=cnt;++i)printf("%d\n",ans[i]);
	exit(0);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值