原始對偶費用流ver2 0 hdu4744

這場比賽我就不多說什麼了。。。幾乎整場都在寫計算幾何，恰好三維計算幾何又是我的軟肋，後面的題基本都沒看。。。

建圖就不多說了，挺水的，標算不是費用流，但我聽說有人zkw費用流過了，於是我就寫個原始對偶費用流，但是居然超了，本機測試60組極限1050+ms==。。。於是進行優化，lyp跟我說過一次反向spfa就可以dij下去了，但是自己推了一下,有乙個對原圖的邊修改的過程.

考慮i->j,w[i][j]>0,操作應該是w[i][j]-=d[i]-d[j],因為d[i]-d[j]<=w[i][j]，所以w[i][j]-(d[i]-d[j])>=0,但是考慮反邊w[j][i]-=d[j]-d[i],也就是-w[i][j]+d[i]-d[j],因為d[i]-d[j]<=w[i][j]，所以這個反邊貌似還是負的，因此修改過一次後不是所想的全為正權邊，從而不能dij，lyp以前看到的那篇**找不到了，所以也沒法繼續**下去，我自己也不是很理解這個演算法的理論基礎- -，於是另闢蹊徑，用了春瀟最喜歡的優先佇列+spfa，瞬間600+ms過了...

原始對偶費用流對於點度比較大的圖貌似確實有些慢...

#include #include #include #include #include #include #include #define sqr(x) ((x)*(x))

const int oo=1073741819;

using namespace std;

int tail[250],d[250],v[250],flag[250],x[250],y[250],z[250],w[250],b[1050],p[250];

int next[150000],sora[150000],flow[150000],cost[150000],st[150000],po[150000];

int n,ss,s,t,phi,ans,sum,w_time,e,m1;

void change(int x,int w)

bool spfa(int s,int t)

return 1;

}int dfs(int x,int low)

int sum=0,tmp,ne;

flag[x]=w_time;

for (int i=x;next[i];)

bool spfa(int s,int t)

{ int ne,na;

for (int i=0;i<=t;i++) d[i]=p[i]=oo;

p[t]=0;

change(t,p[t]);

for (;d[b[1]]!=oo;) {

ne=b[1];

change(ne,oo);

for (int i=ne;next[i];) {

i=next[i],na=sora[i];

if (flow[po[i]] && p[ne]+cost[po[i]]+p[ne]=oo) return 0;

for (int i=1;i<=t;i++)

if (p[i]

原始對偶費用流ver2 0 hdu4744

原始對偶費用流

學習筆記最小費用流之原始對偶

Primal Dual 原始對偶

原始對偶費用流ver2 0 hdu4744

原始對偶費用流

學習筆記 最小費用流之原始對偶

Primal Dual 原始對偶

相關推薦

學習筆記最小費用流之原始對偶