bzoj3331 压力(圆方树)

题目链接

圆方树

圆方树就是对于联通无向图中的每一个点双新建一个方点,与点双中的每个点连一条边,然后将原来的边删去。将原来的点看作圆点,新建的点看作方点。所以叫做圆方树。

性质

1.圆方树肯定是棵树(废话)。证明显然。
2.圆方树中与圆点相连的点肯定是方点。与方点相连的点肯定是圆点。

算法

根据圆方树的定义就可以知道。构建圆方树的过程实际上就是找点双的过程。本质上就是找割点。所以用tarjan来做就好了。将找出的点双中的点与新建的点连边即可。

思路

这个题就是圆方树的经典应用,先对于原图构建出圆方树。那么对于每个从S到T的路径。在圆方树上两点之间的简单路径上的圆点就是在原图中必须经过的点。然后再圆方树上差分一下就行了。

代码

/*
* @Author: wxyww
* @Date:   2019-01-22 20:03:52
* @Last Modified time: 2019-01-22 20:42:20
*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<bitset>
#include<vector>
using namespace std;
typedef long long ll;
const int N = 200000 + 100,logN = 20;
vector<int>E[N * 2];
ll read() {
   ll x=0,f=1;char c=getchar();
   while(c<'0'||c>'9') {
      if(c=='-') f=-1;
      c=getchar();
   }
   while(c>='0'&&c<='9') {
      x=x*10+c-'0';
      c=getchar();
   }
   return x*f;
}
struct node {
   int u,v,nxt;
}e[N * 4];
int head[N],ejs;
void add(int u,int v) {
   e[++ejs].u = u;e[ejs].v = v;e[ejs].nxt = head[u];head[u] = ejs;
}
void ADD(int u,int v) {
   // printf("!!!%d %d\n",u,v);
   E[u].push_back(v);E[v].push_back(u);
}
int top,sta[N];
int coljs,n,m,Q,cnt,dfn[N],low[N];
void tarjan(int u) {
   dfn[u] = low[u] = ++cnt;
   sta[++top] = u;
   for(int i = head[u];i;i = e[i].nxt) {
      int v = e[i].v;
      // if(v == father) continue;
      if(!dfn[v]) {
         tarjan(v);
         low[u] = min(low[u],low[v]);
         if(low[v] >= dfn[u]) {
            ++n;
            ADD(n,u);
            while(sta[top + 1] != v) {
               ADD(n,sta[top--]);
            }
         }
      }
      else low[u] = min(low[u],dfn[v]);
   }
}
int dep[N];
int lca[N][logN];
void dfs(int u,int father) {
   dep[u] = dep[father] + 1;
   for(int i = 1;i < logN;++i)
      lca[u][i] = lca[lca[u][i - 1]][i - 1];
   int k = E[u].size();
   for(int i = 0;i < k;++i) {
      int v = E[u][i];
      if(v == father) continue;
      lca[v][0] = u;
      dfs(v,u);
   }
}
int LCA(int x,int y) {
   if(dep[x] < dep[y]) swap(x,y);
   for(int i = logN - 1;i >= 0;--i) {
      if(dep[lca[x][i]] >= dep[y]) x = lca[x][i];
   }
   for(int i = logN - 1;i >= 0;--i) {
      if(lca[x][i] != lca[y][i]) {
         x = lca[x][i];y = lca[y][i];
      }
   }
   if(x != y) x = lca[x][0];
   return x;
}
int sum[N];
void dfs2(int u,int father) {
   int k = E[u].size();
   for(int i = 0;i < k;++i) {
      int v = E[u][i];
      if(v == father) continue;
      dfs2(v,u);
      sum[u] += sum[v];
   }
}
int main() {
   n = read(),m = read(),Q = read();
   int nn = n;
   for(int i = 1;i <= m;++i) {
      int u = read(),v = read();
      add(u,v);add(v,u);
   }
   tarjan(1);
   dfs(1,0);
   for(int i = 1;i <= Q;++i) {
      int u = read(),v = read();
      int L = LCA(u,v);
      // printf("%d\n",L);
      sum[u]++;sum[v]++;sum[L]--;sum[lca[L][0]]--;
   }
   dfs2(1,0);
   for(int i = 1;i <= nn;++i) printf("%d\n",sum[i]);
   return 0;
}
全部评论

相关推荐

点赞 评论 收藏
分享
程序员饺子:正常 我沟通了200多个 15个要简历 面试2个 全投的成都的小厂。很多看我是27直接不会了😅
点赞 评论 收藏
分享
05-11 11:48
河南大学 Java
程序员牛肉:我是26届的双非。目前有两段实习经历,大三上去的美团,现在来字节了,做的是国际电商的营销业务。希望我的经历对你有用。 1.好好做你的CSDN,最好是直接转微信公众号。因为这本质上是一个很好的展示自己技术热情的证据。我当时也是烂大街项目(网盘+鱼皮的一个项目)+零实习去面试美团,但是当时我的CSDN阅读量超百万,微信公众号阅读量40万。面试的时候面试官就告诉我说觉得我对技术挺有激情的。可以看看我主页的美团面试面经。 因此花点时间好好做这个知识分享,最好是单拉出来搞一个板块。各大公司都极其看中知识落地的能力。 可以看看我的简历对于博客的描述。这个帖子里面有:https://www.nowcoder.com/discuss/745348200596324352?sourceSSR=users 2.实习经历有一些东西删除了,目前看来你的产出其实很少。有些内容其实很扯淡,最好不要保留。有一些点你可能觉得很牛逼,但是面试官眼里是减分的。 你还能负责数据库表的设计?这个公司得垃圾成啥样子,才能让一个实习生介入数据库表的设计,不要写这种东西。 一个公司的财务审批系统应该是很稳定的吧?为什么你去了才有RBAC权限设计?那这个公司之前是怎么处理权限分离的?这些东西看着都有点扯淡了。 还有就是使用Redis实现轻量级的消息队列?那为什么这一块不使用专业的MQ呢?为什么要使用redis,这些一定要清楚, 就目前看来,其实你的这个实习技术还不错。不要太焦虑。就是有一些内容有点虚了。可以考虑从PR中再投一点产出
点赞 评论 收藏
分享
今天投了小鹏,收到了AI面,大概会问哪些啊?
期末一定及格:总共4个部分,心理测评、行测、然后就是问岗位、对岗位的理解、过往遇到了哪些难点怎么解决,很简单,没有什么特别专业的问题,都是一些综合素质相关的
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务