把博客园图标替换成自己的图标
把博客园图标替换成自己的图标end

【规律】ABC179 E Sequence Sum

题目链接

题目解析

刚开始把题目看错了,看成了乘法,然后直接上指数循环节了qwq

虽然N很大,但是M却很小,Ai的范围被限定在[0,M1],而这个f()是类似于函数一样的东西,一旦第一个数给定,后面的数就唯一确定(拉普拉斯信条。而根据抽屉原理,至少m+1个数之后就会重复,那么重复之后就会产生循环节!那么就做完了!


►Code View

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define LL long long
#define N 100005
#define DEL 100000
#define INF 0x3f3f3f3f
LL rd()
{
	LL x=0,f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-')f=-1; c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48); c=getchar();}
	return f*x;
}
LL n,x,m,ans;
LL a[N<<1],s[N<<1];
int pos[N];
int main()
{
	n=rd(),x=rd(),m=rd();
	a[1]=x,s[1]=x;
	for(int i=2;i<=min(100000ll,n);i++)
	{
		a[i]=a[i-1]*a[i-1]%m;
		if(pos[a[i]])
		{
			int l=pos[a[i]];
			LL len=i-l,sum=s[i-1]-s[l-1];
			ans=s[l-1];
			LL tmp=n-(l-1);
			LL p=tmp/len,q=tmp%len;
			ans=ans+p*sum;
			ans=ans+s[l+q-1]-s[l-1];
			printf("%lld\n",ans);
			return 0;
		}
		s[i]=s[i-1]+a[i];
		pos[a[i]]=i;
	}
	printf("%lld\n",s[n]);
	return 0;
} 
/*LL ksm(LL a,LL b,LL MOD)
{
	LL res=1ll;
	while(b)
	{
		if(b&1) res=res*a%MOD;
		a=a*a%MOD;
		b>>=1;
	}
	return res;
}
//phi(x)=x*(1-1/p1)*(1-1/p2)*...*(1-1/pn)
LL Phi(LL a)
{
	LL res=a;
	for(LL i=2;i*i<=a;i++)
		if(a%i==0)
		{
			res=res/i*(i-1);
			while(a%i==0) a/=i;
		}
	if(a>1) res=res/a*(a-1);
	return res;
}
//x^{2^n-1}%MOD=x^{(2^n-1)%phi+phi}%MOD
int main()
{
	n=rd(),x=rd(),m=rd();
	LL phi=Phi(m);
	LL ind=(ksm(2,n,phi)-1+phi)%phi+phi;
	LL ans=ksm(x,ind,m);
	printf("%lld\n",ans);
	return 0;
}*///想成乘法了 qwq 

posted @   Starlight_Glimmer  阅读(152)  评论(0编辑  收藏  举报
编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 使用C#创建一个MCP客户端
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
历史上的今天:
2019-11-15 CSP考试策略
浏览器标题切换
浏览器标题切换end
点击右上角即可分享
微信分享提示