평행 사각형 최적화 (HDOJ 3506)

2451 단어 공부 하 다.
사각형 부등식 은 비교적 흔히 볼 수 있 는 동적 계획 을 최적화 하 는 방법 이다.
증명:http://baike.baidu.com/view/1985058.htm?fr=aladdin
이런 문 제 를 해결 하 는 대략적인 절 차 는:
상태 이동 방정식 dp [i] [j] = min {dp [i] [k - 1] + dp [k] [j]} + w [i] [j] (i<=k<=j) 
평행 사각형 으로 최적화 하려 면 w [i] [j], dp [i] [j] 가 사각형 부등식 을 만족 시 키 는 지 증명 해 야 한다.
w [a, c] + w [b, d] < = w [b, c] + w [a, d] (a 볼록 사각형 부등식
또는 시 계 를 쳐 서 w [i] [j + 1] - w [i] i 에 관 한 표현 식 을 관찰 합 니 다. i 체감 에 관 해 서 는 w 가 돌출 사각형 의 부등식 을 만족 시 킵 니 다.
하면, 만약, 만약... w[i'][j]<=w[i][j']  i<=i'<=j<=j' 구간 포함 관계 단조
만약 w 가 사각형 부등식 과 구간 단조 로 운 관 계 를 동시에 만족 시 키 면 dp 도 사각형 부등식 을 만족시킨다.
일반적인 동적 계획 의 복잡 도 는 O (n ^ 3) 이 고 사각형 부등식 프로그램 에서 한 번 뛰 는 i 는 j 만 한 번 뛰 기 때문에 O (n ^ 2) 로 최적화 할 수 있 습 니 다.
상세 증명:
http://baike.baidu.com/view/1985058.htm?fr=aladdin
http://blog.csdn.net/lmyclever/article/details/6677683
전형 적 인 돌 합병 문제:http://blog.csdn.net/acdreamers/article/details/18039073
HDOJ2829
제목 의 대의: 길이 가 n 인 서열 을 정 하고 기껏해야 서열 을 m 단 으로 나 누 며 각 단락 의 서열 은 모두 가중치 가 있 으 며, 가중치 는 서열 내 두 개의 수 와 두 개의 곱 하기 합 이다. m < = n < 1000. 가중치 가 가장 작다.
상태 이동 방정식:
dp[c][i]=min(dp[c][i],dp[c-1][j]+w[j+1][i])
#include 
#include 
#include 
using namespace std;
const int INF=1<<30;
const int MAXN=1000+10;
typedef long long LL;
LL dp[MAXN][MAXN];//dp[c][j]   j    c       
int val[MAXN];
int w[MAXN][MAXN];//w[i][j]  i j      
int s[MAXN][MAXN];//s[c][j]   j     c    
int sum[MAXN];
int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		if(n==0&&m==0)break;
		memset(s,0,sizeof(s));
		memset(w,0,sizeof(w));
		memset(dp,0,sizeof(dp));
		memset(sum,0,sizeof(sum));
		for(int i=1;i<=n;++i)
		{
			scanf("%d",&val[i]);
			sum[i]+=sum[i-1]+val[i];
		}
		for(int i=1;i<=n;++i)
		{
			w[i][i]=0;
			for(int j=i+1;j<=n;++j)
			{
				w[i][j]=w[i][j-1]+val[j]*(sum[j-1]-sum[i-1]);
			}
		}
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=m;++j)
			{
				dp[j][i]=INF;
			}
		}
		for(int i=1;i<=n;++i)
		{
			dp[0][i]=w[1][i];
			s[0][i]=0;
		}
		for(int c=1;c<=m;++c)
		{
			s[c][n+1]=n;//    
			for(int i=n;i>c;--i)
			{
				int tmp=INF,k;
				for(int j=s[c-1][i];j<=s[c][i+1];++j)
				{
					if(dp[c-1][j]+w[j+1][i]

좋은 웹페이지 즐겨찾기