日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

spoj839 Optimal Marks(最小割,dinic)

發(fā)布時(shí)間:2024/7/5 编程问答 79 豆豆
生活随笔 收集整理的這篇文章主要介紹了 spoj839 Optimal Marks(最小割,dinic) 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

題目大意:

給你一個(gè)無向圖\(G(V,E)\)。 每個(gè)頂點(diǎn)都有一個(gè)int范圍內(nèi)的整數(shù)的標(biāo)記。 不同的頂點(diǎn)可能有相同的標(biāo)記。
對(duì)于邊\((u,v)\),我們定義\(Cost(u,v)=mark [u]\ \ xor\ \ mark [v]\)
現(xiàn)在我們知道某些節(jié)點(diǎn)的標(biāo)記了。你需要確定其他節(jié)點(diǎn)的標(biāo)記,以使邊的總成本盡可能小。
最后要求輸出的每個(gè)點(diǎn)的標(biāo)號(hào)

QwQ一看到這種跟位運(yùn)算有關(guān)題目,就會(huì)想到按位來處理

仔細(xì)考慮,發(fā)現(xiàn)這個(gè)題滿足最小割的模型,對(duì)于每一位,當(dāng)時(shí)將所有點(diǎn)的對(duì)應(yīng)位分成0,或者是1

那么,我們按位來,假設(shè)當(dāng)前位是\(i\),對(duì)于已經(jīng)知道編號(hào)的點(diǎn)\(x\),如果當(dāng)前位是1的話,我們\(insert(s,x,inf)\),否則\(insert(x,t,inf)\)表示,這個(gè)點(diǎn)是0還是1,同時(shí)inf的原因是給定的點(diǎn)的編號(hào)的不能改的
同時(shí)對(duì)于原圖的邊\(u->v\),我們只需要\(insert(u,v,1),insert(v,u,1)\) 表示這兩個(gè)點(diǎn)的當(dāng)前位是否相同,最后跑\(dinic\),剩下的殘余網(wǎng)絡(luò)中,與s相連,且沿途流量\(>0\)的,就是1,否則就是0

大致就是這樣,最后千萬別忘記:
1.編號(hào)可能是0
2.初始化數(shù)組

#include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> #include<queue>using namespace std;inline int read() {int x=0,f=1;char ch=getchar();while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f; }const int maxn = 1010; const int maxm = 200010; const int inf = 1e9;int point[maxn],nxt[maxm],to[maxm],val[maxm]; int h[maxn]; int num[maxn]; int ans[maxn]; int x[maxm],y[maxm]; int n,m,cnt=1; int s,t; int vis[maxn]; queue<int> q;void addedge(int x,int y,int w){nxt[++cnt]=point[x];to[cnt]=y;val[cnt]=w;point[x]=cnt; }void init() {cnt=1;memset(point,0,sizeof(point));memset(vis,0,sizeof(vis)); }void insert(int x,int y,int w) {addedge(x,y,w);addedge(y,x,0); }bool bfs(int s) {memset(h,-1,sizeof(h));h[s]=0;q.push(s);while (!q.empty()){int x = q.front();q.pop();for (int i=point[x];i;i=nxt[i]){int p = to[i];if (val[i]>0 && h[p]==-1){h[p]=h[x]+1;q.push(p);}}}if (h[t]==-1) return false;else return true; }int dfs(int x,int low) {if (x==t || low==0) return low;int totflow=0;for (int i=point[x];i;i=nxt[i]){int p = to[i];if (val[i]>0 && h[p]==h[x]+1){int tmp = dfs(p,min(val[i],low));val[i]-=tmp;val[i^1]+=tmp;low-=tmp;totflow+=tmp;if (low==0) return totflow;}}if (low>0) h[x]=-1;return totflow; }int dinic(){int ans=0;while (bfs(s)){ans+=dfs(s,inf);} }void dfs1(int x,int d) {vis[x]=1;ans[x]|=(1 << d); for (int i=point[x];i;i=nxt[i]){int p = to[i];if (!vis[p] &&val[i]>0){dfs1(p,d);}} } void build(int xx) {init();s=n+10;t=s+1; for (int i=1;i<=n;i++){if (num[i]!=-1){if (num[i] & (1<<xx)) insert(s,i,inf);else insert(i,t,inf);}}for (int i=1;i<=m;i++){insert(x[i],y[i],1);insert(y[i],x[i],1);}dinic();dfs1(s,xx); }int T; int main() {scanf("%d",&T);while (T--){memset(num,-1,sizeof(num));memset(ans,0,sizeof(ans));init();n=read(),m=read();for (int i=1;i<=m;i++) x[i]=read(),y[i]=read();int k;k=read();for (int i=1;i<=k;i++){int oo;oo=read();num[oo]=read();}for (int i=0;i<=32;i++){build(i);}for (int i=1;i<=n;i++) printf("%d\n",ans[i]);}return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/yimmortal/p/10160835.html

總結(jié)

以上是生活随笔為你收集整理的spoj839 Optimal Marks(最小割,dinic)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。