Loj#3130-「COCI 2018.12」Praktični【线性基】
生活随笔
收集整理的這篇文章主要介紹了
Loj#3130-「COCI 2018.12」Praktični【线性基】
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
正題
題目鏈接:https://loj.ac/p/3130
題目大意
給出nnn個點(diǎn)mmm條邊的一張有權(quán)無向圖,你每次可以選擇一個邊集異或上一個值,要求最少次數(shù)使得所有簡單環(huán)異或和都為000。
1≤n,m≤1051\leq n,m\leq 10^51≤n,m≤105
解題思路
先找一棵生成樹,然后每條非樹邊都會產(chǎn)生一個簡單環(huán),顯然這些簡單環(huán)合法了其他的也一定合法。
而肯定存在一種最優(yōu)的方案是只改非樹邊,因為如果該樹邊首先我們可以一次改一個集合所有必須改一條樹邊會對多個簡單環(huán)產(chǎn)生不同的影響,而如果我們異或的是valvalval,那么我們能做到的只是讓某個簡單環(huán)異或上valvalval(奇數(shù)條樹邊操作),或者不異或上valvalval(偶數(shù)條樹邊操作),所以是和我們操作非樹邊能做到的是相同的。
那么問題就變?yōu)橐阎恍?shù)要異或多少,求最少操作次數(shù)了,我們用線性基求出最小的線性空間就好了。
時間復(fù)雜度:O(n+mlog?m)O(n+m\log m)O(n+mlogm)
code
#include<cstdio> #include<cstring> #include<algorithm> #include<vector> #define mp(x,y) make_pair(x,y) using namespace std; const int N=1e5+10; struct node{int to,next,w; }a[N<<1]; int n,m,tot,k,w[N],d[N],ls[N],v[N]; bool vis[N]; void addl(int x,int y,int w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return; } void dfs(int x,int fa){vis[x]=1;for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(y==fa)continue;if(vis[y]){if(i&1)v[(i+1)/2]=w[x]^w[y]^a[i].w;}else{w[y]=w[x]^a[i].w;dfs(y,x);}}return; } void ins(int x){for(int i=29;i>=0;i--)if((x>>i)&1){if(d[i])x^=d[i];else {d[i]=x;k+=(x!=0);break;}}return; } int main() {scanf("%d%d",&n,&m);for(int i=1,x,y,w;i<=m;i++){scanf("%d%d%d",&x,&y,&w);addl(x,y,w);addl(y,x,w);}dfs(1,0);for(int i=1;i<=m;i++)ins(v[i]);printf("%d\n",k);for(int i=0;i<=29;i++)if(d[i]){int cnt=0;for(int j=1;j<=m;j++)cnt+=((v[j]>>i)&1);printf("%d %d ",d[i],cnt);for(int j=1;j<=m;j++)if((v[j]>>i)&1)printf("%d ",j);putchar('\n');}return 0; }總結(jié)
以上是生活随笔為你收集整理的Loj#3130-「COCI 2018.12」Praktični【线性基】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 菜鸟与速卖通升级中东、北美、巴西等国际快
- 下一篇: P5299-[PKUWC2018]Sla