Hatena::ブログ(Diary)

tobiasの日記

2013-08-13

2415:Sashimi

問題文
http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=2415











































Monge DP。
週刊spaghetti_sourceさんにわかりやすい説明があったので参考にしながら解いた。
http://topcoder.g.hatena.ne.jp/spaghetti_source/20120915

dp[ 左端 ][ 右端 ] := その部分を全て1cmに切るときの最小コスト
とすると、すでに1cmのときはコスト0だから
dp[ i ][ i ] = 0
区間[ i, j ]を[ i, r ]と[ r + 1, j ]に切るときは
dp[ i ][ j ] = min (dp[ i ][ r ]+dp[ r + 1 ][ j ]+cost[ i ][ j ])
となる。
このままだとO(n^3)で通らないが、cost[ i ][ j ]にはMonge性があってO(n^2)に落ちる。
K[ i ][ j ]=argmin (dp[ i ][ r ]+dp[ r + 1 ][ j ])
とすると
dp[ i ][ j ]={ min ( dp[ i ][ r ]+dp[ r + 1 ][ j ]) | K[ i ][ j - 1] <= r <= K[ i + 1 ][ j ] }

#include<iostream>
#include<vector>
#include<algorithm>
#define INF (1LL<<61)

using namespace std;

typedef unsigned long long ull;

int main(void){
	
	int n;
	cin >> n;
	
	ull cost[4001];
	for(int i=0;i<n;i++){
		cin >> cost[i];
		if(i>0)cost[i]+=cost[i-1];
	}
	
	static ull dp[4001][4001],k[4001][4001];
	fill(dp[0],dp[4001],INF);
	
	for(int i=0;i<n;i++)dp[i][i]=0,k[i][i]=i;
	
	for(int w=1;w<=n;w++){
		for(int i=0,j=i+w;j<n;i++,j++){
			for(int r=k[i][j-1];r<=k[i+1][j];r++){
				ull c=dp[i][r]+dp[r+1][j]+cost[j]-((i>0)?cost[i-1]:0);
				if(dp[i][j]>c)dp[i][j]=c,k[i][j]=r;
			}
		}
	}
	cout << dp[0][n-1] << endl;
	
	return 0;
}

スパム対策のためのダミーです。もし見えても何も入力しないでください
ゲスト


画像認証

トラックバック - http://d.hatena.ne.jp/TobiasGSmollett/20130813/1376371704