bzoj 1001 [BeiJing2006]狼抓兔子——最小割转最短路

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1001

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
const int N=2e6+,M=6e6+;
int n,m,hd[N],xnt,to[M],nxt[M],w[M];
ll dis[N];bool vis[N];
priority_queue<pair<ll,int> > q;
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return fx?ret:-ret;
}
void add(int x,int y,int z)
{
  to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;w[xnt]=z;
}
void dj()
{
  memset(dis,0x3f,sizeof dis);dis[]=;
  q.push(make_pair(,));
  while(q.size())
    {
      int k=q.top().second;q.pop();
      if(vis[k])continue;
      for(int i=hd[k],v;i;i=nxt[i])
    if(dis[v=to[i]]>dis[k]+w[i])
      dis[v]=dis[k]+w[i],q.push(make_pair(-dis[v],v));
    }
}
int main()
{
  n=rdn();m=rdn();int t=(n-)*(m-)*+,pl=m-,ml=pl<<;
  for(int i=,d,lj=;i<=n;i++,lj+=ml)
    {
      int x=lj+;
      for(int j=,y;j<m;j++,x++)
    {
      d=rdn();if(!lj)y=;else y=x-pl;
      add(x>t?t:x,y,d);add(y,x>t?t:x,d);//not change x,for y
    }
    }
  for(int i=,d,lj=m-;i<n;i++,lj+=ml)
    {
      int x=lj+;
      for(int j=,y;j<=m;j++,x++)
    {
      d=rdn();if(x==lj+)y=t;else y=x-m;
      if(x>lj+pl)x=;//change x is ok for its the last x
      add(x,y,d);add(y,x,d);
    }
    }
  for(int i=,d,lj=;i<n;i++,lj+=ml)
    {
      int x=lj+;
      for(int j=,y;j<m;j++,x++)
    {
      d=rdn();y=x+pl;
      add(x,y,d);add(y,x,d);
    }
    }
  dj();printf("%lld\n",dis[t]);
  return ;
}