我的第一篇文章～

题外话：

　　我是江苏人。

　　高考刚刚考完。初步了解了和计算机有关的专业，发现人最多便是计算机科学与技术和软件工程了。

（就我目前粗浅的了解来看，似乎软件工程的人干的都是体力活==！）

 

　　我买了《The C++ Programming Language》的译文版（英文阅读也是个很头痛的事。。），准备暑假里学好C++的语法并能灵活运用STL，以应付ACM/ICPC的训练和比赛。我也买了《编译原理》《操作系统》和《人工智能》，这可能是以后一段时间内（没有变故的话）想看的书。

 

现在还是以学算法为主，比如——

 

回归主题：

 

      今天刚看了匈牙利算法，事实上二分图就一退化得不行的流网络（因此套用我前几天的dinic也可以的吧。。）。所以匈牙利算法就是一退化得不行的深搜增广路，退化到连流都不要存了，只要维护每个右点（我觉得这么说比较直观。。就是从左点出发）和哪个左点相连。我做的是USACO的stall4，因为数据太小了就用的数组邻接表图个方便。。。没事贴个源代码玩。。

 #include<fstream>int n,m;int outd[201];int g[201][201];int coved[201];int vis[201];bool dfs(int u){ for(int v,i=1;i<=outd[u];i++) { v=g[u][i]; if (vis[v]) continue; vis[v]=1; if (!coved[v] || dfs(coved[v])) {coved[v]=u;return true;} } return false;}int main(){ freopen("stall4.in","r",stdin); scanf("%d%d/n",&n,&m); for(int r,i=1;i<=n;i++) { scanf("%d",&outd[i]); for(int j=1;j<=outd[i];scanf("%d",&g[i][j++])); scanf("/n"); } int ans=0; for(int i=1;i<=m;coved[i++]=0); for(int i=1;i<=n;i++) { for(int j=1;j<=n;vis[j++]=0); if (dfs(i)) ans++; } freopen("stall4.out","w",stdout); printf("%d/n",ans); return 0;}

毕竟是退化的网络流算法，也非常短。。

其实上次写的dinic也很短啊。。就是USACO的ditch。顺便测了下oibh上dd牛很早前发的1000个点的稠密图的ditch，效果不错～我又用的数组邻接表。。存流还用的矩阵。。==！

#include<stdio.h>const int oo=32767;FILE *in,*out;int n,m;int g[1001][1001],ct[1001],d[1001];long c[1001][1001];int min(int a,int b){return (a<b)?a:b;}int level(){ int q[1001]; int i,u,v,t1=0,t2=1; q[1]=1; d[1]=0; for(i=2;i<=n;d[i++]=oo); do { u=q[++t1]; for(i=1;i<=ct[u];i++) { v=g[u][i]; if (c[u][v]>0 && d[v]==oo) { d[v]=d[u]+1; q[++t2]=v; } } } while(t1!=t2); return(d[n]);}long maxflow(int u,long ava){ if (u==n) return ava; int i,now=0,v,fl=0; for(i=1;i<=ct[u] && now<ava;i++) { v=g[u][i]; if (c[u][v]>0 && d[v]==d[u]+1) if(fl=maxflow(v,min(c[u][v],ava-now))) { now+=fl; c[u][v]-=fl; c[v][u]+=fl; } } if (!now) d[u]=oo; return now;}int main(){ in=fopen("ditch.in","r"); fscanf(in,"%d%d/n",&m,&n); for(int i=1;i<=n;i++) for(int j=1;j<=n;g[i][j]=0,c[i][j]=0,j++); { int l,r,cc; for(int i=0;i<m;i++) { fscanf(in,"%d%d%d/n",&l,&r,&cc); c[l][r]+=cc; } } fclose(in); for(int i=1;i<=n;i++) { ct[i]=0; for(int j=1;j<=n;j++) if (c[i][j] || c[j][i]) g[i][++ct[i]]=j; } long res=0,fl; while (level()!=oo) while (fl=maxflow(1,0x7fffffff))res+=fl; out=fopen("ditch.out","w"); fprintf(out,"%d/n",res); fclose(out); return 0;}

今天还看了下费用流，没写。我想不通为什么zkw的深搜+最短路在稀疏图上比连续spfa慢。。明天去写道employee好了。。话说那个挺抽象。。上次看了下还没理解过来。。暂时不深究了。。把USACO刷完先。。