hdu 2389 Rain on your Parade 二分图匹配+HK算法

原创 2015年11月17日 21:45:43
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
const int MAXN=3000+5;
const int INF=1<<28;
int g[MAXN][MAXN],Mx[MAXN],My[MAXN],Nx,Ny;
int dx[MAXN],dy[MAXN],dis;
bool vst[MAXN];
bool searchP()
{
    queue<int>Q;
    dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=Nx;i++)
        if(Mx[i]==-1)
        {
            Q.push(i);
            dx[i]=0;
        }
    while(!Q.empty())
    {
        int u=Q.front();
        Q.pop();
        if(dx[u]>dis)  break;
        for(int v=1;v<=Ny;v++)
            if(g[u][v]&&dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(My[v]==-1)  dis=dy[v];
                else
                {
                    dx[My[v]]=dy[v]+1;
                    Q.push(My[v]);
                }
            }
    }
    return dis!=INF;
}
bool DFS(int u)
{
    for(int v=1;v<=Ny;v++)
       if(!vst[v]&&g[u][v]&&dy[v]==dx[u]+1)
       {
           vst[v]=1;
           if(My[v]!=-1&&dy[v]==dis) continue;
           if(My[v]==-1||DFS(My[v]))
           {
               My[v]=u;
               Mx[u]=v;
               return 1;
           }
       }
    return 0;
}
int MaxMatch()
{
    int res=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(My));
    while(searchP())
    {
        memset(vst,0,sizeof(vst));
        for(int i=1;i<=Nx;i++)
          if(Mx[i]==-1&&DFS(i))  res++;
    }
    return res;
}


int main()
{
    int _,n,m,i,j,cas=0;
    double t,tp,x1[MAXN],y1[MAXN],t1[MAXN],x[MAXN],y[MAXN];
    scanf("%d",&_);
    while(_--)
    {
        memset(g,0,sizeof(g));
        scanf("%lf",&t);
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            scanf("%lf%lf%lf",&x1[i],&y1[i],&t1[i]);

        scanf("%d",&m);
        for(i=1;i<=m;i++)
            scanf("%lf%lf",&x[i],&y[i]);

        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
        {
            tp=(x1[i]-x[j])*(x1[i]-x[j])+(y1[i]-y[j])*(y1[i]-y[j]);
            tp=sqrt(tp);
            if(tp<=t1[i]*t)
                g[i][j]=1;
        }
        Nx=n;
        Ny=m;
        printf("Scenario #%d:\n%d\n\n",++cas,MaxMatch());
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。http://www.xiang578.com

HDU 2389 Rain on your Parade(HK算法)

看这数据量,显然匈牙利算法要超时,就拿了HK的板子上了。#include #include #include #include using namespace std; const int maxn=...
  • qq_33901573
  • qq_33901573
  • 2017年01月06日 13:35
  • 67

HDU 2389 Rain on your Parade (HK二分图)

大意:给定时间t,并且n个人的坐标,然后有伞的坐标求最少有多少人不被淋湿。hk算法模板。#include #include #include #include #include #include #i...
  • Grit_ICPC
  • Grit_ICPC
  • 2016年04月24日 14:08
  • 187

hdu 2389 Rain on your Parade (HK)

http://acm.hdu.edu.cn/showproblem.php?pid=2389 Rain on your Parade Problem Description Yo...
  • w144215160044
  • w144215160044
  • 2015年08月15日 15:02
  • 233

HDU2389 Rain on your Parade(二分图匹配 ,HK算法 )

匈牙利算法果断超时,只有用HC模板测试一下 #include #include #include #include #include #include #include #include #inc...
  • u013167299
  • u013167299
  • 2015年08月24日 20:48
  • 897

hdu 2389 Rain on your Parade 二分匹配 HK算法

题目链接 匈牙利算法
  • u011742541
  • u011742541
  • 2014年07月17日 20:43
  • 280

hdu 2389 Rain on your Parade(二分最大匹配HK算法)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2389
  • u013903504
  • u013903504
  • 2014年08月20日 16:41
  • 429

【二分匹配】 HDU 2389 Rain on your Parade HK算法

注意建图 原来不是曼哈顿距离。。。 #include #include #include #include #include #include #include #include ...
  • u012749539
  • u012749539
  • 2014年09月15日 22:26
  • 538

二分图最大匹配(HK)

二分图最大匹配经典的算法就是匈牙利算法,但是本文并不是讲述匈牙利算法,而是说一个时间复杂度更为优的HK算法。...
  • qq_35776579
  • qq_35776579
  • 2017年02月09日 11:10
  • 651

Rain on your Parade HDU - 2389

HK算法,匈牙利算法超时 题目描述:t时间后会有一场大雨,有m个人和n个雨伞,求在t时间内,使得最多的人不被雨淋。求不被雨淋的人数。 解题分析:显然是匹配问题,第一次写用匈牙利算法超时了,时间复杂...
  • Mrhanice
  • Mrhanice
  • 2017年07月18日 16:36
  • 116

hdoj 2389 Rain on your Parade 【HK】

Rain on your Parade Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 655350/165535 K (Java/Ot...
  • chenzhenyu123456
  • chenzhenyu123456
  • 2015年10月17日 11:00
  • 231
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hdu 2389 Rain on your Parade 二分图匹配+HK算法
举报原因:
原因补充:

(最多只允许输入30个字)