HNOI(数列)

本文深入探讨了一种数学推导方法,通过定义差分数组并利用数学公式简化复杂问题,实现了算法的优化。文章详细解释了如何通过枚举和求和操作减少计算复杂度,给出了一段C++代码实现,展示了理论到实践的转化过程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

这题没有信息思想,其实就是数学推式子
——————————————————————————————————————
定义差分数组ai定义差分数组a_iai
∑a1=1m∑a2=1m.......∑ak−1=1m(n−∑i=1k−1)=mk−1n−∑a1=1m......∑ak−1=1m∑i=1k−1a[i]\sum _{{a_1}=1}^m\sum _{a_2=1}^m .......\sum _{a_{k-1}=1}^m(n-\sum _{i=1}^{k-1})=m^{k-1}n-\sum_{a_{1}=1}^m......\sum_{a_{k-1}=1}^{m}\sum_{i=1}^{k-1}a[i]a1=1ma2=1m.......ak1=1m(ni=1k1)=mk1na1=1m......ak1=1mi=1k1a[i]

因为后一个式子∑a1=1m......∑ak−1=1m∑i=1k−1a[i]\sum_{a_{1}=1}^m......\sum_{a_{k-1}=1}^{m}\sum_{i=1}^{k-1}a[i]a1=1m......ak1=1mi=1k1a[i]的意义是所有可能的aia_iai值求和,总共有(k−1)∗mk−1个数贡献答案(k-1)*m^{k-1}个数贡献答案k1mk1,因为枚举每个aia_iai的次数相同,所以这个式子可以装化为(k−1)∗mk−1∗∑i=1mi(k-1)*m^{k-1}* \sum_{i=1}^{m}i(k1)mk1i=1mi


#include<iostream>
#include<cstdio>

using namespace std;
typedef long long ll;
ll mod,n,k,m;
ll ksm(ll x,ll pow){
	if(pow==0) return 1; 
	ll ans=1%mod,res=x%mod;
	while(pow){
		if(pow&1) ans=ans*res%mod;
		res=res*res%mod;
	    pow>>=1;
	}
	return ans%mod;
}
int main(){
	scanf("%lld%lld%lld%lld",&n,&k,&m,&mod);
	if(k==1){
		printf("%lld",n%mod);
		return 0;
	}
	n%=mod;// very important 不取余会爆精度
	ll tmp=ksm(m,k-2);
	ll ans=((m*n%mod-(k-1)%mod*((((m+1)*m)/2)%mod))%mod)%mod;
	if(ans<0) ans+=mod;
	printf("%lld",ans*tmp%mod); 
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值