ICode9

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

AT3962 [AGC024E] Sequence Growing Hard 题解

2021-10-14 08:00:07  阅读:212  来源: 互联网

标签:AT3962 Sequence int 题解 305 read sm cdot dp


ATcoder
Luogu

Description.

给定 \(n,k,m\),统计序列序列数 \(\{A_i(i\in[0,n])\}\),使得

  • \(\text{size}(A_i)=i\)
  • \(\forall i\in[1,n],j\in[1,i],A_i(j)\in[1,k]\)
  • \(\forall i\in[1,n]\),\(A_{i-1}\) 是 \(A_i\) 的子序列且字典序小于 \(A_i\)

Solution.

自己的想法:
题意转化成了求不同操作序列数

  1. 每次末尾插入一个任意元素
  2. 每次在第 \(k\) 个位置插入一个权值 \(>a_k\) 的元素

题解:
发现在开头插个 \(0\),上面的 \(1\) 就和 \(2\) 合并了,两者是本质相同的。
然后考虑建树,一个点向它插入位置连边,还要记录一个当前是第几次操作。
本质上来说,一个点当前第几次操作其实在树的形态确定时是基本确定的。
只能在所有儿子改变顺序时才能变。

那我们考虑记录 \(dp_{x,y}\) 表示大小为 \(x\) 根权为 \(y\) 的子树的方案数。
我们枚举第一个子树是什么,然后删掉后 \(dp\) 状态是已知的。
枚举第一个子树的大小和根权,则有

\[\begin{aligned} dp_{x,y}&=\sum_{a=1}^{x-1}\sum_{v=y+1}^k\dbinom{x-2}{a-1}\cdot dp_{x-a,y}\cdot dp_{a,v}\\ &=\sum_{a=1}^{x-1}\dbinom{x-2}{a-1}\cdot dp_{x-a,y}\cdot \sum_{v=y+1}^k dp_{a,v} \end{aligned} \]

Coding.

点击查看代码
//是啊,你就是那只鬼了,所以被你碰到以后,就轮到我变成鬼了{{{
#include<bits/stdc++.h>
using namespace std;typedef long long ll;
template<typename T>inline void read(T &x)
{
	x=0;char c=getchar(),f=0;
	for(;c<48||c>57;c=getchar()) if(!(c^45)) f=1;
	for(;c>=48&&c<=57;c=getchar()) x=(x<<1)+(x<<3)+(c^48);
	f?x=-x:x;
}
template<typename T,typename...L>inline void read(T &x,L&...l) {read(x),read(l...);}//}}}
int n,K,P,dp[305][305],sm[305][305],C[305][305];
inline void upd(int &a,int b) {(a+=b)>=P?a-=P:a;}
int main()
{
	read(n,K,P);for(int i=K;i>=0;i--) sm[1][i]=sm[1][i+1]+(dp[1][i]=1);
	for(int i=0;i<=n;i++) {C[i][0]=1;for(int j=1;j<=i;j++) upd(C[i][j]=C[i-1][j-1],C[i-1][j]);}
	for(int i=2;i<=n+1;i++) for(int j=K;j>=0;upd(sm[i][j]=sm[i][j+1],dp[i][j]),j--)
		for(int l=1;l<i;l++) upd(dp[i][j],1ll*C[i-2][l-1]*dp[i-l][j]%P*sm[l][j+1]%P);
	return printf("%d\n",dp[n+1][0]),0;
}

标签:AT3962,Sequence,int,题解,305,read,sm,cdot,dp
来源: https://www.cnblogs.com/pealfrog/p/15405024.html

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

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

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

ICode9版权所有