codevs 1227 方格取数 网络流

题目描述 Description

给出一个n*n的矩阵,每一格有一个非负整数Aij,(Aij <= 1000)现在从(1,1)出发,可以往右或者往下走,最后到达(n,n),每达到一格,把该格子的数取出来,该格子的数就变成0,这样一共走K次,现在要求K次所达到的方格的数的和最大

输入描述 InputDescription

第一行两个数n,k(1<=n<=50, 0<=k<=10)

接下来n行,每行n个数,分别表示矩阵的每个格子的数

输出描述 OutputDescription

一个数,为最大和

样例输入 SampleInput

3 1

1 2 3

0 2 1

1 4 2

样例输出 SampleOutput

11

数据范围及提示 Data Size &Hint

1<=n<=50, 0<=k<=10

    嗯题目说的很清楚,经过分析我们会知道,如果经过了一个格子已经被取过,那么此时继续走这个格子肯定不是最优方案。那么我们可以由此入手。总共有k次机会,那么也就是说这个网络的流量就是k,那么我们把第一行的所有节点与S'连接,且最大流量均为1再将S'与S相连,S -》S'的最大流量即为k。然后将最后一行的所有节点与E相连,最大流量为1。这里的通过每个点的代价即为这个格子的值,然后我们需要把有代价(代价为a)的节点拆成两个无代价点,两个点之间用一条代价为 a的边。当我们构图时需要注意,除S'到S的最大流量为k外,其他边的最大流量始终为1(除反平行边)。然后再建普通边 格子【i】【j】到格子【i】【j-1】和到
格子【i-1】【j-1】的边。

    我的程序里多加了一个点mid,其实没什么用。
  #include
#include
#include
#include
#include
using namespace std;
int now=1,start=1,end,n,m;
intnext[100001],head[100001],other[100001],flow[100001],cost[100001],cap[100001];
inline void add_side(int start,int end,int floww,int capp,intcostt)
{
    next[++now]=head[start];
    other[now]=end;
    head[start]=now;
    flow [now] = floww;
    cap [now] = capp;
    cost[now]=costt;
    return ;
}

struct point
{
int mid,in,out;
};
point pointnumber[100][100];


inline point point_number(int x,int y)
  {
  int k=0;
 
      k=n*(x-1)+y;
  point res;
res.in=(k-1)*3+2;
res.out=k*3+1;
res.mid=res.out-1;
return res;    
  }

inline void point_add(int x,int y)
  {
  int j;
  scanf("%d",&j);
      add_side(point_number(x,y).in,point_number(x,y).out,0,1,j),
  add_side(point_number(x,y).out,point_number(x,y).in,0,0,-j);
  if(x!=1)
    if(y!=1)
    add_side(point_number(x-1,y).out,point_number(x,y).in,0,1,0),
    add_side(point_number(x,y).in,point_number(x-1,y).out,0,0,0),
add_side(point_number(x,y-1).out,point_number(x,y).in,0,1,0),
add_side(point_number(x,y).in,point_number(x,y-1).out,0,0,0);
     
  else
    add_side(point_number(x-1,y).out,point_number(x,y).in,0,1,0),
    add_side(point_number(x,y).in,point_number(x-1,y).out,0,0,0);   
    else
    if(y!=1)
    add_side(0,point_number(x,y).in,0,1,0),
      add_side(point_number(x,y).in,0,0,0,0),
          add_side(point_number(x,y-1).out,point_number(x,y).in,0,1,0),
add_side(point_number(x,y).in,point_number(x,y-1).out,0,0,0);
    else
          add_side(0,point_number(x,y).in,0,1,0),
      add_side(point_number(x,y).in,0,0,0,0);
       
    if(x==n)
        add_side(point_number(x,y).out,end,0,1,0),
  add_side(end,point_number(x,y).out,0,0,0);
   
 
  }

void read()
{
int k;
      scanf("%d%d",&n,&k);
      end=1;
      add_side(1,0,0,k,0),add_side(0,1,0,0,0);
      for(int i=1;i<=n;i++)
      for(int j=1;j<=n;j++)
        end+=3;
        end++;
      for(int i=1;i<=n;i++)
      for(int j=1;j<=n;j++)
        point_add(i,j);
}
int res;
int line[100001],pre[100001],dis[100001];
bool visited[100001];


inline bool spfa()
  {
  res=0;
  memset(line,0,sizeof(line));
  memset(visited,false,sizeof(visited));
  memset(dis,0,sizeof(dis));
  memset(pre,-1,sizeof(pre));
line[1]=start;
visited[start]=true;
  int l,r,left,right,j;
  l=r=1;
  while(l<=r)
    {
    left=l%end;
    j=head[line[left]];
        while(j!=0)
        {
        if(flow[j]dis[other[j]]||
              (pre[other[j]]==-1&&dis[line[left]]+cost[j]==dis[other[j]])))
          {
          dis[other[j]]=dis[line[left]]+cost[j],
pre[other[j]]=line[left];
          if(!visited[other[j]])
    line[(++r)%end]=other[j],visited[other[j]]=true;
          }
        j=next[j];
        }
    if(line[left]!=1)visited[line[left]]=false;
    l++;
    }
    res=dis[end];
    if(res==0)
      returnfalse;
    int i=end;
    do
    {
    j=head[pre[i]];
    while(other[j]!=i)
    {
    j=next[j];
    }
    flow[j]+=1;
flow[j^1]-=1; 
    i=pre[i];
 
    }
    while(i!=1);
return res!=0;
  }
int main()
{
read();
int ans=0;
while(spfa())
ans+=res;
printf("%d\n",ans);
return 0;
}






你可能感兴趣的:(codevs 1227 方格取数 网络流)