导弹防御塔 Page428 二分图多重匹配

导弹防御塔 Page428 二分图多重匹配

• 二分时间
• 对于每一个导弹，在限定的时间内可以打出多少弹药是固定的，根据距离可以判断出每一枚子弹是否能够打中敌人，然后连线
  我是把左边防御塔，右边敌人，也可以变成左边敌人，右边防御塔，由于是拆开来做的多重匹配，所以原理是一样的
• 需要注意这边的T1是秒数，一开始看成了分钟怎么做都错
• 题目范围给的不是很明确，所以开数组的时候开了很大，这时候就不能用memset了，会超时

代码:

const int maxn=2e6+7;
const int INF=0x3f3f3f3f;
const double eps=1e-7;
const ll INFF=1e18;
int n1,tot,head[maxn],linker[maxn],n,m;
bool vis[maxn];
double T1,T2,V;
void init()
{
    n1=tot=0;
    mem(head,-1);
    mem(linker,0);
}
struct node
{
    double x,y;
}p[maxn],pp[maxn];
struct
{
    int to,next;
}edge[maxn];
void add(int u,int v){edge[++tot].to=v;edge[tot].next=head[u];head[u]=tot;}
double dis(node a,node b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
bool dfs(int x)
{
    for (int i=head[x];i!=-1;i=edge[i].next)
    {
        int y=edge[i].to;
        if (vis[y])continue;
        vis[y]=true;
        if (!linker[y]||dfs(linker[y]))
        {
            linker[y]=x;
            return true;
        }
    }
    return false;
}
int count_(int n1)
{
    int ans=0;
    rep(i,1,n1)
    {
        rep(j,1,n1)vis[j]=false;
        if (dfs(i))ans++;
    }
    return ans;
}
bool solve(double TotalTime)
{
    init();
    rep(i,1,n)
    {
        for (double k1=1.0;;k1+=1.0)
        {
            if (k1*T1+(k1-1)*T2-TotalTime>eps)break;//刚打出子弹的时间就已经超过了总时间
            n1++;
            rep(j,1,m)if (TotalTime-(k1*T1+(k1-1)*T2+dis(p[i],pp[j])/V)>eps)add(n1,j);
        }
    }
    return count_(n1)==m;//最大匹配是m，说明敌人被全部杀完了
}
int main()
{
    scanf("%d%d%lf%lf%lf",&n,&m,&T1,&T2,&V);T1/=60;
    rep(i,1,m)scanf("%lf%lf",&pp[i].x,&pp[i].y);
    rep(i,1,n)scanf("%lf%lf",&p[i].x,&p[i].y);
    double l=0,r=100000;
    while(l+eps<r)
    {
        double mid=(l+r)/2;
        if (solve(mid))r=mid;
        else l=mid;
    }
    printf("%.6lf\n",l);
    return 0;
}