「一本通 4.4 例 2」暗的连锁（LOJ 10131）

看这里

原题来自：POJ 3417(点不开的)

Dark 是一张无向图，图中有 N 个节点和两类边，一类边被称为主要边，而另一类被称为附加边。Dark 有 N–1 条主要边，并且 Dark 的任意两个节点之间都存在一条只由主要边构成的路径。另外，Dark 还有 M 条附加边。

你的任务是把 Dark 斩为不连通的两部分。一开始 Dark 的附加边都处于无敌状态，你只能选择一条主要边切断。一旦你切断了一条主要边，Dark 就会进入防御模式，主要边会变为无敌的而附加边可以被切断。但是你的能力只能再切断 Dark 的一条附加边。

现在你想要知道，一共有多少种方案可以击败 Dark。注意，就算你第一步切断主要边之后就已经把 Dark 斩为两截，你也需要切断一条附加边才算击败了 Dark。

输入格式

第一行包含两个整数 N 和 M；

之后 N–1 行，每行包括两个整数 A 和 B，表示 A 和 B 之间有一条主要边；

之后 M行以同样的格式给出附加边。

输出格式

输出一个整数表示答案。

样例

样例输入

4 1 
1 2 
2 3 
1 4 
3 4

样例输出

3

 

嗯，根据题目条件，应该能看出来，这张图是有n个点的树+m条边。

熟知，一个树，每多加一条边就多出一个环。固有m个环。（附边x-y对应的环是x---......---lca(x,y)---......---y）

显然目前（破坏之前）这张图可以用主边联通，而我们只能破坏一条主边

假设我们破坏的主边是边u-fa（以下如没有特别说明，u-fa都代表fa是u的父结点   p.s.  我们默认1是根，剩下的边父子关系随便）

所以如果我们想成功，必须破坏u-fa之间的连边关系

1.u-fa不在任何一个环中，那么已经成功了，随便删一条附边即可

2.u-fa恰在一个环中。那么这个环中必恰有一个附边。设为x-y

因为边是无向的，所以此时u-fa仍是联通的。则只能删去附边x-y

3.u-fa在至少两个环中。那么这这几个环所对应的附边必然不相同，无法同时删去

所以如果将每个环所经过的边的权值加1（初始为0）

那么只需求每条边的覆盖次数

利用差分的思想，连边x-y时，w[x]++,w[y]++,w[lca(x,y)]减2。w初始值为0

（x---......---lca(x,y)时w[x]++,w[lca(x,y)]--     y---......---lca(x,y)时w[x]++,w[lca(x,y)]--） 

p.s.  有没有想到树状数组的差分？路径上每一个点的p值++（初始为0），w其实就是差分。类似的，w的和可以表示某一个值。

记边u-fa覆盖次数为f[u]，则f[u]即为u及其子树的w值之和

然后就做完了

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
using namespace std;
const int N=100010;
int n,m;
int ans;
int head[N],cnt;
struct edge{int to,next;}e[N<<1];
void add(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=head[u];head[u]=cnt;
    e[++cnt].to=u;e[cnt].next=head[v];head[v]=cnt;
}
int d[N],anc[N][20],f[N],w[N];
void dfs1(int u,int fa)
{
    for(int i=head[u];i;i=e[i].next)
    {
        int v=e[i].to;if(v==fa)continue;
        d[v]=d[u]+1;anc[v][0]=u;dfs1(v,u);
    }
}
void init()
{
    for(int j=1;j<=18;j++)
    for(int i=1;i<=n;i++)anc[i][j]=anc[anc[i][j-1]][j-1];
}
int LCA(int u,int v)
{
    if(d[u]<d[v])swap(u,v);
    for(int i=18;i>=0;i--)if(d[anc[u][i]]>=d[v])u=anc[u][i];
    if(u==v)return u;
    for(int i=18;i>=0;i--)if(anc[u][i]!=anc[v][i])u=anc[u][i],v=anc[v][i];
    return anc[u][0];
}
void dfs2(int u,int fa)
{
    for(int i=head[u];i;i=e[i].next)
    {
        int v=e[i].to;if(v==fa)continue;
        dfs2(v,u);f[u]+=f[v];//注意到这里是先dfs  把我坑了一下
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int u,v;scanf("%d%d",&u,&v);
        add(u,v);
    }
    d[1]=1;dfs1(1,0);init();
    for(int i=1;i<=m;i++)
    {
    	int u,v;scanf("%d%d",&u,&v);
    	w[u]++;w[v]++;w[LCA(u,v)]-=2;
    }
    for(int i=1;i<=n;i++)f[i]=w[i];
    dfs2(1,0);
    for(int i=2;i<=n;i++)
	{
		if(f[i]==0)ans+=m;
		else if(f[i]==1)ans++;
	}
	//for(int i=1;i<=n;i++)printf("%d\n",f[i]);
	printf("%d",ans);
    return 0;
}
/*
9 2
1 2
1 3
1 4
2 5
2 6
4 7
4 8
7 9
6 7
8 9


9
*/