PKU2060 Taxi Cab Scheme - 最小路径覆盖

题目描述：

出租车公司接到N（N<500）个订单，每个订单的描述是：t时刻需要从[x,y]坐标出发去[tx,ty]坐标。出租车行驶需要的时间是|x-tx|+|y-ty|。出租车公司希望派出最少的车完成所有的订单。派出一个出租车，它可以去接任何一个订单；若一个出租车在完成了一个订单之后，能在下一个订单的时刻之前赶到出发地点，那么它可以继续接下一个订单。输出最少需要多少出租车。

分析：

将一个订单作为有向无环图的一个节点，若订单i完成后能够赶到订单j，那么在有向图i到j之间连有向边。然后求这个有向图的最小路径覆盖。

最小路径覆盖问题可以转换为二分图最大匹配问题：若节点i到节点j有连边，那么在二部图的左边节点i到右边节点j连一条边。（事实上就是原来的有向图，只是转换为二分图来考虑）最小路径覆盖数=节点数N-最大匹配数。

关于这个题注意的是，不仅仅是前一个点的终点等于后一个点的起点才能连边，只要出租车在前一个目的地能开车赶到后一个起点也是可以连边的。

 

1. /*
2. PKU2060 Taxi Cab Scheme
3. */
5. #include <stdio.h>
6. #include <memory.h>
8. #define clr(a) memset(a,0,sizeof(a))
9. #define ABS(a) ((a)>0?(a):-(a))
10. #define N 1005
12. typedef struct{
13.     int h,m;
14. }Time;
15. typedef struct{
16.     int x,y,tx,ty;
17. }Ride;
19. int find(int i,int m,int g[][N],int mat[],int tmat[]){
20.     int v,j;
21.     for(j=0;j<m;j++) if(g[i][j]&&tmat[j]==0){
22.         tmat[j]=1;v=mat[j];mat[j]=i;
23.         if(v==-1||find(v,m,g,mat,tmat)) return 1;
24.         mat[j]=v;
25.     }
26.     return 0;
27. }
29. int match(int g[][N],int n,int m,int mat[]){
30.     int i,j,v,k=0;
31.     int tmat[N];
32.     for(i=0;i<m;i++) mat[i]=-1;
33.     for(i=0;i<n;i++){clr(tmat); k+=find(i,m,g,mat,tmat);}
34.     return k;
35. }
37. //最小路径覆盖
38. int PathCover(int g[][N],int n){
39.     int mat[N];
40.     return n-match(g,n,n,mat);
41. }
43. int Dist(int x,int y,int tx,int ty){
44.     return ABS(x-tx)+ABS(y-ty);
45. }
47. int Early(Time a,Ride p,Time b,Ride q){
48.     int dist=0;
49.     dist+=Dist(p.x,p.y,p.tx,p.ty);
50.     dist+=Dist(p.tx,p.ty,q.x,q.y);
52.     a.m+=dist;
53.     a.h+=a.m/60;
54.     a.m%=60;
55.     if(a.h<b.h) return 1;
56.     if(a.h>b.h) return 0;
57.     if(a.m<b.m) return 1;
58.     return 0;
59. }
61. Time t[N];
62. Ride r[N];
63. int a[N][N];
65. int main()
66. {
67.     int i,j,k;
68.     int T;
69.     int n;
70.     
71.     scanf("%d",&T);
72.     while(T--){
73.         //input
74.         scanf("%d",&n);
75.         for(k=0;k<n;k++){
76.             scanf("%d:%d",&t[k].h,&t[k].m);
77.             scanf("%d%d%d%d",&r[k].x,&r[k].y,&r[k].tx,&r[k].ty);
78.         }
79.         //create map
80.         clr(a);
81.         for(i=0;i<n;i++)
82.         for(j=i+1;j<n;j++)
83.             if( Early(t[i],r[i],t[j],r[j]) )
84.                 a[i][j]=1;
85.         //output
86.         printf("%d/n",PathCover(a,n));
87.     }
89.     return 0;
90. }