ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

「HNOI 2013」数列

2019-02-18 15:39:58  阅读:245  来源: 互联网

标签:数列 read HNOI 个数 int 数组 sum 2013 式子


题目链接

戳我

\(Solution\)

这道题貌似并不难的样子\(QAQ\)

我们发现这个因为有首项的关系所以有点不太好弄.所以我们要将这个首项对答案的影响给去掉.

我们可以构建一个差分数组,我们令他等于\(a[1],a[2]...a[k-1]\)

则一个查分数组对答案的贡献为:

\[\sum_{i=1}^{k-1}n-a[i]\]

然后我们一共有\(m^(k-1)\)个这样的查分数组,所以总贡献为:

\[\sum_{j=1}^{m^{k-1}}\sum_{i=1}^{k-1}n-a[j][i]\]

我们将\(n\)提出来,式子变为:

\[n*m^{k-1}-\sum_{j=1}^{m^{k-1}}\sum_{i=1}^{k-1}a[j][i]\]

所以现在只需要化简后面的式子了.

枚举一些数发现(实际上是我不会证明)
发现在区间\([1,m]的数每个数出现的个数相同\)
至于怎么发现的,打表找规律啊.

这样的话,每个数出现的次数就可以确定了:
\(m^{k-1}\)个数组,每个数组\((k-1)\)个数,
则每个数的个数为:
\[m^{k-1}*(k-1)/m\]
\[=m^{k-2}*(k-1)\]
然后后面式子的值就只需要用这个数乘上\(1+2+3+...+m的值了\)

所以后面式子实际上就是:
\[m^{k-2}*(k-1)*((1+m)*m)/2\]

所以最终答案为:
\[n*m^{k-1}-m^{k-2}*(k-1)*((1+m)*m)/2\]

注意取模的问题啊,好坑!!!

\(Code\)

#include<bits/stdc++.h>
using namespace std;
#define int long long
int read(){
    int x=0,f=1;char c=getchar();
    while(c<'0'||c>'9') f=(c=='-')?-1:1,c=getchar();
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*f;
}
int ksm(int a,int b,int mod){
    int ans=1;
    while(b){
    if(b&1)
        ans=a*ans%mod;
    a=a*a%mod;
    b>>=1;
    }
    return ans%mod;
}
main(){
    int n=read(),k=read(),m=read(),p=read();
    printf("%lld",(n%p*ksm(m,k-1,p)%p-ksm(m,k-2,p)*(k-1)%p*((1+m)*m/2%p)%p+p)%p);
}

标签:数列,read,HNOI,个数,int,数组,sum,2013,式子
来源: https://www.cnblogs.com/hbxblog/p/10395890.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有