[HNOI2011] 卡农

题目传送门

众所周知卡农是一种复调音乐的写作技法,小余在听卡农音乐时灵感大发,发明了一种新的音乐谱写规则。

我还真不知道

解法

做计数DP的题先观察性质;

性质1
由于一段音乐每种音符被奏响的次数为偶数,
那么我们如果确定了前 m − 1 m-1 m1个集合那么第 m m m个集合就随之确定了(注意第 m m m个集合是最后确定的集合)
性质2
集合不可为 ∅ \empty 且对于任意两个选出的集合 A A A B B B,有 A ≠ B A\ne B A=B

设计状态

f i : 前 i 个集合合法可得到的合法方案数 f_i:前i个集合合法可得到的合法方案数 fi:i个集合合法可得到的合法方案数

状态转移

之后我们考虑容斥
显然 f i = C 2 n − 1 m f_i=C_{2^{n}-1}^{m} fi=C2n1m
根据性质1 f i = C 2 n − 1 m − 1 f_i=C_{2^{n}-1}^{m-1} fi=C2n1m1
考虑集合为空的情况:
显然为 f i − 1 f_{i-1} fi1
而后是重复的情况:
f i − 2 ∗ ( 2 n − 1 − ( i − 2 ) ) f_{i-2}*(2^n-1-(i-2)) fi2(2n1(i2))
再仔细想想,每种方案会被重算i次,比如:
{ { 1 , 2 } , { 2 , 3 } } \{\{1,2\},\{2,3\}\} {{1,2},{2,3}}可以由 { 1 , 2 } \{1,2\} {1,2} { 2 , 3 } \{2,3\} {2,3}转移而来
最后:
f i = C 2 n − 1 m − 1 − f i − 1 − f i − 2 ∗ ( 2 n − 1 − ( i − 2 ) ) i f_i=\frac{C_{2^{n}-1}^{m-1}-f_{i-1}-f_{i-2}*(2^n-1-(i-2))}{i} fi=iC2n1m1fi1fi2(2n1(i2))

code:
#include<bits/stdc++.h>
using namespace std;
template<class T>inline void read(T &x){
	x=0;
	char c=getchar();
	while(!isdigit(c))c=getchar();
	while(isdigit(c))x=x*10+(c&15),c=getchar();
}
typedef long long ll;
const int N=1e6+7,mod=1e8+7;
int n,m,x,f[N],C;
int qpow(int a,int b=mod-2){
	int s=1;
	for(;b;b>>=1,a=(ll)a*a%mod)if(b&1)s=(ll)s*a%mod;
	return s;
}
int main(){
	read(n),read(m),x=qpow(2,n)-1;
	if(x<0)x+=mod;
	C=(ll)x*(x+mod-1)%mod*qpow(2)%mod;
	for(int i=3;i<=m;i++)
		f[i]=((C-f[i-1]-(ll)f[i-2]*((x-i+2+mod)%mod)%mod)%mod+mod)%mod*qpow(i)%mod,C=(ll)C*(x-i+1+mod)%mod*qpow(i)%mod;
	printf("%d\n",f[m]);
}

TXL

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
根据引用\[1\]和引用\[2\]的描述,题目中的影魔拥有n个灵魂,每个灵魂有一个战斗力ki。对于任意一对灵魂对i,j (i<j),如果不存在ks (i<s<j)大于ki或者kj,则会为影魔提供p1的攻击力。另一种情况是,如果存在一个位置k,满足ki<c<kj或者kj<c<ki,则会为影魔提供p2的攻击力。其他情况下的灵魂对不会为影魔提供攻击力。 根据引用\[3\]的描述,我们可以从左到右进行枚举。对于情况1,当扫到r\[i\]时,更新l\[i\]的贡献。对于情况2.1,当扫到l\[i\]时,更新区间\[i+1,r\[i\]-1\]的贡献。对于情况2.2,当扫到r\[i\]时,更新区间\[l\[i\]+1,i-1\]的贡献。 因此,对于给定的区间\[l,r\],我们可以根据上述方法计算出区间内所有下标二元组i,j (l<=i<j<=r)的贡献之和。 #### 引用[.reference_title] - *1* *3* [P3722 [AH2017/HNOI2017]影魔(树状数组)](https://blog.csdn.net/li_wen_zhuo/article/details/115446022)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [洛谷3722 AH2017/HNOI2017 影魔 线段树 单调栈](https://blog.csdn.net/forever_shi/article/details/119649910)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值