POJ 3216 Repairing Company(FLOYD+DAG最小路径覆盖)

POJ 3216 Repairing Company(FLOYD+DAG最小路径覆盖)

http://poj.org/problem?id=3216

题意:

       给出Q的街道和M个任务

       然后给出i*j的矩阵..表示第i个地点到第j个地点的距离 其中-1表示不可到达.

       然后接下来M行有 p t d 表示 任务在p地点, 开始时间是t, 完成工作花费时间是d.

       问最少派出多少人可以完成M个任务

分析:

       本题与刘汝佳<<训练指南>>P356例题29基本一样.

       其实如果一个工人做完第一个任务时,还有任务没人做的话且他能及时赶到这个任务的话,那么他可以立马跑到这个新任务处去做这个新任务. 那么这样就可以节省人力.

       把每个任务看成一个节点,如果一个人做i任务结束后还能赶到j任务处去做任务,那么就连一条从i到j的有向边. 最终我们得到了一个DAG图(时间是天然的序,所以该图不会存在环). 注意判断任务i是否存在到j的有向边, 需要判断是否 i的结束时间+i到j的最短路径时间<= j的开始时间. 而i到j的最短路径时间需要用Floyd算法求出.

       对于该DAG,我们要找出其最小路径覆盖(因为最小路径覆盖就是我们需要派遣的最少工人数目).

       DAG的最小路径覆盖=n-二分图的最大匹配数. 即本题转化为求二分图的最大匹配问题了.

AC代码:

#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=306+5;
#define inf 0x3f3f3f3f
struct Max_Match
{
    int n,m;
    vector<int>g[maxn];
    bool vis[maxn];
    int left[maxn];

    void init(int n,int m)
    {
        this->n=n;
        this->m=m;
        for(int i=1;i<=n;i++)
			g[i].clear();
        memset(left,-1,sizeof(left));
    }

    bool match(int u)
    {
        for(int i=0;i<g[u].size();i++)
		{
			int v=g[u][i];
			if(!vis[v])
			{
				vis[v]=true;
				if(left[v]==-1 || match(left[v]))
				{
					left[v]=u;
					return true;
				}
			}
		}

        return false;
    }

    int solve()
    {
        int ans=0;
        for(int i=1;i<=n;i++)
        {
            memset(vis,0,sizeof(vis));
            if(match(i)) ans++;
        }
        return ans;
    }
}MM;

int maps[maxn][maxn];
struct node
{
	int place,start,during;
}s[maxn];
void floyd(int m)
{
	for(int k=1;k<=m;k++)
	for(int i=1;i<=m;i++)
			for(int j=1;j<=m;j++)
	{
		if(maps[i][k]!=inf&&maps[k][j]!=inf)
		{
			maps[i][j]=min(maps[i][j],maps[i][k]+maps[k][j]);
		}
	}
}
int main()
{
	int n,m;
    while(~scanf("%d%d",&n,&m))
	{
		if(n==0&&m==0)break;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
			{
				scanf("%d",&maps[i][j]);
				if(maps[i][j]==-1)
					maps[i][j]=inf;
			}
		floyd(n);
		MM.init(m,m);
		for(int i=1;i<=m;i++)
			scanf("%d%d%d",&s[i].place,&s[i].start,&s[i].during);
		for(int i=1;i<=m;i++)
		{
			for(int j=1;j<=m;j++)
			{
				if(i!=j&&maps[s[i].place][s[j].place]!=-1)
				{
					if(s[i].start+s[i].during+maps[s[i].place][s[j].place]<=s[j].start)
					{
						MM.g[i].push_back(j);
					}
				}
			}
		}
		printf("%d\n",m-MM.solve());
	}
    return 0;
}