POJ - 3417 Network LCA+树上差分
生活随笔
收集整理的這篇文章主要介紹了
POJ - 3417 Network LCA+树上差分
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
思路:首先來說,給一個樹加一條邊肯定要構成一個環,我們假設加了該邊后產生的環上的每一條邊都累計加一。
假設這條邊是a到b,那么其實就是原樹a->lca(a,b)到b->lca(a,b)上的邊累計加一。
這個有什么用呢?
第一:當原樹上某一條邊的累加值為0時,也就是這一條邊與所有的假邊都沒有構成環,那么你只要把這一條邊刪除了,再任意刪除一條假邊,肯定把圖分成兩份。答案+M
第二:當原樹上某一條邊的累加值為1時,也就是說原樹上的這條邊,有且僅有一條假邊與他構成了環,那么此時將這條邊與他對應的假邊刪除即可,原樹上的這條邊也只有這一種情況可以分成兩份。答案+1
第三:當原樹上某一條邊的累加值大于等于2時,因為你始終只能刪除一條假邊,所以刪除該原邊對應的情況無論如何都無法將他分成兩部分。答案+0
所以將M條邊跑一邊找到lca并且進行累加即可求出答案,但是問題就是怎么來進行這個累加,這里就用到了類似前綴和的方式或者dp的方式來求:用dif數組記錄每一次累加的段,然后再用過一次dfs來將狀態進行轉換。
dif[x]++,dif[y]++,dif[lca(x,y)]-=2;
sta數組表示累加后的結果:
int dfs(int u)
{sta[u]=dif[u];for(int i=head[u];~i;i=edge[i].p){int v=edge[i].to;if(v==f[u][0]) continue;sta[u]+=dfs(v);//遞歸下去}return sta[u];
}
AC代碼:
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define me(a,b) memset(a,b,sizeof(a))
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
const int dir[4][2]= {0,-1,-1,0,0,1,1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=1e5+5;
int f[maxn][33],d[maxn],dif[maxn],head[maxn],sta[maxn];
int sign,n,m;
struct node
{int to,p;
}edge[maxn<<1];
void add(int u,int v)
{edge[sign]=node{v,head[u]};head[u]=sign++;
}
void init()
{memset(dif,0,sizeof(dif));sign=0;for(int i=0;i<=n;i++)head[i]=-1;
}
void lca_dfs(int u)
{for(int i=1;(1<<i)<=n;i++)f[u][i]=f[f[u][i-1]][i-1];for(int i=head[u];~i;i=edge[i].p){int v=edge[i].to;if(v==f[u][0]) continue;f[v][0]=u;d[v]=d[u]+1;lca_dfs(v);}
}
int dfs(int u)
{sta[u]=dif[u];for(int i=head[u];~i;i=edge[i].p){int v=edge[i].to;if(v==f[u][0]) continue;sta[u]+=dfs(v);}return sta[u];
}
int lca(int a,int b)
{if(d[a]<d[b]) swap(a,b);int x=d[a]-d[b];for(int i=0;(1<<i)<=x;i++)if((1<<i)&x) a=f[a][i];if(a!=b){for(int i=(int)log2(n);i>=0;i--)if(f[a][i]!=f[b][i]) a=f[a][i],b=f[b][i];a=f[a][0];}return a;
}
int main()
{int x,y;while(cin>>n>>m){init();for(int i=1;i<n;i++) scanf("%d%d",&x,&y),add(x,y),add(y,x);d[1]=0,f[1][0]=1;lca_dfs(1);for(int i=1;i<=m;i++){scanf("%d %d",&x,&y);dif[x]++,dif[y]++,dif[lca(x,y)]-=2;}dfs(1);int ans=0;for(int i=1;i<=n;i++){if(sta[i]==0&&i!=1) ans+=m;if(sta[i]==1) ans+=1;}printf("%d\n",ans);}return 0;
}
?
總結
以上是生活随笔為你收集整理的POJ - 3417 Network LCA+树上差分的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: POJ - 2763 Housewife
- 下一篇: A and B and Lecture