BZOJ 2711: [Violet 2]After 17

Description

Input

Output

Sample Input

4
4 5
1 2
3 3
4 1

Sample Output

-38.00

HINT

题解

一道很巧妙的dp。

要我们最小化sigma(xi*xj+yi*yj)。原式可化为：[(x1+x2+...+xn)^2-(x1^2+x2^2+...xn^2)+(y1+y2+...+yn)^2-(y1^2+y2^2+...+yn^2)]/2应为要最小化，x1^2+x2^2+...xn^2一定是取题目所给的那个“边界值”时最大，y1^2+y2^2+...+yn^2同理。此时就要求我们要得到绝对值最小的x1+x2+...+xn与y1+y2+...+yn。一个简单的dp就可以了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
#define ll long long
#define MAXN 40000
using namespace std;
int n;
int x[202],y[202],a[202],f[202][80002];
ll ans;
void init()
{
	scanf("%d",&n);
	int i;
	for(i=1;i<=n;i++)
	   {scanf("%d%d",&x[i],&y[i]);
	    ans-=x[i]*x[i]+y[i]*y[i];
	   }
}
void dp()
{
	int i,j;
	for(i=1;i<=n;i++) a[i]=x[i];
	f[0][MAXN]=1;
	for(i=1;i<=n;i++)
	for(j=-MAXN;j<=MAXN;j++)
	   {if(f[i-1][j+MAXN])
	       {f[i][j+MAXN+a[i]]=f[i][j+MAXN-a[i]]=1;}
	   }
	for(i=0;i<=MAXN;i++)
	   {if(f[n][MAXN+i]||f[n][MAXN-i]) {ans+=i*i; break;}
	   }
	memset(f,0,sizeof(f));
	for(i=1;i<=n;i++) a[i]=y[i];
	f[0][MAXN]=1;
	for(i=1;i<=n;i++)
	for(j=-MAXN;j<=MAXN;j++)
	   {if(f[i-1][j+MAXN])
	       {f[i][j+MAXN+a[i]]=f[i][j+MAXN-a[i]]=1;}
	   }
	 for(i=0;i<=MAXN;i++)
	   {if(f[n][MAXN+i]||f[n][MAXN-i]) {ans+=i*i; break;}
	   }  
	printf("%.2lf\n",double(ans/2));
}
int main()
{
	init(); dp();
	return 0;
}