【お昼寝G】リングdpテンプレート質問

[USACO05JAN]お昼寝G

質問の説明

ベッシーはとても退廃的な牛です。彼女の 1 日は N N に均等に分割されます。N セクション( 3 ≤ N ≤ 3830 3 \leq N \leq 3830 3N3830),但是她要用其中的 B B B フェーズ時間( 2 ≤ B <; N 2 \leq B \lt N 2B<;N) は台無しになります。各期間には利用価値があります U i U_i i 0 ≤ U i ≤ 2 × 1 0 5 0 \leq U_i \leq 2 \times 10^5 0i 2×105)、壊れた場合にのみ効果があります。

目覚まし時計の助けを借りて、ベッシーはいつでも振り子を開く時間を選択できますが、もちろん、彼女は時間の境界でのみ振り子と本を開くことができます。

ベッシーは、すべての廃棄物の有用性の合計を最大化したいと考えています。残念ながら、各スイング期間の最初の期間は「オープンスイング」段階であり、利用価値は記録されません。

タイム ステージは、Bessie が間に合えば継続的に循環する円です (毎日がサイクル) N N N 和时间 1 1 1 失敗した場合は時間を与えられます 1 1 1 の利用価値。

入力フォーマット

最初の行の 2 つの整数 N , B N,B NB

接下来 N N N 行。各行には i i i 期間の効用価値。

出力フォーマット

最大の利用価値を出力します。

例1

サンプル入力 #1

5 3
2
0
3
1
4

サンプル出力 #1

6

ヒント

从第 4 4 4 期間が開始し、 1 1 この本は 1 期間後に開かれます。

アイデア

被害が 2 日間にわたる状況、つまり、
ステータスf[i,j,k]: 前< を考慮せずに、単純化した問題を考えてみましょう。 a i =3> i i i 个小时中一共睡了 j j j 个小时,且第 i i i 時間の廃墟状態は k k k (1 はスリープし、0 はスリープしません)。
状態遷移:

F [ i , j , 0 ] = MAX ( F [ i − 1 , j , 0 ] , F [ i − 1 , j , 1 ] ) F[i,j,0]=MAX(F[i-1, j,0],F[i-1,j,1])F[i,j0]=MAX(F [i1j0],F[i1j1])
F [ i , j , 1 ] = M A X ( F [ i − 1 , j − 1 , 0 ] , F [ i − 1 , j − 1 , 1 ] + U i ) F[i,j,1]=MAX(F[i-1,j-1,0],F [i-1,j-1,1]+U_i) F[i,j1]=MAX(F [i1j10],F[i1j11]+i )

わかりやすいのですが、この場合問題があります。スイングが始まるのは毎日早くても1時で、このときの利用価値は集計に含めることができません。

次に、別の状況を考えて、強制的に No にします。 N N N 个小时和当天第 1 1 1个小时摆烂,那么第 1 1 1 時間の利用価値を計算できます。

分割した 2 つの部分を dp に分割し、最後に M A X MAX を取ります。MAX 即为答案。

コード:

#include<iostream>
#include<cstring>
using namespace std;
int n,b;
int val[3850];
int dp[3850][3850][2],ans;

int main(){
    
    
	cin>>n>>b;
	for(int i=1;i<=n;i++)cin>>val[i];
	
	memset(dp,-0x3f,sizeof dp);
	dp[1][1][1]=dp[1][0][0]=0;
	
	for(int i=2;i<=n;i++){
    
    
		dp[i][0][0]=dp[i-1][0][0];
		for(int j=1;j<=b;j++){
    
    
			dp[i][j][0]=max(dp[i-1][j][0],dp[i-1][j][1]);
			dp[i][j][1]=max(dp[i-1][j-1][0],dp[i-1][j-1][1]+val[i]);
		}
	}
	ans=max(dp[n][b][0],dp[n][b][1]);
	
	memset(dp,-0x3f,sizeof dp);
	dp[1][1][1]=val[1];//第1个小时也记上了!
	dp[1][0][0]=0;
	
	for(int i=2;i<=n;i++){
    
    
		dp[i][0][0]=dp[i-1][0][0];
		for(int j=1;j<=b;j++){
    
    
			dp[i][j][0]=max(dp[i-1][j][0],dp[i-1][j][1]);
			dp[i][j][1]=max(dp[i-1][j-1][0],dp[i-1][j-1][1]+val[i]);
		}
	}
	ans=max(ans,dp[n][b][1]);//第N个小时必须得摆!
	
	cout<<ans;
	
	return 0;
}

上記に基づいて、ローリング配列を最適化して空間を複素数にします n 2 → n n^2→n n2n

#include<iostream>
#include<cstring>
using namespace std;
int n, b;
int val[3850];
int dp[2][3850][2], ans;

int main() {
    
    
	cin >> n >> b;
	for (int i = 1; i <= n; i++)cin >> val[i];

	memset(dp, -0x3f, sizeof dp);
	dp[1][1][1] = dp[1][0][0] = 0;
	
	for (int i = 2; i <= n; i++) {
    
    
		for (int j = 0; j <= b; j++) {
    
    
			dp[i & 1][j][0] = max(dp[(i - 1) & 1][j][0], dp[(i - 1) & 1][j][1]);
			if(j >= 1)dp[i & 1][j][1] = max(dp[(i - 1) & 1][j - 1][0], dp[(i - 1) & 1][j - 1][1] + val[i]);
		}
	}
	ans = dp[n & 1][b][0];

	memset(dp, -0x3f, sizeof dp);
	dp[1][1][1] = val[1]; //第1个小时也记上了!
	dp[1][0][0] = 0;
	
	for (int i = 2; i <= n; i++) {
    
    
		for (int j = 0; j <= b; j++) {
    
    
			dp[i & 1][j][0] = max(dp[(i - 1) & 1][j][0], dp[(i - 1) & 1][j][1]);
			if(j >= 1)dp[i & 1][j][1] = max(dp[(i - 1) & 1][j - 1][0], dp[(i - 1) & 1][j - 1][1] + val[i]);
		}
	}
	ans = max(ans, dp[n & 1][b][1]); //第N个小时必须得摆!

	cout << ans;

	return 0;
}

おすすめ

転載: blog.csdn.net/skyflying266/article/details/125343020