首页 诗词 字典 板报 句子 名言 友答 励志 学校 网站地图
当前位置: 首页 > 教程频道 > 开发语言 > 编程 >

[2分+匈牙利]zoj 3460:Missile

2012-09-04 
[二分+匈牙利]zoj 3460:Missile大致题意:? ? 用n个导弹发射塔攻击m个目标。每个发射架在某个时刻只能为一颗

[二分+匈牙利]zoj 3460:Missile

大致题意:
? ? 用n个导弹发射塔攻击m个目标。每个发射架在某个时刻只能为一颗导弹服务,发射一颗导弹需要准备t1的时间,一颗导弹从发射到击中目标的时间与目标到发射架的距离有关。每颗导弹发射完成之后发射架需要t2的时间进入下个发射流程。现在问最少需要多少时间可以击毁所有m个目标。

?

大致思路:
? ? 二分枚举这个最大时间的最小值,每次按照这个枚举的时间构出二分图,求最大匹配来判定枚举值是否符合要求。

?

#include<iostream>#include<cstring>#include<cmath>#include<cstdio>using namespace std;const int inf=1<<30;const double eps=1e-8;int N,M;double tt1,tt2,v;int mis[100][2],tar[100][2];double dis[3000][60];double getdis(int a,int b){    double tmp=(mis[a][0]-tar[b][0])*(mis[a][0]-tar[b][0])+(mis[a][1]-tar[b][1])*(mis[a][1]-tar[b][1]);    double res=sqrt(tmp);    return res;}void getinit(){    int i,j,k;    for(i=1;i<=N;i++)  // 枚举第i个发射塔    {        for(j=1;j<=M;j++)  //第i个发射塔的第j次发射        {            for(k=1;k<=M;k++)  //第j次发射,攻击第k个目标            {                dis[(i-1)*M+j][k]=(j-1)*tt2+getdis(i,k)/v+j*tt1;            }        }    }}int map[3000][60];bool vis[3000];int linkk[3000];int dfs(int s){    for(int i=1;i<=M;i++){        if(!vis[i]&&map[s][i]){            vis[i]=1;            if(linkk[i]==-1||dfs(linkk[i])){                linkk[i]=s;                return 1;            }        }    }    return 0;}bool check(double mid){//    cout<<"check"<<mid<<endl;    int i,j,res=0;    memset(map,0,sizeof(map));    for(i=1;i<=N*M;i++)    {        for(j=1;j<=M;j++)        {            if(dis[i][j]<=mid)            {                map[i][j]=1;            }        }    }    memset(linkk,-1,sizeof(linkk));    for(i=1;i<=N*M;i++){        memset(vis,0,sizeof(vis));        if(dfs(i))res++;    }    if(res==M)return 1;    return 0;}int main(){    int i;    while(scanf("%d%d%lf%lf%lf",&N,&M,&tt1,&tt2,&v)!=EOF)    {        tt1/=60.0;        for(i=1;i<=M;i++)        {            scanf("%d%d",&tar[i][0],&tar[i][1]);        }        for(i=1;i<=N;i++)        {            scanf("%d%d",&mis[i][0],&mis[i][1]);        }        getinit();        double left=0,right=200000000,mid;        while(right-left>=eps)        {            mid=(left+right)/2.0;            if(check(mid))right=mid;            else left= mid;        }        printf("%.6f\n",left);    }    return 0;}
?

热点排行