洛谷$P3324\ [SDOI2015]$星际战争 网络流+二分

正解:网络流+二分

解题报告:

传送门$QwQ$

其实我第一反应是费用流来着,,,但是仔细想了下发现我不会实现各个武器之间独立同时?而且攻击是连续的答案可能是小数嘛$QwQ$.

所以显然不是递推就二分呗$QwQ$.然后依然是因为小数的约束,所以最后选择的二分$QwQ$

$umm$然后就傻逼题了呗?

$check$的话直接武器建一排点机器人建一排点,$S$和武器$i$连$B_{i}\cdot time$,机器人$j$和$T$连$A_{i}$,然后武器和能攻击的机器人之间连$inf$,欧克做完辣!

$over$,感觉也没什么要注意的细节,那就放下代码就成$QwQ$

#include<bits/stdc++.h>
using namespace std;
#define il inline
#define lf double
#define gc getchar()
#define t(i) edge[i].to
#define n(i) edge[i].nxt
#define w(i) edge[i].wei
#define ri register int
#define rb register bool
#define rc register char
#define rp(i,x,y) for(ri i=x;i<=y;++i)
#define e(i,x) for(ri i=head[x];~i;i=n(i))

const int N=300000+10,M=60;
const lf eps=1e-5,inf=1e9;
int S,T,ed_cnt=-1,head[N],cur[N],a[N],b[N],dep[N],m,n;
bool d[M][M];
lf l,r=1e6,sum;
struct ed{int to,nxt;lf wei;}edge[N<<2];

il int read()
{
    rc ch=gc;ri x=0;rb y=1;
    while(ch!='-' && (ch>'9' || ch<'0'))ch=gc;
    if(ch=='-')ch=gc,y=0;
    while(ch>='0' && ch<='9')x=(x<<1)+(x<<3)+(ch^'0'),ch=gc;
    return y?x:-x;
}
il void ad(ri x,ri y,lf z)
{edge[++ed_cnt]=(ed){x,head[y],z};head[y]=ed_cnt;edge[++ed_cnt]=(ed){y,head[x],(lf)0};head[x]=ed_cnt;}
il bool bfs()
{
    queue<int>Q;Q.push(S);memset(dep,0,sizeof(dep));dep[S]=1;
    while(!Q.empty()){ri nw=Q.front();Q.pop();e(i,nw)if(w(i)-eps>0 && !dep[t(i)])dep[t(i)]=dep[nw]+1,Q.push(t(i));}
    return dep[T];
}
il lf dfs(ri nw,lf flow)
{
    if(nw==T || !flow)return flow;lf ret=0;
    for(ri &i=cur[nw];~i;i=n(i))
        if(w(i) && dep[t(i)]==dep[nw]+1)
        {lf tmp=dfs(t(i),min(flow,w(i)));ret+=tmp,w(i)-=tmp;w(i^1)+=tmp,flow-=tmp;}
    return ret;
}
il lf dinic(){lf ret=0;while(bfs()){rp(i,S,T)cur[i]=head[i];while(lf d=dfs(S,inf))ret+=d;}return ret;}
il bool check(lf tim)
{
    ed_cnt=-1;memset(head,-1,sizeof(head));
    rp(i,1,m)rp(j,1,n)if(d[i][j])ad(j+m,i,inf);
    rp(i,1,n)ad(T,i+m,a[i]);rp(i,1,m)ad(i,S,b[i]*tim);
    //printf("dinic=%.6lf sum=%.1lf\n",dinic(),sum);
    return dinic()<sum;
}

int main()
{
    //freopen("3324.in","r",stdin);freopen("3324.out","w",stdout);
    n=read();m=read();S=0;T=n+m+1;rp(i,1,n)sum+=a[i]=read();rp(i,1,m)b[i]=read();rp(i,1,m)rp(j,1,n)d[i][j]=read();
    while(r-l>eps){lf mid=(l+r)/2;if(check(mid))l=mid;else r=mid;}printf("%.6lf\n",r);
    //printf("%d",check(82));
    //check(82);
    return 0;
}

View Code