bzoj3042 Acting Cute （DP）

bzoj3042 Acting Cute

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=3042

题意：
(权限题贴题面…)
正在rainbow的城堡游玩的freda恰好看见了在地毯上跳舞卖萌的水叮当……于是……
freda：“呜咕>_< 我也要卖萌T_T！”
rainbow给了freda N秒的自由活动时间，不过由于刚刚游览城堡有些累了，freda只想花B秒的时间来卖萌，剩下的时间她要在rainbow的城堡里睡个好觉好好休息一下。
rainbow给这N秒每秒定义了一个值Ui，如果第i秒钟freda在卖萌，那么她可以获得Ui点卖萌指数lala~
freda开始卖萌后可以随时停止，休息一会儿之后再开始。不过每次freda开始卖萌时，都需要1秒来准备= =，这一秒是不能获得卖萌指数的。当然，freda卖萌和准备的总时间不能超过B。
更特殊的是，这N秒钟时间是环形的。也就是freda可以从任意时间开始她的自由活动并持续N秒。
为了使自己表现得比水叮当更萌，现在freda想知道，她最多能获得多少卖萌指数呢？

数据范围
0<=B<=N<=3600，0<=Ui<=200000。

题解：

0<=B<=N<=3600，
对于一般的给出的N的区间选j个点来取得分值这样的问题，容易想到：
dp[i][j]前i秒卖萌j秒且第i秒卖萌的最大值
N*N*B->TLE B*N -> MLE

容易想到的处理：
dp[i][j][0/1]前i秒卖萌j秒且第i秒不卖萌/卖萌的最大值
dp[i][j][1]=max(dp[i-1][j-1][1]+w[i],dp[i-1][j-1][0])
dp[i][j][0]=max(dp[i-1][j][0],dp[i-1][j][1])
发现只由i-1转移过来，可滚动

但是这是一个环？在后面复制不好处理。

实际上这个处理环的方法有些不同，
由题意知，如果第一秒卖了萌且不是准备，那么n秒必须卖萌或准备，如果第一秒没卖萌或第一秒准备，则无所谓。
最后只需要只赋相应的初值，分别跑一次即可。

（提供了一种思路，就是对于环，不一定要像一般的复制一个接在后面，这样由于头的情况不止一个（卖萌/准备），反而把问题复杂化。因为连接的情况就是头尾一致/不一致，可以分别处理）

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
int inf;
const int N=3605;
int n,m,dp[2][N][2],w[N];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    scanf("%d",&w[i]);
    int now=1; int pre;
    memset(dp,-63,sizeof(dp));
    dp[now][1][1]=w[1];//第1秒卖萌
    now=now^1;
    int ans=0;
    for(int i=2;i<=n;i++)
    {
        pre=now^1;
        for(int j=0;j<=min(i,m);j++)
        {
            if(j) dp[now][j][1]=max(dp[pre][j-1][1]+w[i],dp[pre][j-1][0]);
            dp[now][j][0]=max(dp[pre][j][0],dp[pre][j][1]);
        }   
        now=now^1;

    }
    ans=max(ans,dp[now^1][m][1]);
    memset(dp,-63,sizeof(dp));
    now=1; dp[now][0][0]=dp[now][1][1]=0;//第1秒不卖萌或准备 
    now=now^1;
    for(int i=2;i<=n;i++)
    {
        pre=now^1;
        for(int j=0;j<=min(i,m);j++)
        {
            if(j) dp[now][j][1]=max(dp[pre][j-1][1]+w[i],dp[pre][j-1][0]);
            dp[now][j][0]=max(dp[pre][j][0],dp[pre][j][1]);
        }   
        now=now^1;
    }
    ans=max(ans,max(dp[now^1][m][1],dp[now^1][m][0]));
    printf("%d\n",ans); 
    return 0;
}