[bzoj4621] Tc605 解题报告
2016-09-23 08:24
176 查看
。。感觉这题跟zj的那题好像。。然后就在想区间dp什么的。。然后突然发现顺着dp就可以了。。
f(i,j,k)表示位置i最终变成了位置j上的数操作了k次的方案数,那么转移就是f(i,j,k)−>f(i+1,j,k)[aj>ai+1]−>f(i+1,j′,k+1)[j′>j][aj′≥max{amin(j′,i+1) max(j′,i+1)}]
但是还有一些蛋疼的问题,就是如果i变成了i,那么它有可能没花费操作,也有可能花费了,所以这个要单独再开一个状态考虑一下。然后关于这个奇怪的状态的转移就需要单独讨论一下。。
代码:
总结:
一定要静下心来慢慢读题慢慢想,一定要把细节都想清楚了再写。
f(i,j,k)表示位置i最终变成了位置j上的数操作了k次的方案数,那么转移就是f(i,j,k)−>f(i+1,j,k)[aj>ai+1]−>f(i+1,j′,k+1)[j′>j][aj′≥max{amin(j′,i+1) max(j′,i+1)}]
但是还有一些蛋疼的问题,就是如果i变成了i,那么它有可能没花费操作,也有可能花费了,所以这个要单独再开一个状态考虑一下。然后关于这个奇怪的状态的转移就需要单独讨论一下。。
代码:
#include<cstdio> #include<cstring> #include<iostream> using namespace std; #include<algorithm> const int Mod=1e9+7; const int N=500+5,K=500+5; int a ,Max ; int f[2][K] ,s[K] ; int main() { freopen("bzoj4621.in","r",stdin); freopen("bzoj4621.out","w",stdout); int n,kmax; scanf("%d%d",&n,&kmax); for(int i=1;i<=n;++i)scanf("%d",a+i); for(int i=1;i<=n;++i) { Max[i][i]=a[i]; for(int j=i+1;j<=n;++j)Max[i][j]=max(Max[i][j-1],a[j]); } f[0][0][0]=1; int now=0,next=1; for(int i=0;i<n;++i) { //printf("-----%d-----\n",i); memset(f[next],0,sizeof(f[next])); for(int k=0;k<=min(i,kmax);++k) { int s=0; for(int j=1;j<=i;++j) { if(a[j]>a[i+1]&&f[now][k][j])(f[next][k][j]+=f[now][k][j])%=Mod; if(a[j]>Max[j+1][i+1])(f[next][k+1][j]+=s)%=Mod; if(f[now][k][j])(s+=f[now][k][j])%=Mod; } if(f[now][k][0])(s+=f[now][k][0])%=Mod; if(a[i]>a[i+1]&&f[now][k][0])(f[next][k+1][i]+=f[now][k][0])%=Mod; (f[next][k][0]+=s)%=Mod; if(f[now][k][i+1]) { (f[next][k][i+1]+=f[now][k][i+1])%=Mod; (s+=f[now][k][i+1])%=Mod; } for(int j=i+2;j<=n;++j) { if(f[now][k][j])(f[next][k][j]+=f[now][k][j])%=Mod; if(a[j]>Max[i+1][j-1])(f[next][k+1][j]+=s)%=Mod; if(f[now][k][j])(s+=f[now][k][j])%=Mod; } } /*for(int k=0;k<=min(i,kmax);++k) for(int j=0;j<=n;++j) if(f[next][k][j]) { printf("f(%d,%d,%d)=%d\n",i+1,k,j,f[next][k][j]); }*/ swap(now,next); } int ans=0; for(int k=0;k<=kmax;++k) for(int j=0;j<=n;++j) if(f[now][k][j]) (ans+=f[now][k][j])%=Mod; printf("%d\n",ans); }
总结:
一定要静下心来慢慢读题慢慢想,一定要把细节都想清楚了再写。
相关文章推荐
- BZOJ 1367 [Baltic2004]sequence 解题报告
- [BZOJ2754] [SCOI2012]喵星球上的点名解题报告|后缀数组
- BZOJ 2351 [BeiJing 2011] 字符串Hash 解题报告
- BZOJ4010: [HNOI2015]菜肴制作 解题报告
- BZOJ3790 神奇项链 解题报告【字符串】【Manacher】【树状数组】【数据结构优化DP】
- BZOJ 4318 期望DP 解题报告
- BZOJ 1379 [Baltic 2001] 解题报告
- BZOJ 4145 [AMPPZ2014] The Prices 解题报告
- BZOJ 4127 Abs 解题报告
- BZOJ 2733 [HNOI 2012] 线段树合并 解题报告
- BZOJ 1093 [ZJOI 2007] Tarjan+DAG拓扑排序DP 解题报告
- [bzoj4314] 倍数?倍数! 解题报告
- [BZOJ1031][JSOI2007]字符加密Cipher 解题报告
- [BZOJ1036][ZJOI2008]树的统计Count 解题报告|树链剖分
- bzoj 1052 覆盖问题 贪心 解题报告
- BZOJ 1011 数学 解题报告
- BZOJ 3996 [TJOI 2015] 线性代数 解题报告
- BZOJ1028[JS0I2007] 麻将 解题报告【模拟+贪心(?)】
- BZOJ 2527 [Poi 2011] 整体二分 解题报告
- BZOJ4444: [Scoi2015]国旗计划 解题报告