AcWing374 导弹防御塔

标签：stdin   二分答案   代码   for   double   inf   return   pen   c++   

二分图

此题看书的时候觉得特别难，实际的代码却非常简单

现在分析是什么让代码如此简单的：

1. 首先预处理出第i个防御塔发射第j个导弹的时间（计算发射时间，不计冷却时间）

2. 二分答案，判断时间mid内能否解决问题

3. 利用vector建边，不用管标号的冲突，在二分图中十分方便（一般网络流可能就没办法了）

时间复杂度：\(O(n^4*log(T))\) 实际更快

#include<bits/stdc++.h>
using namespace std;

#define go(i,a,b) for(int i=a;i<=b;++i)
#define com(i,a,b) for(int i=a;i>=b;--i)
#define mem(a,b) memset(a,b,sizeof(a))

const int inf=0x3f3f3f3f,N=60,M=2510;
const double eps=1e-8;
typedef long long ll;

int n,m,v,match[M],tot;
bool vis[M];
double t1,t2;
vector<int>e[N];
struct node{
    double x,y;
}a[N],b[N];
struct new_node{
    int id;double t;
}c[M];

bool dfs(int u){
    for(int i=0;i<e[u].size();++i){
        int v=e[u][i];
        if(vis[v]) continue;
        vis[v]=1;
        if(!match[v]||dfs(match[v])){
            match[v]=u;
            return 1;
        }
    }
    return 0;
}

double dis(const node &x,const node &y){
    return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
}

bool pd(double mid){
    mem(match,0);
    go(i,1,m){
        e[i].clear();
        go(j,1,tot)//代表n*m颗导弹（不一定每个都发出，但最多这么多） 
            if(c[j].t+dis(a[i],b[c[j].id])/v<=mid)
                e[i].push_back(j);
    }
    go(i,1,m){
        mem(vis,0);
        if(!dfs(i)) return 0;
    }
    return 1;
}

signed main(){
    //freopen("input.txt","r",stdin);
    scanf("%d%d%lf%lf%d",&n,&m,&t1,&t2,&v);
    tot=n*m;
    t1/=60;
    go(i,1,m)
        scanf("%lf%lf",&a[i].x,&a[i].y);
    go(i,1,n)
        scanf("%lf%lf",&b[i].x,&b[i].y);
    int k;
    go(i,1,m)
        go(j,1,n){
            k=(i-1)*n+j;
            c[k].id=j;
            c[k].t=(i-1)*(t1+t2)+t1;
        }
    double l=t1,r=100000;
    while(r-l>eps){
        double mid=(l+r)*0.5;
        if(pd(mid)) r=mid;
        else l=mid;
    }
    printf("%.6lf",r);
    return 0;
}

AcWing374 导弹防御塔

标签：stdin   二分答案   代码   for   double   inf   return   pen   c++   

原文地址：https://www.cnblogs.com/White-star/p/11526325.html