简化一下题意,我们先看成一副强连通的图,这时候应该是最简单了,去点任意点都是其他的乘积。那再加强一点难度,改为两个强连通图连接的非强连通图呢?那应该就是找出关键的那个点,并求出两边的乘积。但是一个一个去找是不可能的。

假设如图中的非绿色线是题目给的图。然后我们根据强连通分量去新建一副如图中绿色线条的图,那么这时候我们就把原图转化为以可树了。。对于每一个点我们求的是该点以及以下的乘积。然后我们从A出发这时候我们发现A点的值刚好就是整幅图的乘积。这时候如果我们需要求删除3这个点的得到的结果应该就是整一副图去除以3点及一下的乘积得到1,2的乘积,再加上3点的子树的乘积和也就是4、5 和 6、7的乘积和。

这道题目的难点就是转化建图的那一个步骤,应该说是最核心的部分。

#include<bits/stdc++.h>
#define LL long long
using namespace std; const int maxn = 1e5 + ;
const LL mod = 1e9 + ; LL n, w[maxn], vis[maxn << 1], sum[maxn << 1], pro[maxn << 1]; struct oldEdge{
int v, next;
bool flag;
};
int oldHead[maxn], ocnt;
oldEdge oedge[maxn << ]; void addOldEdge(int u, int v){
oedge[ocnt].v = v;
oedge[ocnt].flag = false;
oedge[ocnt].next = oldHead[u];
oldHead[u] = ocnt ++;
} struct newEdge{
int v, next;
};
int newHead[maxn << ], ncnt;
newEdge nedge[maxn << ]; void addNewEdge(int u, int v){
nedge[ncnt].v = v;
nedge[ncnt].next = newHead[u];
newHead[u] = ncnt ++;
} int dfn[maxn], low[maxn], root[maxn], rn, cnt;
stack<int>sta; void tarjan(int u, int fa){
dfn[u] = low[u] = ++cnt;
for(int i = oldHead[u]; i != -; i = oedge[i].next){
if(oedge[i].flag) continue;
oedge[i].flag = oedge[i ^ ].flag = true;
sta.push(i);
int v = oedge[i].v; if(dfn[v]){
low[u] = min(low[u], dfn[v]);
continue;
}
tarjan(v, fa);
low[u] = min(low[u], low[v]); if(low[v] >= dfn[u]){
rn ++;
int ek;
do{
ek = sta.top();sta.pop();
root[oedge[ek].v] = root[oedge[ek ^ ].v] = fa;
addNewEdge(rn, oedge[ek].v); addNewEdge(oedge[ek].v, rn);
addNewEdge(rn, oedge[ek ^ ].v); addNewEdge(oedge[ek ^ ].v, rn);
}while(oedge[ek ^ ].v != u);
}
}
} void dfs(int u){
vis[u] = true;
sum[u] = ;
pro[u] = (u <= n) ? w[u] : ;
for(int i = newHead[u]; i != -; i = nedge[i].next){
int v = nedge[i].v;
if(vis[v]) continue;
dfs(v);
if(u <= n)
sum[u] = (sum[u] + pro[v]) % mod;
pro[u] = pro[u] * pro[v] % mod;
}
} LL inv(LL a){
int p = mod - ;
LL ret = ;
while(p){
if(p & )ret = ret * a % mod;
a = a * a % mod;
p >>= ;
}
return ret;
} void init(){
memset(root, , sizeof(root));
memset(newHead, -, sizeof(newHead));
memset(oldHead, -, sizeof(oldHead));
memset(dfn, , sizeof(dfn));
memset(vis, false, sizeof(vis));
ncnt = ocnt = cnt = ;
} int main(){
int T, m, a, b;scanf("%d",&T);
while(T --){
scanf("%lld%d",&n,&m);
init();rn = n;
for(int i = ; i <= n; i ++)scanf("%lld",&w[i]);
for(int i = ; i < m; i ++){
scanf("%d%d",&a,&b);
addOldEdge(a,b);
addOldEdge(b,a);
}
for(int i = ; i <= n; i ++)
if(!dfn[i])tarjan(i, rn + ); LL tot = ;
for(int i = ; i <= n; i ++){
if(vis[i]) continue;
if(root[i]){
dfs(root[i]);
tot = (tot + pro[root[i]]) %mod;
}else
tot = (tot + w[i]) % mod;
}
LL ans = ;
for(int i = ; i <= n; i ++){
if(root[i]){
LL temp = ((tot - pro[root[i]] + pro[root[i]] * inv(pro[i]) + sum[i]) % mod + mod) * i % mod;
ans = (ans + temp) % mod;
}
else
ans = ((ans + (tot - w[i]) * i) % mod + mod) % mod;
}
printf("%lld\n",ans);
}
return ;
}

Fantasia (点强连通分量建图 + 树形DP)的更多相关文章

  1. bzoj 4871: [Shoi2017]摧毁“树状图” [树形DP]

    4871: [Shoi2017]摧毁"树状图" 题意:一颗无向树,选两条边不重复的路径,删去选择的点和路径剩下一些cc,求最多cc数. update 5.1 : 刚刚发现bzoj上 ...

  2. 【HDU5934】Bomb——有向图强连通分量+重建图

    题目大意 二维平面上有 n 个爆炸桶,i−thi-thi−th爆炸桶位置为 (xi,yi)(x_i, y_i)(xi​,yi​) 爆炸范围为 rir_iri​ ,且需要 cic_ici​ 的价格引爆, ...

  3. 强连通 反向建图 hdu3639

    Hawk-and-Chicken Time Limit: 6000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...

  4. poj1949Chores(建图或者dp)

    /* 题意:n个任务,有某些任务要在一些任务之前完成才能开始做! 第k个任务的约束只能是1...k-1个任务!问最终需要最少的时间完成全部的 任务! 思路:第i个任务要在第j个任务之前做,就在i,j之 ...

  5. 深探树形dp

    看到同学在写一道树形dp,好奇直接拿来写,发现很不简单. 如图,看上去是不是很像选课,没错这不是选课,升级版吧,多加了点东西罢了.简单却调了一晚上和一上午. 思路:很简单强联通分量+缩点+树形dp.直 ...

  6. Tarjan算法 求 有向图的强连通分量

    百度百科 https://baike.baidu.com/item/tarjan%E7%AE%97%E6%B3%95/10687825?fr=aladdin 参考博文 http://blog.csdn ...

  7. poj 2186 tarjan求强连通分量

    蕾姐讲过的例题..玩了两天后才想起来做 貌似省赛之后确实变得好懒了...再努力两天就可以去北京玩了! 顺便借这个题记录一下求强连通分量的算法 1 只需要一次dfs 依靠stack来实现的tarjan算 ...

  8. POJ 3107.Godfather 树形dp

    Godfather Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 7536   Accepted: 2659 Descrip ...

  9. Tarjan算法初探 (1):Tarjan如何求有向图的强连通分量

    在此大概讲一下初学Tarjan算法的领悟( QwQ) Tarjan算法 是图论的非常经典的算法 可以用来寻找有向图中的强连通分量 与此同时也可以通过寻找图中的强连通分量来进行缩点 首先给出强连通分量的 ...

随机推荐

  1. IDEA集成的 Thrift 插件进行 thrift 编译

    IDEA集成的 Thrift 插件进行 thrift 编译 注意 1.thrift文件要放在source目录,才有compile选项 2.generator list的output path不用加包名 ...

  2. QT4.5.3移植到hi3536

    QT版本:qt-embedded-linux-opensource-src-4.5.31.添加交叉编译(1)copy qt-embedded-linux-opensource-src-4.5.3/mk ...

  3. JavaScript学习(三)

  4. laravel出现No hint path defined for [sudosu]的解决方法

    今天ytkah在部署laravel项目时出现了No hint path defined for [sudosu]的问题,大概意思是没有定义sudosu的提示路径,那我们找一下配置文件有没相关设置,看到 ...

  5. mysql开发总结

    一.理解MySQL基本概念 1.MySQL软件:MySQL实际上就是一软件,是一工具,是关系型数据库管理系统软件 2.MySQL数据库:就是按照数据结构来组织.存储和管理数据的仓库 3.MySQL数据 ...

  6. Mysql事务原理介绍

    事务 一个事务会涉及到大量的cpu计算和IO操作,这些操作被打包成一个执行单元,要么同时都完成,要么同时都不完成. 事务是一组原子性的sql命令或者说是一个独立的工作单元,如果数据库引擎能够成功的对数 ...

  7. Redis入门到高可用(九)——有序集合zset

    一.数据结构 集合与有序集合,列表与有序集合的对比 二.主要API zadd 将一个或多个 member 元素及其 score 值加入到有序集 key 当中. zrem 移除有序集 key 中的一个或 ...

  8. js计算常见操作

    如何实现数字相加 var a = 1, var b = 2, var c = a + b 这样c得出来的解果是12, 使用Number()函数可以解决这个问题,如下 var c = Number(a) ...

  9. 解决PHP5.6版本“No input file specified”的问题

    问题描述:使用TP框架做项目时,在启用REWRITE的伪静态功能的时候,首页可以访问,但是访问其它页面的时候,就提示:“No input file specified.”原因在于使用的PHP5.6是f ...

  10. lnmp/nginx系统真正有效的图片防盗链完整设置详解

    http://www.it300.com/article-15345.html 关于nginx防盗链的方法网上有很多教程,都可以用,但是我发现很多教程并不完整,所做的防盗链并不是真正的彻底的防盗链! ...