E. Phoenix and Computers（组合+计数dp)

https://codeforces.com/problemset/problem/1515/E

---

思路：

首先会发现结果的序列方案数会是一段连续的，一个空的，一段连续的，那我们要先知道一段连续的方案数是多少。

枚举从哪一台电脑开始打开

如果从 1 开始，那么剩下的必须按照 2,3,…,n的顺序打开。（如果不按照这个顺序，就会有电脑自动开启） 方案数为1

如果从 2开始，那么 2 的右边必须按从左往右的顺序开启，即 3,4,…,n而 1 可以在这其中任意时刻开启。   后面那么多数的位置插一个空 C(k-1,1)

如果从3开始，那么2的右边必须按左到右的顺序开启，即 4 5 6,,,n。但是1 2可以随意插进去，但必须保证1 和2 的相对顺序不变。

那么此时的方案数是多少？

举个例子1 2 3 4 5，现在枚举3开始， 右边还有4 5，左边要1 2。方案数一共有4个位置待选，我们任意选择两个位置，C(4,2)，代表给1 2定两个位置。此时发现，4 5只有一种顺序，所以剩下的坑的位置是已经固定了。C的定义，位置和顺序无关，也就是4 5或者5 4在这两个位置的顺序是只算一遍的，和题意符合，因为4 5满足的答案只有一种顺序。

那么对应这个方案数就是C(K-1,2)

总和起来就是：C(n−1,0)​+C(n−1,1)​+C(n−1,2​)+…+C(n−1,n−1)​=2^n−1  就是高中的二项式展开。

 

那么这时候再来考虑一段一段的方案数怎么统计。

dp[len][cnt]：表示当前开了len个灯，其中第len个是自动开的，有 cnt个是手动开的。

那么顺推的转移就是 dp[len+k+1] [cnt+k] = Σdp[len][k]*2^(k-1)*C(cnt+k,k)

再枚举开k盏灯，此时len+k+1按照定义就是自动开的，对于开出的连续k盏，刚才的结论就是2^(k-1)。然后总共开了cnt+k个，里面k个新开的随意放哪个顺序最后答案都是符合的。

由于n^3logn太大，组合数部分预处理。

O(n^3)

不过我还看到了O(n^2)的插头dp.....（太菜了不会阿

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=4e2+100;
typedef long long LL;
inline LL read(){LL x=0,f=1;char ch=getchar();	while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL fac[maxn];
LL dp[maxn][maxn];
LL c[maxn][maxn];
LL n,mod;
void pre(){
     c[0][0]=1;
     for(LL i=1;i<maxn;i++){
         for(LL j=0;j<=i;j++){
             c[i][j]=(c[i-1][j-1]%mod+c[i-1][j]%mod)%mod;
         }
     }
}
int main(void){
   cin.tie(0);std::ios::sync_with_stdio(false);
   cin>>n>>mod;
   fac[0]=1;
   for(LL i=1;i<maxn;i++) fac[i]=fac[i-1]*2%mod;
   pre();
   dp[0][0]=1;
   for(LL len=0;len<=n;len++){
       for(LL cnt=0;cnt<=len;cnt++){
           for(LL k=1;len+k<=n;k++){
               dp[len+k+1][cnt+k]+=( ( (dp[len][cnt]%mod*fac[k-1]%mod)%mod*c[cnt+k][k])%mod)%mod;
           }
       }
   }
   LL ans=0;
   for(LL cnt=0;cnt<=n;cnt++){
       ans=(ans%mod+dp[n+1][cnt]%mod)%mod;
   }
   cout<<ans<<"\n";
   return 0;
}