Description
m
种工人,
Input
第一行一整数
T
表示用例组数,每组用例输入两整数
Output
对于每组用例,输出完成所有工作至少需要的工人数量
Sample Input
2
4 1
0 0
0 1 1 1 3
1 1 3 3 4
1 0 10 1 5
4 1
0 0
0 1 1 1 3
1 1 3 3 4
1 0 3 1 5
Sample Output
5
9
Solution1
可以看出每种工人之间相互独立,所以可以单独分析,考虑完成所有工作所需第
k
种工人的数量,考虑把工人看作费用,把一项工作所需工人看作流量,那么问题转化为一个费用流,由于工人可以做多项工作,把每个工作拆点
1.源点向i连容量为无穷,费用为 1 的边,体现出每个工人是一个费用
2.
3.源点向
4.
跑一边费用流,最小费用即为所需第
Code1
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
#define maxn 333
#define maxm 66666
#define INF 0x3f3f3f3f
int head[maxn],d[maxn],s,e,no,vis[maxn],pre[maxn];//s为源点,e为汇点
struct point
{
int u,v,flow,next,cost;
point(){};
point(int _u,int _v,int _next,int _flow,int _cost)
{
u=_u,v=_v,next=_next,flow=_flow,cost=_cost;
}
}p[maxm];
void add(int x,int y,int z,int c)//从x到y建一条容量为z,花费为c的边
{
p[no]=point(x,y,head[x],z,c);
head[x]=no++;
p[no]=point(y,x,head[y],0,-c);
head[y]=no++;
}
void init()//初始化
{
memset(head,-1,sizeof(head));
no=0;
}
bool spfa()
{
int i,x,y;
queue<int>q;
memset(d,0x3f,sizeof(d));
memset(vis,false,sizeof(vis));
memset(pre,-1,sizeof(pre));
d[s]=0;
vis[s]=true;
q.push(s);
while(!q.empty())
{
x=q.front();
q.pop();
vis[x]=false;
for(i=head[x];i!=-1;i=p[i].next)
{
if(p[i].flow&&d[y=p[i].v]>d[x]+p[i].cost)
{
d[y]=d[x]+p[i].cost;
pre[y]=i;
if(vis[y])
continue;
vis[y]=true;
q.push(y);
}
}
}
return d[e]!=d[e+1];
}
int mcmf()//最小费用最大流
{
int mincost=0,maxflow=0,minflow,i;
while(spfa())
{
//if(d[e]>=0)break;//可行流
minflow=INF;
for(i=pre[e];i!=-1;i=pre[p[i].u])
minflow=min(minflow,p[i].flow);
for(i=pre[e];i!=-1;i=pre[p[i].u])
{
p[i].flow-=minflow;
p[i^1].flow+=minflow;
}
mincost+=d[e]*minflow;
maxflow+=minflow;
}
return mincost;//最小费用,最大流为maxflow
}
struct node
{
double x,y,b,p;
int num[6];
}a[maxn],o;
double dis(node a,node b)
{
double x=a.x-b.x,y=a.y-b.y;
return sqrt(x*x+y*y);
}
int T,n,m;
int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
scanf("%lf%lf",&o.x,&o.y);
for(int i=1;i<n;i++)
{
scanf("%lf%lf%lf%lf",&a[i].x,&a[i].y,&a[i].b,&a[i].p);
for(int j=1;j<=m;j++)scanf("%d",&a[i].num[j]);
}
s=0,e=2*n-1;
int ans=0;
for(int k=1;k<=m;k++)
{
init();
for(int i=1;i<n;i++)
{
add(s,i,INF,1);
add(i,e,a[i].num[k],0);
add(s,i+n-1,a[i].num[k],0);
for(int j=1;j<n;j++)
if(i!=j&&a[i].b+a[i].p+dis(a[i],a[j])<=a[j].b)
add(i+n-1,j,a[i].num[k],0);
}
ans+=mcmf();
}
printf("%d\n",ans);
}
return 0;
}
Solution2
同样考虑完成所有工作所需第
k
种工人的数量,如果去掉工人可以完成多项工作的条件,那么所需人数应为
1.源点向
i
建容量为
2.
3.
i
向
Code2
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
#define maxn 333
#define maxm 66666
#define INF 0x3f3f3f3f
int head[maxn],cur[maxn],d[maxn],st[maxm],s,e,no;//s为源点,e为汇点,n为点数,no为边数
struct point
{
int u,v,flow,next;
point(){};
point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)//从x到y建容量为z的边
{
p[no]=point(x,y,head[x],z);//前向弧,标号为偶
head[x]=no++;
p[no]=point(y,x,head[y],0);//后向弧,标号为奇
head[y]=no++;
}
void init()//初始化
{
memset(head,-1,sizeof(head));
no=0;
}
bool bfs()
{
int i,x,y;
queue<int>q;
memset(d,-1,sizeof(d));
d[s]=0;
q.push(s);
while(!q.empty())
{
x=q.front();
q.pop();
for(i=head[x];i!=-1;i=p[i].next)
{
if(p[i].flow&& d[y=p[i].v]<0)
{
d[y]=d[x]+1;
if(y==e)
return true;
q.push(y);
}
}
}
return false;
}
int dinic()//最大流
{
int i,loc,top,x=s,nowflow,maxflow=0;
while(bfs())
{
memcpy(cur,head,sizeof(head));
top=0;
while(true)
{
if(x==e)
{
nowflow=INF;
for(i=0;i<top;i++)
{
if(nowflow>p[st[i]].flow)
{
nowflow=p[st[i]].flow;
loc=i;
}
}
for(i=0;i<top;i++)
{
p[st[i]].flow-=nowflow;
p[st[i]^1].flow+=nowflow;
}
maxflow+=nowflow;
top=loc;
x=p[st[top]].u;
}
for(i=cur[x];i!=-1;i=p[i].next)
if(p[i].flow&&d[p[i].v]==d[x]+1)
break;
cur[x]=i;
if(i!=-1)
{
st[top++]=i;
x=p[i].v;
}
else
{
if(!top)
break;
d[x]=-1;
x=p[st[--top]].u;
}
}
}
return maxflow;
}
struct node
{
double x,y,b,p;
int num[6];
}a[maxn],o;
double dis(node a,node b)
{
double x=a.x-b.x,y=a.y-b.y;
return sqrt(x*x+y*y);
}
int T,n,m;
int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
scanf("%lf%lf",&o.x,&o.y);
int ans=0;
for(int i=1;i<n;i++)
{
scanf("%lf%lf%lf%lf",&a[i].x,&a[i].y,&a[i].b,&a[i].p);
for(int j=1;j<=m;j++)scanf("%d",&a[i].num[j]),ans+=a[i].num[j];
}
s=0,e=2*n-1;
for(int k=1;k<=m;k++)
{
init();
for(int i=1;i<n;i++)
{
add(s,i,a[i].num[k]);
add(i+n-1,e,a[i].num[k]);
for(int j=1;j<n;j++)
if(i!=j&&a[i].b+a[i].p+dis(a[i],a[j])<=a[j].b)
add(i,j+n-1,a[i].num[k]);
}
ans-=dinic();
}
printf("%d\n",ans);
}
return 0;
}