ZOJ3156_Taxi(二分圖/二分構圖)

ZOJ3156_Taxi(二分圖/二分構圖)

解題報告

題意：

n個人，m輛車，給出人和車的坐標，還有人的速度，求全部人都坐上車的最小時間。（一輛車只能做一個人）

思路：

原本以為在二分圖上求最小的時間就變成了求二分圖的最佳匹配，其實可以二分時間，滿足時間的人和車連線構圖，這樣二分的求最小時間。

可能寫挫了，我竟然這樣二分時間，，，sad

其實只要把所有的可能時間存下在二分時間更快。

#include 
#include 
#include 
#include 
#define N 1000
#define M 1000
#define inf 99999999
using namespace std;
int n,m,mmap[N][N],vis[N],pre[N];
struct point
{
    double x,y;
} g[N],h[N];
double _time[110][110];
double dis(point p1,point p2)
{
    return  sqrt((p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
int dfs(int x)
{
    for(int i=1; i<=m; i++)
    {
        if(!vis[i]&&mmap[x][i])
        {
            vis[i]=1;
            if(pre[i]==-1||dfs(pre[i]))
            {
                pre[i]=x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int i,j;
    double v;
    while(~scanf("%d%d",&n,&m))
    {
        memset(mmap,0,sizeof(mmap));
        for(i=1; i<=n; i++)
        {
            scanf("%lf%lf",&g[i].x,&g[i].y);
        }
        for(i=1; i<=m; i++)
        {
            scanf("%lf%lf",&h[i].x,&h[i].y);
        }
        scanf("%lf",&v);
        double tmax=-1;
        for(i=1; i<=n; i++)
        {
            for(j=1; j<=m; j++)
            {
                _time[i][j]=dis(g[i],h[j])/v;
                if(_time[i][j]>tmax)
                {
                    tmax=_time[i][j];
                }
            }
        }
        double l=0,r=tmax;
        double minn=0;
        while(l<=r)
        {
            double mid=(l+r)/2;
            for(i=1; i<=n; i++)
            {
                for(j=1; j<=m; j++)
                {
                    if(_time[i][j]<=mid)
                        mmap[i][j]=1;
                    else mmap[i][j]=0;
                }
            }
            int ans=0;
            memset(pre,-1,sizeof(pre));
            for(i=1; i<=n; i++)
            {
                memset(vis,0,sizeof(vis));
                ans+=dfs(i);
            }
            if(ans>=n)
            {
                minn=mid;
                r=mid-0.00001;
            }
            else
            {
                l=mid+0.00001;
            }
        }
        printf("%.2lf\n",minn);
    }
    return 0;
}