bzoj 5120 [2017国家集训队测试]无限之环——网络流

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

旋转的话相当于去掉一个插头、新增一个插头，所以在这两个插头之间连边并带上费用即可。

网格图可以黑白染色，转化为相邻格子间插头的匹配问题。

注意：

1.黑白染色不是移动一格就 fx = ! fx ；每换一行，开头位置的颜色应该和上一行的开头不一样！不然有偶数列的话自己原来写的那个染色就崩了；

2. L 形的判断不是 d&(d>>1) 判断是否有两个相邻的1，如果是第一个位置和最后一个位置是1的话（9）就判断不出来了！可以判断 d != 5 && d != 10 ；

3.对于无解的判断，比较好的是判断插头个数的是不是流量的两倍。

然后就能很慢地 A 了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=,K=,M=K<<,INF=N;//<<4:*6 for trans *2 for between
int n,m,t,bh[K][K],hd[N],xnt=,to[M],nxt[M],cap[M],w[M],bin[];//=1!
int dis[N],pre[N],info[N],ans,val,flow; bool ins[N];
queue<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;
}
int Mn(int a,int b){return a<b?a:b;}
void upd(int &x){x>=?x-=:;}
void add(int x,int y,int z)
{
  if(!x||y==t)val++;//
  to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;cap[xnt]=;w[xnt]=z;
  to[++xnt]=x;nxt[xnt]=hd[y];hd[y]=xnt;cap[xnt]=;w[xnt]=-z;
}
void init()
{
  int bin[];
  bin[]=;for(int i=;i<;i++)bin[i]=bin[i-]<<;
  bool Fx=,fx=;
  for(int i=;i<=n;i++,Fx=!Fx)/////Fx!!!!!
    for(int j=,fx=Fx;j<=m;j++,fx=!fx)
      {
    int d=rdn(),ct=;

    if(!fx){for(int k=;k<;k++)if(d&bin[k])ct++,add(,k+bh[i][j],);}
    else {for(int k=;k<;k++)if(d&bin[k])ct++,add(k+bh[i][j],t,);}
    if(ct==||ct==)///ct==1:not a chain but 1,1,2
      {
        int x;//caution for ct=1|3 the diff
        if(ct==) {for(int k=;k<;k++)if(d&bin[k]){x=k;break;}}
        else {for(int k=;k<;k++)if(!(d&bin[k])){x=k;break;}}
        int y=x^,cr=x+bh[i][j];
        for(int k=;k<;k++)
          if(k!=x)(fx^(ct==))?add(k+bh[i][j],cr,+(k==y)):add(cr,k+bh[i][j],+(k==y));
      }
    //    else if(ct==2&&(d&(d>>1)))//L
    else if(ct==&&(d!=&&d!=))/////caution 9!!!
      {
        for(int k=;k<;k++)
          if(d&bin[k])fx?add((k^)+bh[i][j],k+bh[i][j],):add(k+bh[i][j],(k^)+bh[i][j],);
      }
    if(!fx)
      {
        if(i>)add(bh[i][j],+bh[i-][j],);
        if(j<m)add(+bh[i][j],+bh[i][j+],);
        if(i<n)add(+bh[i][j],bh[i+][j],);
        if(j>)add(+bh[i][j],+bh[i][j-],);
      }
      }
}
bool spfa()
{
  memset(dis,0x3f,sizeof dis);dis[]=;
  info[]=INF;info[t]=; q.push();ins[]=;
  while(q.size())
    {
      int k=q.front();q.pop();ins[k]=;
      for(int i=hd[k],v;i;i=nxt[i])
    if(cap[i]&&dis[v=to[i]]>dis[k]+w[i])
      {
        dis[v]=dis[k]+w[i];
        pre[v]=i;info[v]=Mn(info[k],cap[i]);
        if(!ins[v])q.push(v),ins[v]=;
      }
    }
  return info[t];
}
void ek()
{
  int s=info[t]; flow+=s; ans+=dis[t]*s;//dis[t]!
  for(int i=pre[t];i;i=pre[to[i^]])
    {cap[i]-=s;cap[i^]+=s;}
}
int main()
{
  n=rdn();m=rdn();t=;
  for(int i=;i<=n;i++)
    for(int j=;j<=m;j++,t+=)bh[i][j]=t;
  init();  if(val&){puts("-1");return ;}//
  while(spfa())ek();  printf("%d\n",flow<<==val?ans:-);///
  return ;
}

还学习了多路增广SPFA算法。大概适用于边权变化范围很小的那种吧。

因为边权相近，所以期望做一次SPFA之后能不仅更新出一条路，还能更多地更新，因为算出来的 dis[ ] 可能对很多条路都适用。

所以就在SPFA之后按照 dis[ ] 来 dfs ，找找所有合法的路径。在这道题上真的能快好多。

注意 dfs 要打 vis 标记，走过的就不再走了。因为网络里可能有环，这个又是按 dis 来走的，所以比如有很多0权边的话，就可能死循环。dinic 没有这个问题是因为它是按 dfn 来走的，dfn 弄出来的一定是 DAG 。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=,K=,M=K<<,INF=N;//<<4:*6 for trans *2 for between
int n,m,t,bh[K][K],hd[N],xnt=,cur[N],to[M],nxt[M],cap[M],w[M],bin[];//=1!
int dis[N],pre[N],info[N],ans,val,flow; bool ins[N];
queue<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;
}
int Mn(int a,int b){return a<b?a:b;}
void upd(int &x){x>=?x-=:;}
void add(int x,int y,int z)
{
  if(!x||y==t)val++;//
  to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;cap[xnt]=;w[xnt]=z;
  to[++xnt]=x;nxt[xnt]=hd[y];hd[y]=xnt;cap[xnt]=;w[xnt]=-z;
}
void init()
{
  int bin[];
  bin[]=;for(int i=;i<;i++)bin[i]=bin[i-]<<;
  bool Fx=,fx=;
  for(int i=;i<=n;i++,Fx=!Fx)/////Fx!!!!!
    for(int j=,fx=Fx;j<=m;j++,fx=!fx)
      {
    int d=rdn(),ct=;

    if(!fx){for(int k=;k<;k++)if(d&bin[k])ct++,add(,k+bh[i][j],);}
    else {for(int k=;k<;k++)if(d&bin[k])ct++,add(k+bh[i][j],t,);}
    if(ct==||ct==)///ct==1:not a chain but 1,1,2
      {
        int x;//caution for ct=1|3 the diff
        if(ct==) {for(int k=;k<;k++)if(d&bin[k]){x=k;break;}}
        else {for(int k=;k<;k++)if(!(d&bin[k])){x=k;break;}}
        int y=x^,cr=x+bh[i][j];
        for(int k=;k<;k++)
          if(k!=x)(fx^(ct==))?add(k+bh[i][j],cr,+(k==y)):add(cr,k+bh[i][j],+(k==y));
      }
    //    else if(ct==2&&(d&(d>>1)))//L
    else if(ct==&&(d!=&&d!=))/////caution 9!!!
      {
        for(int k=;k<;k++)
          if(d&bin[k])fx?add((k^)+bh[i][j],k+bh[i][j],):add(k+bh[i][j],(k^)+bh[i][j],);
      }
    if(!fx)
      {
        if(i>)add(bh[i][j],+bh[i-][j],);
        if(j<m)add(+bh[i][j],+bh[i][j+],);
        if(i<n)add(+bh[i][j],bh[i+][j],);
        if(j>)add(+bh[i][j],+bh[i][j-],);
      }
      }
}
bool spfa()
{
  memset(ins,,sizeof ins);
  memset(dis,0x3f,sizeof dis);dis[]=;
  info[]=INF;info[t]=; q.push();ins[]=;
  while(q.size())
    {
      int k=q.front();q.pop();ins[k]=;
      for(int i=hd[k],v;i;i=nxt[i])
    if(cap[i]&&dis[v=to[i]]>dis[k]+w[i])
      {
        dis[v]=dis[k]+w[i];
        pre[v]=i;info[v]=Mn(info[k],cap[i]);
        if(!ins[v])q.push(v),ins[v]=;
      }
    }
  return info[t];
}
int dfs(int cr)//for cr!=0 flow must be 1
{
  if(cr==t)return ; ins[cr]=;
  for(int& i=cur[cr],v,tmp;i;i=nxt[i])
    if(!ins[v=to[i]]&&cap[i]&&dis[v=to[i]]==dis[cr]+w[i])
      if(tmp=dfs(v)){ans+=w[i];cap[i]--;cap[i^]++;return ;}
  return ;
}
void MCMF()
{
  int tmp;
  while(spfa())
    {
      memcpy(cur,hd,sizeof hd);
      while(tmp=dfs())flow+=tmp;
    }
}
int main()
{
  n=rdn();m=rdn();t=;
  for(int i=;i<=n;i++)
    for(int j=;j<=m;j++,t+=)bh[i][j]=t;
  init();  if(val&){puts("-1");return ;}//
  MCMF();  printf("%d\n",flow<<==val?ans:-);///
  return ;
}