吐槽 isap/ 转dinic

  

       经典问题：方格取数，二分图的最大权独立集，直接把我的isap卡的天崩地裂。。。。。。

      

        >.<!话说当时写现在用的程序是以前看了某某的博客，加上YY了一顿，no 当前弧优化， 单路增广，递归版。。。。。。仅仅加了个gap优化>.<

        就这样还是挤过了秋哥的“最大流强数据”，还踩掉了hyc的多路增广dinic = 。 =！结果这次一个普通的二分图就把我卡了。。。。。。

        

         以前那个丑程序就不贴了=。=！

 

         热烈祝贺 从递归版单路增广无当前弧isap  转变为 非递归版多路增广加当前弧dinic（转变真大）， 感谢lzn 提供程序，以及lyp，hyc 提供不同的代码对比长度。

  

         3个过程一共37 行，无奇葩缩行，很好理解。 

     速度快，空间小，代码短，真可谓是居家旅行之必备呀.

# include <cstdlib>
# include <cstdio>
# include <cmath>
#include <ctime>
#include <cstring>
using namespace std;

const int maxn = 110, maxm = 105000, oo=1073741819;
int top, next[maxm*2], point[maxm*2], linke[maxm*2], wis[maxm*2], a[maxn][maxn];
int tot, p, tab[maxm], prev[maxm], shift[maxm], flow, s, t, maxflow, n, m;
int mx[10], my[10];
int l, r, que[maxm], step[maxm], adm[maxm];
bool flag;
 
inline int min(int x, int y){return x<y?x:y;};
void link(int x, int y, int z)
{
	++top; next[top] = linke[x]; linke[x] = top; point[top] = y; wis[top] = z;
	++top; next[top] = linke[y]; linke[y] = top; point[top] = x; wis[top] = 0;
}

bool bfs()
{
	memset(tab, -1, sizeof(tab));
	int l = 1, r = 1; que[l] = s; tab[s] = 1;
	for (;l <= r;l++)
	  for (int ke = linke[que[l]]; ke; ke = next[ke])
      {
	    if (tab[point[ke]]==-1 && wis[ke])
	      tab[que[++r] = point[ke]] = tab[que[l]] + 1;
	    if (que[r] == t) return true;
	  }
	return false;
}

void improve()
{
	flow = oo;
	for (int x = prev[t]; x ;x = prev[x])
	  if (wis[shift[x]] < flow) p = x, flow = wis[shift[x]];
	for (int x = prev[t]; x ;x = prev[x])
	  wis[shift[x]] -= flow, wis[shift[x]^1] += flow;
	maxflow += flow;
}
void dfs()
{
	memcpy(shift, linke, sizeof(shift));
	for (int x = s, y;x ;)
	{
		int flag = false;
		for (int ke = shift[x]; ke; ke = next[ke])
		  if (wis[ke] && (tab[y = point[ke]] == tab[x] + 1))
		  {
				prev[y] = x; shift[x] = ke; x = y; flag = true;
				if (x == t) improve(), x = p;
				if (flag) break;
		  }
		if (!flag) tab[x] = -1,x = prev[x]; 
	}
}
int main()
{
	int i, j, k, ii, jj;
	freopen("input.txt", "r", stdin);
	freopen("output.txt", "w", stdout);
	scanf("%d%d",&n, &m); top = 1;
	mx[1]=1; mx[2]=0; mx[3]=-1; mx[4]=0;
	my[1]=0; my[2]=1; my[3]=0; my[4]=-1;
	for (i = 1; i <= n; i++)
	  for (j = 1; j <= m; j++)
	    scanf("%d", &a[i][j]), tot+= a[i][j];
	s = n*m+1; t = s+1;
	for (i = 1; i <= n; i++)
	  for (j = 1; j <= m; j++)
	    if (i+j&1) 
		   link(s, (i-1)*m+j, a[i][j]);
	    else 
		   link((i-1)*m+j, t, a[i][j]);
	for (i = 1; i <= n; i++)
	  for (j = 1; j <= m; j++)
	   if (i+j&1)
	    for (k = 1; k <= 4; k++)
	    {
			ii = i+mx[k]; jj= j+my[k];
			if (ii>=1&&ii<=n&&jj>=1&&jj<=m) link((i-1)*m+j, (ii-1)*m+jj, oo);
		}
	memcpy (adm, linke, sizeof(adm));
	while ( bfs()) 
	  dfs();
	printf("%d", tot-maxflow);
	return 0;
}

          话说考场上看出来问题是一个二分图最大权独立集的模型，但是不会算法=。=！ 下星期强力强力恶补图论！！！！！！