1. 程式人生 > >BZOJ - 2115 獨立回路 線性基

BZOJ - 2115 獨立回路 線性基

不能 block AD 重復 int #define mat 其他 ++

題意:給定一個圖集\((V,E)\),求路徑\(1...n\)的最大異或和,其中重復經過的部分也會重復異或

所求既任意一條\(1...n\)的路徑的異或和,再異或上任意獨立回路的組合的異或和(仔細想想,異或的過程是不是不斷抵消並選取更優異或路徑的過程?)
因此dfs返向邊把環的異或值弄出來丟入線性基中貪心選取即可

順便轉載一下菊苣的獨立回路小姿勢(電路課似乎講過但我摸魚了XD)

首先有個結論:一個無向連通圖G中有且僅有M-N+1個獨立回路。
獨立回路是指任意一個都不能由其他回路構成。
引用一段數學歸納法證明:
“M=N-1時,樹,結論成立
設M=K時結論成立,當M=K+1時,任取G中一條邊e,G-e中有K-N+1個獨立回路,且
任取一個包含e的回路C,顯然獨立於之前的回路
任意兩個包含e的回路C1與C2,C12=C1+C2是G-e的回路,C2不獨立
故能且僅能增加一個包含e的獨立回路
從而G中恰有(K+1)-N+1個獨立回路,證畢”

另外,樹只要添加一條邊肯定成環,這是歸納最開始使用到的地方

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<string>
#include<vector>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<bitset> #define rep(i,j,k) for(register int i=j;i<=k;i++) #define rrep(i,j,k) for(register int i=j;i>=k;i--) #define erep(i,u) for(register int i=head[u];~i;i=nxt[i]) #define iin(a) scanf("%d",&a) #define lin(a) scanf("%lld",&a) #define din(a) scanf("%lf",&a) #define s0(a) scanf("%s",a)
#define s1(a) scanf("%s",a+1) #define print(a) printf("%lld",(ll)a) #define enter putchar(‘\n‘) #define blank putchar(‘ ‘) #define println(a) printf("%lld\n",(ll)a) #define IOS ios::sync_with_stdio(0) using namespace std; const int MAXN = 3e5+11; const double EPS = 1e-7; typedef long long ll; typedef unsigned long long ull; const ll MOD = 1e9+7; unsigned int SEED = 17; const ll INF = 1ll<<60; ll read(){ ll x=0,f=1;register char ch=getchar(); while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } int to[MAXN<<1],nxt[MAXN<<1];ll cost[MAXN<<1]; int head[MAXN<<1],tot; void init(){ memset(head,-1,sizeof head); tot=0; } void add(int u,int v,ll w){ to[tot]=v; cost[tot]=w; nxt[tot]=head[u]; head[u]=tot++; swap(u,v); to[tot]=v; cost[tot]=w; nxt[tot]=head[u]; head[u]=tot++; } bool vis[MAXN]; ll dis[MAXN],tot2; ll a[MAXN],b[66],n; void dfs(int u,int fa){ vis[u]=1; erep(i,u){ ll v=to[i],w=cost[i]; if(v==fa)continue; if(!vis[v]){ dis[v]=dis[u]^w; dfs(v,u); }else{ a[++tot2]=dis[u]^dis[v]^w;//環的異或值 } } } void cal(int n){ memset(b,0,sizeof b); rep(i,1,n){ rrep(j,62,0){ if(a[i]>>j&1){ if(b[j]) a[i]^=b[j]; else{ b[j]=a[i]; rrep(k,j-1,0) if(b[k]&&(b[j]>>k&1))b[j]^=b[k]; rep(k,j+1,62) if(b[k]>>j&1) b[k]^=b[j]; break; } } } } } int main(){ int n,m; while(cin>>n>>m){ init(); rep(i,1,m){ int u=read(); int v=read(); ll w=read(); add(u,v,w); } tot2=0; memset(vis,0,sizeof vis); memset(dis,0,sizeof dis); dfs(1,-1); cal(tot2); ll ans=dis[n]; rep(i,0,62) ans=max(ans,ans^b[i]); println(ans); } return 0; }

BZOJ - 2115 獨立回路 線性基