去年时还是渣渣中的渣渣。。只记得比赛的时候过这个题的不过,今天拿来做一做发现构图还是挺简单的。
首先要突破的就是5种工人是相互独立的,因为每一个任务的起始和终止时间是固定的,我们用费用流的话又可以保证每个任务都完成(满足最大流),那么每种工人怎么跑是互不影响的,所以做m次费用流,结果为之和。
我的构图有2*n个点,起点为0,终点为2*n-1
流量表示工人数量。
每个维修点拆成两个点,但是和大多数网络流不同,左边的点是用来流那些已经执行完该点任务的工人。
首先0点与左边的n-1个点连流量为各点工人数,费用为0的边。
然后右边的n-1个点与终点连流量为各点工人数,费用为0的边,来确保最大流。
如果0点到某个维修点可达(符合时间要求),那么从0到右边的对应点连流量为无穷,费用为1的边。
如果i点到j点可达(符合时间要求),那么从左边的i点 到 右边的 j点连流量为无穷,费用为0的边。
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<stdlib.h>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define ll __int64
int cal[6][160];
int x[160],y[160],t1[160],t2[160];
double dist(double x1,double y1,double x2,double y2)
{
return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
const int N = 500;//点
const int M = 2 * 30010;//边
const int inf = 1000000000;
struct Node{//边,点f到点t,流量为c,费用为w
int f, t, c, w;
}e[M];
int next1[M], point[N], dis[N], q[N], pre[N], ne;//ne为已添加的边数,next,point为邻接表,dis为花费,pre为父亲节点
bool u[N];
void init(){
memset(point, -1, sizeof(point));
ne = 0;
}
void add_edge(int f, int t, int d1, int d2, int w){//f到t的一条边,流量为d1,反向流量d2,花费w,反向边花费-w(可以反悔)
e[ne].f = f, e[ne].t = t, e[ne].c = d1, e[ne].w = w;
next1[ne] = point[f], point[f] = ne++;
e[ne].f = t, e[ne].t = f, e[ne].c = d2, e[ne].w = -w;
next1[ne] = point[t], point[t] = ne++;
}
bool spfa(int s, int t, int n){
int i, tmp, l, r;
memset(pre, -1, sizeof(pre));
for(i = 0; i < n; ++i)
dis[i] = inf;
dis[s] = 0;
q[0] = s;
l = 0, r = 1;
u[s] = true;
while(l != r) {
tmp = q[l];
l = (l + 1) % (n + 1);
u[tmp] = false;
for(i = point[tmp]; i != -1; i = next1[i]) {
if(e[i].c && dis[e[i].t] > dis[tmp] + e[i].w) {
dis[e[i].t] = dis[tmp] + e[i].w;
pre[e[i].t] = i;
if(!u[e[i].t]) {
u[e[i].t] = true;
q[r] = e[i].t;
r = (r + 1) % (n + 1);
}
}
}
}
if(pre[t] == -1)
return false;
return true;
}
void MCMF(int s, int t, int n, int &flow, int &cost){//起点s,终点t,点数n,最大流flow,最小花费cost
int tmp, arg;
flow = cost = 0;
while(spfa(s, t, n)) {
arg = inf, tmp = t;
while(tmp != s) {
arg = min(arg, e[pre[tmp]].c);
tmp = e[pre[tmp]].f;
}
tmp = t;
while(tmp != s) {
e[pre[tmp]].c -= arg;
e[pre[tmp] ^ 1].c += arg;
tmp = e[pre[tmp]].f;
}
flow += arg;
cost += arg * dis[t];
}
}
//建图前运行init()
//节点下标从0开始
//加边时运行add_edge(a,b,c,0,d)表示加一条a到b的流量为c花费为d的边(注意花费为单位流量花费)
//特别注意双向边,运行add_edge(a,b,c,0,d),add_edge(b,a,c,0,d)较好,不要只运行一次add_edge(a,b,c,c,d),费用会不对。
//求解时代入MCMF(s,t,n,v1,v2),表示起点为s,终点为t,点数为n的图中,最大流为v1,最大花费为v2
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
int n,m,i,j,k,ans=0,ta,tb;
scanf("%d%d",&n,&m);
scanf("%d%d",&x[0],&y[0]);
for(i=1;i<=n-1;i++)
{
scanf("%d%d%d%d",&x[i],&y[i],&t1[i],&t2[i]);
for(j=1;j<=m;j++)
scanf("%d",&cal[j][i]);
}
for(k=1;k<=m;k++)
{
init();
for(i=1;i<=n-1;i++)
{
if(dist((double)x[0],(double)y[0],(double)x[i],(double)y[i]) <=t1[i] )
add_edge(0,i+n-1,inf,0,1);
add_edge(0,i,cal[k][i],0,0);
add_edge(n-1+i,n+n-1,cal[k][i],0,0);
}
for(i=1;i<=n-1;i++)
for(j=1;j<=n-1;j++)
{
if(i!=j)
{
if(t1[i]+t2[i]+dist((double)x[j],(double)y[j],(double)x[i],(double)y[i]) <=t1[j])
add_edge(i,n-1+j,inf,0,0);
}
}
MCMF(0,n+n-1,n+n,ta,tb);
ans+=tb;
}
printf("%d\n",ans);
}
return 0;
}