【JZOJ4784】【NOIP2016提高A组模拟9.15】Map
题目描述
输入
输出
样例输入
4 4 2
1 2
2 3
3 2
3 4
1 2
1 4
样例输出
14
数据范围
样例解释
upd:保证原图连通。
“不相交路径”的定义为不存在相同的边。可以存在相同的点。重边视为不同的边。
对于样例:
原图有2个安全点对为(2,3),(3,2)
询问1答案为4,新增的安全点对为(1,2),(1,3),(2,1)(3,1)
询问2答案为10,新增的安全点对为(1,2),(1,3),(1,4),(2,1)(2,4),(3,1),(3,4),(4,1),(4,2),(4,3)
因此输出14
解法
利用边双连通分量缩点,再利用lca统计答案。
代码
#include<iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
#define ll long long
#define ln(x,y) (ll)(log(x)/log(y))
#define sqr(x) ((x)*(x))
using namespace std;
const char* fin="aP3.in";
const char* fout="aP3.out";
const ll inf=0x7fffffff;
const ll maxn=400007,maxm=maxn*4,maxk=22;
ll n,m,t,i,j,k,l,tot,limit,tmp,cnt;
ll fi[maxn],ne[maxm],la[maxm],de[maxn];
ll fa[maxn][maxk],sum[maxn][maxk],qsum[maxn][maxk];
ll dfn[maxn],low[maxn],num,stack[maxn],blg[maxn];
ll az[maxn];
ll ans;
bool bz[maxn],cz[maxn],dz[maxn];
ll read(){
ll x=0;
char ch=getchar();
while (ch<'0' || ch>'9') ch=getchar();
while (ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
return x;
}
void add_line(ll a,ll b){
tot++;
ne[tot]=fi[a];
la[tot]=b;
fi[a]=tot;
}
void dfs(ll v,ll from){
ll i,j,k;
cz[v]=true;
fa[v][0]=from;
qsum[v][0]=sqr(sum[v][0]);
de[v]=de[from]+1;
for (i=1,j=ln(de[v],2);i<=j;i++){
k=fa[v][i-1];
fa[v][i]=fa[k][i-1];
sum[v][i]=sum[v][i-1]+sum[k][i-1];
qsum[v][i]=qsum[v][i-1]+qsum[k][i-1];
}
for (k=fi[v];k;k=ne[k]){
if (!cz[az[blg[la[k]]]]){
dfs(az[blg[la[k]]],v);
}
}
}
void up(ll &a,ll i,ll &j,ll &k){
j+=sum[a][i];
k+=qsum[a][i];
a=fa[a][i];
}
ll lca(ll a,ll b){
ll i,j=0,k=0;
if (de[a]<de[b]) swap(a,b);
for (i=ln(de[a]-de[b],2);i>=0;i--) if (de[fa[a][i]]>de[b]) up(a,i,j,k);
if (de[a]!=de[b]) up(a,0,j,k);
for (i=ln(de[a],2);i>=0;i--) if (fa[a][i]!=fa[b][i]) up(a,i,j,k),up(b,i,j,k);
if (a!=b) up(a,0,j,k),up(b,0,j,k);
j+=sum[a][0];
k+=qsum[a][0];
return sqr(j)-k;
}
void tarjan(ll v,ll from){
ll i,j,k;
dfn[v]=low[v]=++num;
bz[j=stack[++stack[0]]=v]=true;
for (k=fi[v];k;k=ne[k])
if (!dz[(k+1)/2]){
dz[(k+1)/2]=true;
if (!dfn[la[k]]){
tarjan(la[k],v);
low[v]=min(low[la[k]],low[v]);
}else if (bz[la[k]]) low[v]=min(low[v],dfn[la[k]]);
}
if (low[v]==dfn[v]){
while (stack[stack[0]]!=j){
blg[stack[stack[0]]]=v;
bz[stack[stack[0]--]]=false;
}
blg[stack[stack[0]]]=v;
bz[stack[stack[0]--]]=false;
}
}
int main(){
cnt=n=read();m=read();t=read();
for (i=1;i<=m;i++){
j=read();k=read();
add_line(j,k);
add_line(k,j);
}
tarjan(1,0);
for (i=1;i<=n;i++){
if (!az[blg[i]]) az[blg[i]]=++cnt;
sum[az[blg[i]]][0]++;
for (k=fi[i];k;k=ne[k]) add_line(az[blg[i]],la[k]);
}
dfs(n+1,0);
for (i=1;i<=t;i++){
j=read();
k=read();
ans+=lca(az[blg[j]],az[blg[k]]);
}
printf("%lld\n",ans);
return 0;
}
启发
边双连通分量在某种程度上等同于强连通分量。
【JZOJ4784】【NOIP2016提高A组模拟9.15】Map的更多相关文章
- NOIP2016提高A组模拟10.15总结
第一题,就是将原有的式子一步步简化,不过有点麻烦,搞了很久. 第二题,枚举上下边界,维护一个单调队列,二分. 比赛上没有想到,只打了个暴力,坑了80分. 第三题,贪心,最后的十多分钟才想到,没有打出来 ...
- 【NOIP2016提高A组模拟10.15】打膈膜
题目 分析 贪心, 先将怪物按生命值从小到大排序(显然按这个顺序打是最优的) 枚举可以发对少次群体攻击, 首先将所有的群体攻击发出去, 然后一个一个怪物打,当当前怪物生命值大于2,如果还有魔法值就放重 ...
- 【NOIP2016提高A组模拟10.15】最大化
题目 分析 枚举两个纵坐标i.j,接着表示枚举区域的上下边界, 设对于每个横坐标区域的前缀和和为\(s_l\),枚举k, 显然当\(s_k>s_l\)时,以(i,k)为左上角,(j,k)为右下角 ...
- 【NOIP2016提高A组模拟10.15】算循环
题目 分析 一步步删掉循环, 首先,原式是\[\sum_{i=1}^n\sum_{j=1}^m\sum_{k=i}^n\sum_{l=j}^m\sum_{p=i}^k\sum_{q=j}^l1\] 删 ...
- 【NOIP2016提高A组模拟9.15】Map
题目 分析 发现,当原图是一棵树的时候,那么新建一条边后,就会变成环套树, 而环内的所有点对都是安全点对,如果环中有k个点,答案就是\(k(k-1)\) 联想到,当把原图做一遍tarjan缩点,每个环 ...
- 【NOIP2016提高A组模拟9.15】Osu
题目 分析 考虑二分答案, 二分小数显然是不可取的,那么我们将所有可能的答案求出来,记录在一个数组上,排个序(C++调用函数很容易超时,手打快排,时间复杂度约为\(O(>8*10^7)\),但相 ...
- 【NOIP2016提高A组模拟9.15】Math
题目 分析 因为\((-1)^2=1\), 所以我们只用看\(\sum_{j=1}^md(i·j)\)的值模2的值就可以了. 易证,一个数x,只有当x是完全平方数时,d(x)才为奇数,否则为偶数. 那 ...
- 【NOIP2016提高A组模拟8.15】Garden
题目 分析 其实原题就是[cqoi2012][bzoj2669]局部极小值. 有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次.如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点) ...
- 【NOIP2016提高A组模拟8.15】Throw
题目 分析 首先对于一个状态(a,b,c),假定a<=b<=c: 现在考虑一下这个状态,的转移方案: \[1,中间向两边跳(a,b,c)-->(a*2-b,a,c).(a,b,c)- ...
随机推荐
- javascript 的原始数据类型
javascript 的原始数据类型有五种,分别是Number.String.Boolean.null.undefined 原始类型数据的特点:值保存在变量本地,且赋值给其他变量后,其它变量改变 ...
- JS钩子的机制与实现
[什么是钩子] 接触过WordPress的朋友都知道,WP的程序中可以执行类似钩子的函数,当然是这PHP实现的钩子.在JavaScript中一样可以实现类似的功能. 用一句话来形容一下:钩子是将需要执 ...
- bash: express: command not found及vue连接数据库调接口
今天在使用express -e . 的命令时,cmd给我报了一段不识别的错误: bash: express: command not found ,在网上查了一下,有人指出是express4的版本将命 ...
- 简单的选项卡制作(原生JS)
<!doctype html> <html> <head> <meta charset="utf-8"> <title> ...
- NGINX模块开发 之 验证URL參数
版权声明:本文为博主原创文章,未经博主同意不得转载. https://blog.csdn.net/RoyalApex/article/details/26404379 作者:邹祁峰 邮箱:Qifeng ...
- WPF 禁用中文
<TextBox InputMethod.IsInputMethodEnabled="False" />
- <init>与<clinit>的区别
在编译生成class文件时,会自动产生两个方法,一个是类的初始化方法<clinit>, 另一个是实例的初始化方法<init> <clinit>:在jvm第一次加载c ...
- Jquery手机下拉刷新,下拉加载数据
一.Jquery手机下拉刷新,下拉加载数据.附加有源码 <!DOCTYPE html> <html> <head> <title>手机</titl ...
- Mac下,spacy配置
pip3 install -U spacy -i http://pypi.douban.com/simple --trusted-host pypi.douban.com python3 -m spa ...
- springcloud的服务提供者与服务消费者
1.说明 springcloud中由服务消费者调用服务提供者一共有两种方法rest和feign 2.feign (1)使用feign的方式进行服务调,搭建服务提供者. 创建一个web项目(服务提供者) ...