标签:le return LL cases define 杂题 mod
题意
给定\(n,m\),初始序列\(\{a_i\}\)全为\(0\),可以进行任意次操作,选取一段长度为\(m\)的区间,依次赋值为\(1,2,\cdots ,m\)
求能得到多少种任意位置非\(0\)的序列
\(n,m\le 10^6\)
做法
感谢神仙MAOoo的耐心教导
令\(f_i(j)\)为第\(i\)个位置,填\(j\),前\(i\)个位置不同的方案数
考虑若目前所处的位置\(\le n-m+1\),那么这个位置可以赋值为任意\([1,m]\)的数:
对于\(i\le n-m+1\)
观察到\(f_i(j)-f_{i-1}(j-1)=f_{i-1}(m)(j\neq 1)\)
那么就是对于dp数组,滑动窗口,维护一个差值
仅需求出\(f_i(1)\)
对于\(i>n-m+1\)类似
这道题很难讲得很清楚,所以放份代码
Code
#include<bits/stdc++.h>
typedef int LL;
typedef double dl;
#define opt operator
#define pb push_back
#define pii std::pair<LL,LL>
const LL maxn=1e6+9,mod=1e9+7,inf=0x3f3f3f3f;
LL Read(){
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<<3ll)+(x<<1ll)+c-'0'; c=getchar();
}return x*f;
}
void Chkmin(LL &x,LL y){
if(y<x) x=y;
}
void Chkmax(LL &x,LL y){
if(y>x) x=y;
}
LL add(LL x,LL y){
return x+=y,x>=mod?x-mod:x;
}
LL dec(LL x,LL y){
return x-=y,x<0?x+mod:x;
}
LL mul(LL x,LL y){
return 1ll*x*y%mod;
}
LL Pow(LL base,LL b){
LL ret(1); while(b){
if(b&1) ret=mul(ret,base); base=mul(base,base); b>>=1;
}return ret;
}
LL n,m;
LL dp[maxn<<1];
int main(){
LL n(Read()),m(Read());
LL l(1),r(m);
dp[r]=1;
LL delta(0),sum(1);
++l; ++r;
for(LL i=1;i<n;++i,++l,++r){
if(i+1<=n-m+1){
sum=dec(sum,dp[l-1]);
dp[r]=add(1ll*(m-2)*delta%mod,sum);
}
delta=add(delta,add(delta,dp[l-1]));
sum=add(sum,dp[r]);
}
printf("%d\n",add(dp[l-1],delta));
return 0;
}
标签:le,return,LL,cases,define,杂题,mod 来源: https://www.cnblogs.com/Grice/p/14284975.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。