POJ-3216Repairing Company(最小路径覆盖)
题目意思是说有个公司负责地区的维修任务,该公司负责Q个地区。
现在有M个维修任务需要处理,需要派师傅去处理。
每个维修任务由3个特性,P维修任务所在地区,T维修任务开始时间(只能刚刚好这个时间,不可提前),D该维修任务持续时间。
现在问要完成所有维修任务最少要派多少个师傅。
输入数据先给出公司负责的地区数目,然后给出一个Q*Q的矩阵,代表各个城市地区之间的可达情况。如果是-1代表不可达。
然后给出M个任务信息。
如果师傅完成了任务一还有时间完成任务二,那么完成这两个任务只需要派出一个师傅就可以解决。
现在我们构造一个有向无环图(DAG)
如果 一个任务的开始时间+持续时间+到达下一个任务所在地区的时间<=下一个任务开始时间
这样说明这两个任务可以由同一个人完成。
那么在这一个DAG图中,如果选择一条边,那么就会两个顶点被覆盖,完成了两个任务。
那么题目就转变成最小路径覆盖(选择边,把所有顶点覆盖)
因为 最小路径覆盖+最大匹配=顶点个数
那么我们就需要构造一个二分图,求得最大匹配。
现在我们把每个顶点拆分为两分,一个放在集合X中,一个放集合Y中。两个集合如果有前驱和后继的关系就连一条边,也就是代表这两个任务可由同一个人完成。
现在有M个维修任务需要处理,需要派师傅去处理。
每个维修任务由3个特性,P维修任务所在地区,T维修任务开始时间(只能刚刚好这个时间,不可提前),D该维修任务持续时间。
现在问要完成所有维修任务最少要派多少个师傅。
输入数据先给出公司负责的地区数目,然后给出一个Q*Q的矩阵,代表各个城市地区之间的可达情况。如果是-1代表不可达。
然后给出M个任务信息。
如果师傅完成了任务一还有时间完成任务二,那么完成这两个任务只需要派出一个师傅就可以解决。
现在我们构造一个有向无环图(DAG)
如果 一个任务的开始时间+持续时间+到达下一个任务所在地区的时间<=下一个任务开始时间
这样说明这两个任务可以由同一个人完成。
那么在这一个DAG图中,如果选择一条边,那么就会两个顶点被覆盖,完成了两个任务。
那么题目就转变成最小路径覆盖(选择边,把所有顶点覆盖)
因为 最小路径覆盖+最大匹配=顶点个数
那么我们就需要构造一个二分图,求得最大匹配。
现在我们把每个顶点拆分为两分,一个放在集合X中,一个放集合Y中。两个集合如果有前驱和后继的关系就连一条边,也就是代表这两个任务可由同一个人完成。
继而求出最大匹配以后,用V-最大匹配就等于最小路径覆盖
#include<iostream>
#include<cstdio>
#include<cstring>
#include<functional>
#include<algorithm>
#define inf 0x3f3f3f3f
#define maxn 1005
using namespace std;
int maps[maxn][maxn];
int e[maxn][maxn];
int cy[maxn],vis[maxn];
int q,m,ans;
struct node
{
int num,st,ct;
}task[maxn];
void floyd()//求任意两地区之间最短距离
{
int i,j,k;
for(k=1;k<=q;k++)
for(i=1;i<=q;i++)
for(j=1;j<=q;j++)
{
if(maps[i][j]>maps[i][k]+maps[k][j])
maps[i][j]=maps[i][k]+maps[k][j];
}
}
int path(int u)
{
for(int v=1;v<=m;v++)
{
if(e[u][v]&&!vis[v])
{
vis[v]=1;
if(cy[v]==0||path(cy[v]))
{
cy[v]=u;
return 1;
}
}
}
return 0;
}
int MaxMatch()
{
memset(cy,0,sizeof(cy));
int sum=0;
for(int i=1;i<=m;i++)
{
memset(vis,0,sizeof(vis));
sum+=path(i);
}
return sum;
}
int main()
{
int i,j;
int a,b,c;
while(scanf("%d%d",&q,&m),q&&m)
{
memset(e,0,sizeof(e));
for(i=1;i<=q;i++)//读入地区可达信息
for(j=1;j<=q;j++)
{
scanf("%d",&maps[i][j]);
if(maps[i][j]==-1)
maps[i][j]=inf;
}
floyd();
for(i=1;i<=m;i++)//读入任务信息
{
scanf("%d%d%d",&a,&b,&c);
task[i].num=a;
task[i].st=b;
task[i].ct=c;
}
for(i=1;i<=m;i++)//构造二分图
for(j=1;j<=m;j++)
{
if(i==j)
continue;
if(task[i].st+task[i].ct+maps[task[i].num][task[j].num]<=task[j].st)
e[i][j]=1;
}
ans=MaxMatch();//最大匹配
printf("%d\n",m-ans);
}
return 0;
}