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

bzoj千题计划293:bzoj3142: [Hnoi2013]数列

时间:2018-03-20 15:19:53      阅读:103      评论:0      收藏:0      [点我收藏+]

标签:col   lld   href   相同   sig   分数   blog   pac   print   

http://www.lydsy.com/JudgeOnline/problem.php?id=3142

 

如果已知数列的差分数列a[1]~a[k-1]

那么这种差分方式对答案的贡献为 N-Σ a[i],i∈[1,k-1]

差分数列一共有多少种? M^(k-1) 种

所以ans=Σ  (N-Σa[i]) = M^(k-1) * N - Σ Σ a[i] = M^(k-1) *N-(k-1)*M^(k-1) /M * (M+1)*M/2

后面是因为一共 M个数 出现(k-1)*M^(k-1) 次,每个数出现的次数相同

所以每个数出现 (k-1)*M^(k-1) / M 次 

(M+1)*M/2 是所有的数各出现一次的和

 

没做出来的原因:没有展开Σ,没有考虑每个数的出现次数相同

 

#include<cstdio>
 
using namespace std;
 
typedef long long LL;
 
LL Pow(LL a,LL b,LL p)
{
    LL res=1;
    for(;b;a=a*a%p,b>>=1)
        if(b&1) res=res*a%p;
    return res;
}
 
int main()
{
    LL n,k,m,p;
    scanf("%lld%lld%lld%lld",&n,&k,&m,&p);
    printf("%lld",(n%p*Pow(m,k-1,p)%p-(k-1)*Pow(m,k-2,p)%p*(m*(m+1)/2%p)%p+p)%p);
}

 

bzoj千题计划293:bzoj3142: [Hnoi2013]数列

标签:col   lld   href   相同   sig   分数   blog   pac   print   

原文地址:https://www.cnblogs.com/TheRoadToTheGold/p/8608933.html

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