[bzoj4621] Tc605 解题报告

。。感觉这题跟zj的那题好像。。然后就在想区间dp什么的。。然后突然发现顺着dp就可以了。。
f(i,j,k)表示位置i最终变成了位置j上的数操作了k次的方案数，那么转移就是

$$f(i,j,k)\\->f(i+1,j,k)[a_j>a_{i+1}]\\->f(i+1,j',k+1)[j'>j][a_{j'}\ge max\{ a_{min(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[N],Max[N][N];
int f[2][K][N],s[K][N];
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);
}

总结：
一定要静下心来慢慢读题慢慢想，一定要把细节都想清楚了再写。