bzoj3571: [Hnoi2014]画框
来源:互联网 发布:淘宝网二手app 编辑:程序博客网 时间:2024/06/10 22:39
题目链接
bzoj3571
题目描述
Description
小T准备在家里摆放几幅画,为此他买来了N幅画和N个画框。为了体现他的品味,小T希望能合理地搭配画与画框,使得其显得既不过于平庸也不太违和。对于第 幅画与第 个画框的配对,小T都给出了这个配对的平凡度Aij 与违和度Bij 。整个搭配方案的总体不和谐度为每对画与画框平凡度之和与每对画与画框违和度的乘积。具体来说,设搭配方案中第i幅画与第Pi个画框配对,则总体不和谐度为
小T希望知道通过搭配能得到的最小的总体不和谐度是多少。
Input
输入文件第 行是一个正整数T ,表示数据组数,接下来是T组数据。
对于每组数据,第 行是一个正整数N,表示有N对画和画框。
第2到第N+1行,每行有N个非负整数,第i+1 行第j个数表示Aij 。
第N+2到第2*N+1行,每行有N个非负整数,第i+N+1 行第j个数表示Bij 。
Output
包含T行,每行一个整数,表示最小的总体不和谐度
Sample Input
1
3
4 3 2
2 3 4
3 2 1
2 3 2
2 2 4
1 1 3
Sample Output
30
HINT
第1幅画搭配第3个画框,第2幅画搭配第1个画框,第3 幅画搭配第2个画框,则总体不和谐度为30
N<=70,T<=3,Aij<=200,Bij<=200
题解
我们将每个方案看成平面上的一个点
1.找到平面上y最大和x最大的点(可以用KM算法),将他们连线。
2.找到距离此直线最远的点。与左右的点分别递归下去。怎么找这个点?只要那个点与这两个点构成的三角形面积最大即可。用叉积判断。
注意到
#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<cmath>#include<queue>#include<vector>using namespace std;#define N 80typedef pair<int,int>P;P l,r;int tx[N],ty[N],a[N][N],b[N][N],d[N][N],match[N],n,fix,T;bool vx[N],vy[N];bool find(int x){ vx[x]=true; for(int y=1;y<=n;y++) if(!vy[y]){ int tmp=tx[x]+ty[y]-d[x][y]; if(!tmp) { vy[y]=true; if(!match[y]||find(match[y])){ match[y]=x; return true; } } else fix=min(fix,tmp); } return false;}P km(){ memset(tx,192,sizeof(int)*(n+1)); memset(ty,0,sizeof(int)*(n+1)); memset(match,0,sizeof(int)*(n+1)); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) tx[i]=max(tx[i],d[i][j]); for(int x=1;x<=n;x++){ while(true){ memset(vy,false,sizeof(bool)*(n+1)); memset(vx,false,sizeof(bool)*(n+1)); fix=1000000007; if(find(x)) break; for(int i=1;i<=n;i++){ if(vx[i]) tx[i]-=fix; if(vy[i]) ty[i]+=fix; } } } P ans(0,0); for(int i=1;i<=n;i++) ans.first+=a[match[i]][i],ans.second+=b[match[i]][i]; return ans;}int solve(P l,P r){ for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) d[i][j]=a[i][j]*(r.second-l.second)+b[i][j]*(l.first-r.first); P mid=km(); if(mid==l||mid==r) return min(l.first*l.second,r.first*r.second); return min(solve(l,mid),solve(mid,r));}int main(){ scanf("%d",&T); while(T--){ scanf("%d",&n); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&a[i][j]); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&b[i][j]); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) d[i][j]=-a[i][j]; l=km(); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) d[i][j]=-b[i][j]; r=km(); printf("%d\n",solve(l,r)); } return 0;}
- bzoj3571: [Hnoi2014]画框
- BZOJ3571: [Hnoi2014]画框
- BZOJ3571[HNOI2014]画框
- [BZOJ3571][HNOI2014]画框
- bzoj3571: [Hnoi2014]画框
- bzoj3571: [Hnoi2014]画框
- 【BZOJ3571】[Hnoi2014]画框【最小乘积匹配】
- 【KM算法】【最大乘积生成树】[HNOI2014] bzoj3571 画框
- 【BZOJ 3571】 [Hnoi2014]画框
- BZOJ 3571: [Hnoi2014]画框
- 3571: [Hnoi2014]画框
- bzoj 3571: [Hnoi2014]画框 最优乘积匹配
- bzoj3571
- [最小乘积匹配 分治 KM] BZOJ 3571 [Hnoi2014]画框
- 画框
- 画框钉角机
- 【HNOI2014】Jabberwocky
- BZOJ3571【最小乘积最大匹配】
- SQLite清空表并将自增列归零
- QT html解析
- java中初始化与构造器
- java笔记-7
- 实现用画笔画出一个带有渐变效果的线条
- bzoj3571: [Hnoi2014]画框
- MySQL安装版本
- Android之–Activity通过Intent传递Object
- java笔记-8
- Jira平台数据自动化获取(一)--Jira接口请求
- TCP协议中的SO_LINGER选项
- 从零开始学习EasyDarwin(RTSP之TaskThread)
- 将带图片的Java工程用eclipse打包成jar,然后用exe4j将jar打包成exe,并附带jre库
- java笔记-9