四边形优化

匆匆忙忙搞了一下四边形优化，也就是做了几道入门题而已

四边形不等式详解：

http://www.cnblogs.com/zxndgv/archive/2011/08/02/2125242.html

反正我就记住这句话： 判断w是否为凸即判断 w[i,j+1]-w[i,j]的值随着i的增加是否递减 

hdu 2829

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
#define maxn 1111

int dp[maxn][maxn];
int w[maxn][maxn];
int A[maxn], s[maxn];
int t[maxn][maxn];

int main()
{
	int n, m;
	while(scanf("%d %d",&n,&m)!=EOF && (n|| m))
	{
		memset(s, 0, sizeof s);
		memset(w, 0, sizeof w);
		for(int i= 1; i<= n; i++)
		{
			scanf("%d",&A[i]);	
			s[i]= s[i-1]+ A[i];
		}
		for(int i= 1; i<= n; i++)
			for(int j= i+1; j<= n; j++)	
				w[i][j]= w[i][j-1]+ A[j]*(s[j-1]- s[i-1]);
		memset(dp, -1, sizeof dp);		
		for(int i= 1; i<= n; i++)
		{
			t[i][0]= 0;
			t[n+1][i]= n; //这两个是t的边界的处理 
			dp[i][0]= w[1][i];
		}
		for(int j= 1; j<= m; j++)
			for(int i= n; i>= j; i--)
			{
				int a= t[i][j-1];
				int b= t[i+1][j];
				for(int k= a; k<= b; k++)
				if(k> j-1 && k+1<= i) //这里可不要 
				{
					int temp= dp[k][j-1]+ w[k+1][i];
					if(dp[i][j]== -1 || temp< dp[i][j])
					{
						dp[i][j]= temp;
						t[i][j]= k;
					}
				}
			}
		printf("%d\n",dp[n][m]);	
	}
	return 0;
}

hdu 3480 代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 10011
#define maxm 5011

int A[maxn];
int t[maxn][maxm], dp[maxn][maxm];  

int main()
{
//	freopen("in.txt","r",stdin);
	int T;
	scanf("%d",&T);
	for(int C= 1; C<= T; C++)
	{
		int n, m;
		scanf("%d %d",&n,&m);
		for(int i= 1; i<= n; i++)
			scanf("%d",&A[i]);
		sort(A+1, A+n+1);
		memset(dp, -1, sizeof dp);
	//	memset(t, 0, sizeof t); //t数组初始化会TLE，不知道为什么不用初始化 
		for(int i= 1; i<= n; i++)
		{
			t[n+1][i]= n;
			t[i][1]= 1;
			dp[i][1]= (A[i]- A[1])*(A[i]- A[1]);
		}
		for(int j= 2; j<= m; j++)
			for(int i= n; i>= j; i--)
			{
				int a= t[i][j-1];
				int b= t[i+1][j];
				for(int k= a; k<= b; k++)
				{
					int temp= dp[k][j-1]+ (A[i]- A[k+1])*(A[i]- A[k+1]);
					if(dp[i][j]== -1 || temp< dp[i][j])
					{
						dp[i][j]= temp;
						t[i][j]= k;
					}
				}
			}
		int ans;	
		if(m== 0)
			ans= dp[n][1]; //m==0的时候选出0个子集ans不应该为0，我以为结果是0 
		else
			ans= dp[n][m];	
		printf("Case %d: %d\n",C,ans);				
	}
	return 0;
}