码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ4657 : tower

时间:2016-07-31 20:48:21      阅读:181      评论:0      收藏:0      [点我收藏+]

标签:

显然只有横向和纵向的两个炮塔才有可能冲突。

考虑最小割,将每个炮塔所有能攻击到的位置建点,相邻之间连无穷的边,表示前缀和关系,即选了一个点,就必须要选所有比它近的点。

属于横向炮塔的点向$S$连边,容量为前缀最大值的差值;属于纵向炮塔的点向$T$连边,容量为前缀最大值的差值。

对于一个交点,则在两个点之间连无穷边,表示必须舍弃其中一个。

答案$=$总收益$-$最小割。

 

#include<cstdio>
#include<algorithm>
using namespace std;
const int N=126000,inf=~0U>>2;
int dx[4]={0,0,1,-1},dy[4]={1,-1,0,0};
struct E{int t,f;E*nxt,*pair;}*g[N],*d[N],pool[1000000],*cur=pool;
int n,m,i,j,a[55][55],b[55][55],cnt,S,T,h[N],gap[N],ans;
inline void add(int s,int t,int f){
  if(!f)return;
  if(f<inf)ans+=f;
  E*p=cur++;p->t=t;p->f=f;p->nxt=g[s];g[s]=p;
  p=cur++;p->t=s;p->f=0;p->nxt=g[t];g[t]=p;
  g[s]->pair=g[t];g[t]->pair=g[s];
}
int sap(int v,int flow){
  if(v==T)return flow;
  int rec=0;
  for(E*p=d[v];p;p=p->nxt)if(h[v]==h[p->t]+1&&p->f){
    int ret=sap(p->t,min(flow-rec,p->f));
    p->f-=ret;p->pair->f+=ret;d[v]=p;
    if((rec+=ret)==flow)return flow;
  }
  if(!(--gap[h[v]]))h[S]=T;
  gap[++h[v]]++;d[v]=g[v];
  return rec;
}
inline void tag(int x,int y,int k){
  int mx=0,pre=0,now;
  while(1){
    x+=dx[k],y+=dy[k];
    if(x<1||x>n||y<1||y>m)return;
    cnt++;
    if(pre){
      if(k<2)add(cnt,pre,inf);
      else add(pre,cnt,inf);
    }
    if(k<2)add(S,cnt,max(mx,a[x][y])-mx);else add(cnt,T,max(mx,a[x][y])-mx);
    mx=max(mx,a[x][y]);
    pre=cnt;
    if(b[x][y]){
      if(k<2)add(cnt,b[x][y],inf);
      else add(b[x][y],cnt,inf);
    }else b[x][y]=cnt;
  }
}
int main(){
  scanf("%d%d",&n,&m);
  S=n*m*max(n,m)+1;T=S+1;
  for(i=1;i<=n;i++)for(j=1;j<=m;j++)scanf("%d",&a[i][j]);
  for(i=1;i<=n;i++)for(j=1;j<=m;j++)if(a[i][j]<0)tag(i,j,a[i][j]+4);
  for(gap[0]=T,i=1;i<=T;i++)d[i]=g[i];
  while(h[S]<T)ans-=sap(S,inf);
  return printf("%d",ans),0;
}

  

BZOJ4657 : tower

标签:

原文地址:http://www.cnblogs.com/clrs97/p/5723755.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!