「Poetize3」导弹防御塔

描述 Description

Freda控制着N座可以发射导弹的防御塔。每座塔都有足够数量的导弹，但是每座塔每次只能发射一枚。在发射导弹时，导弹需要T1秒才能从防御塔中射出，而在发射导弹后，发射这枚导弹的防御塔需要T2分钟来冷却。
所有导弹都有相同的匀速飞行速度V，并且会沿着距离最短的路径去打击目标。计算防御塔到目标的距离Distance时，你只需要计算水平距离，而忽略导弹飞行的高度。导弹在空中飞行的时间就是 (Distance/V) 分钟，导弹到达目标后可以立即将它击毁。
现在，给出N座导弹防御塔的坐标，M个入侵者的坐标，T1、T2和V，你需要求出至少要多少分钟才能击退所有的入侵者。

题解：
不太好想。

二分时间+二分图匹配

因为每个塔打每个敌人的时间点不同，能否达到也不同，而我们发现这个时间点只取决于这个该塔发射的第几次导弹。

所以我们把每个点拆成m个点，分别表示这个塔的第 i 发导弹，然后就可以根据当前时间+到达时间是不是小于二分的这个时间来建图。

用匈牙利求最大匹配，如果为m说明可行，缩小上界，否则缩小下界。

打的时候出现了出现了几个错误：

1）没有清空v数组

2）题目输入坑爹，t1的单位是s，其他的单位都是分钟。。。

3）。。。

现在还没有A掉，不过实在受不鸟了。。。

代码；

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 200

 #define maxm 200

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007
 #define sqr(x) (x)*(x)
 #define num(x,y) (x-1)*(m)+y

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 int n,m,v[maxn*maxn],p[maxn*maxn];
 double t1,t2,vv,t[maxn][maxn],d[maxn];
 bool f[maxn*maxn][maxn];
 struct rec{double x,y;}a[maxn],b[maxn];
 inline bool find(int x,int y)
 {
     for1(i,n*m)
      if(f[x][i]&&v[i]!=y)
      {
          v[i]=y;
          if(p[i]==||find(p[i],y))
          {
              p[i]=x;
              return ;
          }
      }
     return ;
 }
 bool check(double x)
 {
     int ans=;//cout<<x<<endl;
     memset(f,,sizeof(f));memset(p,,sizeof(p));memset(v,,sizeof(v));
     for1(i,m)
      for1(j,n)
       for1(k,m)
         if(d[k]+t[j][i]<=x)f[i][num(j,k)]=;
     for1(i,m)if(find(i,i))ans++;
     //cout<<ans<<endl;
     return ans==m;
 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();m=read();t1=read();t2=read();vv=read();t1/=;
     for1(i,m)b[i].x=read(),b[i].y=read();
     for1(i,n)a[i].x=read(),a[i].y=read();
     for1(i,n)
      for1(j,m)
       t[i][j]=sqrt((sqr(a[i].x-b[j].x)+sqr(a[i].y-b[j].y)))/vv;
     for1(i,m)d[i]=t1*i+t2*(i-);
     double l=,r=inf,mid;
     while(r-l>1e-)
     {
         mid=(l+r)/;
         //printf("%lf %lf %lf\n",l,mid,r);
         if(check(mid))r=mid;else l=mid;
     }
     printf("%.6lf",l);

     return ;

 }

数组开小查了一天啊！！！！！！！！！！！！！！！！！！！！！！！！！！！1

怎么这么容易反sb错误。。。。。。。

我是不是没救了TAT

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 200

 #define maxm 200

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007
 #define sqr(x) (x)*(x)
 #define num(x,y) (x-1)*(m)+y

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 int n,m,v[maxn*maxn],p[maxn*maxn];
 double t1,t2,vv,t[maxn][maxn],d[maxn];
 bool f[maxn][maxn*maxn];
 struct rec{double x,y;}a[maxn],b[maxn];
 inline bool find(int x,int y)
 {
     for1(i,n*m)
      if(f[x][i]&&v[i]!=y)
      {
          v[i]=y;
          if(p[i]==||find(p[i],y))
          {
              p[i]=x;
              return ;
          }
      }
     return ;
 }
 bool check(double x)
 {
     int ans=;
     memset(f,,sizeof(f));memset(p,,sizeof(p));memset(v,,sizeof(v));
     for1(i,m)
      for1(j,n)
       for1(k,m)
         if(d[k]+t[j][i]<=x)f[i][num(j,k)]=;
     for1(i,m)if(find(i,i))ans++;
     return ans==m;
 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();m=read();t1=read();t2=read();vv=read();t1/=;
     for1(i,m)b[i].x=read(),b[i].y=read();
     for1(i,n)a[i].x=read(),a[i].y=read();
     for1(i,n)
      for1(j,m)
       t[i][j]=sqrt((sqr(a[i].x-b[j].x)+sqr(a[i].y-b[j].y)))/vv;
     for1(i,m)d[i]=(i-)*(t1+t2)+t1;
     double l=,r=inf,mid;
     while(r-l>1e-)
     {
         mid=(l+r)/;
         if(check(mid))r=mid;else l=mid;
     }
     printf("%.6lf",l);

     return ;

 }