题目链接:https://www.luogu.org/problemnew/show/P4180

这个题卡树剖。记得开O2。

这个题inf要到1e18。

定理:次小生成树和最小生成树差距只有在一条边上

非严格次小生成树:枚举每一条不在最小生成树上的边,加入到最小生成树中构成一个环。删去这个环上的最大值。(此最大值有可能与加入生成树中的边相等,故为非严格次小生成树。)重复此操作取min,得到次小生成树。(基于kruskal实现。)

严格次小生成树:与非严格次小生成树类似,不同在于为了避免删去环上的最大值等于加入生成树中的边,需要记录次最大值。恶心点所在。

于是维护最大值和次小值又成了一道数据结构题。

树剖:剖MST,查询加进来的边的两端点编号 = =

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ll long long
#define lson l, mid, rt<<1
#define rson mid+1, r, rt<<1|1
using namespace std;
const int maxn = 300010;
const ll inf = 1e18;
inline ll read()
{
char ch = getchar(); ll u = 0, f = 1;
while (!isdigit(ch)) {if (ch == '-')f = -1; ch = getchar();}
while (isdigit(ch)) {u = u * 10 + ch - 48; ch = getchar();}return u * f;
}
ll n, m, fa[maxn], deep[maxn], size[maxn], son[maxn], top[maxn], seg[maxn], rev[maxn], W[maxn], ans = inf, num;
ll father[maxn], mstans, mstcnt;
bool vis[maxn];
struct EDG{
ll u, v, w;
}G[maxn];//mst
struct edge{
ll to, next, len;
}e[maxn<<2];
ll head[maxn], cnt;
void add(ll u, ll v, ll w)
{
e[++cnt].len = w; e[cnt].next = head[u]; e[cnt].to = v; head[u] = cnt;
e[++cnt].len = w; e[cnt].next = head[v]; e[cnt].to = u; head[v] = cnt;
}
//kruskal
bool cmp(EDG a, EDG b)
{
return a.w < b.w;
}
ll find(ll x)
{
return father[x] == x ? x : father[x] = find(father[x]);
}
void init()
{
for(ll i = 1; i <= n; i++) father[i] = i;
sort(G+1, G+1+m, cmp);
}
void kruskal()
{
init();
for(ll i = 1; i <= m; i++)
{
if(mstcnt == n-1) break;
ll x = find(G[i].u), y = find(G[i].v);
if(x != y)
{
mstans += G[i].w;
mstcnt++;
vis[i] = 1;
add(G[i].u, G[i].v, G[i].w);
father[x] = y;
}
}
}
//Segment_Tree
bool maxcmp(ll a, ll b)
{
return a > b;
}
ll get_sec(ll a, ll b, ll c, ll d)
{
ll z[5] = {a, b, c, d};
sort(z, z+4, maxcmp);
for(ll i = 1; i <= 3; i++)
{
if(z[i] != z[0]) return z[i];
}
}
struct Segment_Tree{
ll fir, sec;
}tree[maxn<<2];
void PushUPfir(ll rt)
{
tree[rt].fir = max(tree[rt<<1].fir, tree[rt<<1|1].fir);
}
void PushUPsec(ll rt)
{
tree[rt].sec = get_sec(tree[rt<<1].fir, tree[rt<<1|1].fir, tree[rt<<1].sec, tree[rt<<1|1].sec);
}
void build(ll l, ll r, ll rt)
{
if(l == r)
{
tree[rt].fir = rev[l];
return;
}
ll mid = (l + r) >> 1;
build(lson);
build(rson);
PushUPfir(rt);
PushUPsec(rt);
}
Segment_Tree query(ll left, ll right, ll l, ll r, ll rt)
{
Segment_Tree max1 = {-inf,-inf}, max2 = {-inf,-inf};
if(left <= l && r <= right)
{
return (Segment_Tree){tree[rt].fir, tree[rt].sec};
}
ll mid = (l + r) >> 1;
if(left <= mid) max1 = query(left, right, lson);
if(right > mid) max2 = query(left, right, rson);
return (Segment_Tree) {max(max1.fir, max2.fir), get_sec(max1.fir, max1.sec, max2.fir, max2.sec)};
}
//Tree_cut
void dfs1(ll u, ll f, ll d)
{
ll maxson = -1;
size[u] = 1;
deep[u] = d;
fa[u] = f;
for(ll i = head[u]; i != -1; i = e[i].next)
{
ll v = e[i].to;
if(f != v)
{
W[v] = W[u] + e[i].len;
dfs1(v, u, d+1);
size[u] += size[v];
if(maxson < size[v]) maxson = size[v], son[u] = v;
}
}
}
void dfs2(ll u, ll t)
{
seg[u] = ++num;
rev[num] = W[u] - W[fa[u]];//前缀和边权上点权
//rev[num] = node[u];
top[u] = t;
if(!son[u]) return;
dfs2(son[u], t);
for(ll i = head[u]; i != -1; i = e[i].next)
{
ll v = e[i].to;
if(v == fa[u] || v == son[u]) continue;
dfs2(v, v);
}
}
ll LCA(ll x, ll y, ll d)//当前边的权值
{
ll res = -inf;
while(top[x] != top[y])
{
if(deep[top[x]] < deep[top[y]]) swap(x, y);
Segment_Tree temp1 = query(seg[top[x]], seg[x], 1, n, 1);
x = fa[top[x]];
res = max(res, (temp1.fir == d) ? temp1.sec : temp1.fir);
}
if(deep[x] < deep[y]) swap(x, y);
Segment_Tree temp2 = query(seg[y] + 1, seg[x], 1, n, 1);
return res = max(res, (temp2.fir == d) ? temp2.sec : temp2.fir);
}
int main()
{
memset(head, -1, sizeof(head));
n = read(); m = read(); //scanf("%lld%lld",&n,&m);
for(ll i = 1; i <= m; i++) {G[i].u = read(); G[i].v = read(); G[i].w = read();}//scanf("%lld%lld%lld",&G[i].u,&G[i].v,&G[i].w);
kruskal();
dfs1(1, 0, 1); dfs2(1, 1);
build(1, n, 1);
for(ll i = 1; i <= m; i++)
{
if(vis[i]) continue;
ll temp = mstans + G[i].w - LCA(G[i].u, G[i].v, G[i].w);
if(ans > temp && temp != mstans + G[i].w && temp > mstans) ans = temp;
}
printf("%lld",ans);
return 0;
}

【luogu P4180 严格次小生成树[BJWC2010]】 模板的更多相关文章

  1. P4180 严格次小生成树[BJWC2010] Kruskal,倍增

    题目链接\(Click\) \(Here\). 题意就是要求一个图的严格次小生成树.以前被题面吓到了没敢做,写了一下发现并不难. 既然要考虑次小我们就先考虑最小.可以感性理解到一定有一种次小生成树,可 ...

  2. Luogu P4180 【模板】严格次小生成树[BJWC2010]

    P4180 [模板]严格次小生成树[BJWC2010] 题意 题目描述 小\(C\)最近学了很多最小生成树的算法,\(Prim\)算法.\(Kurskal\)算法.消圈算法等等.正当小\(C\)洋洋得 ...

  3. P4180 【模板】严格次小生成树[BJWC2010]

    P4180 [模板]严格次小生成树[BJWC2010] 倍增(LCA)+最小生成树 施工队挖断学校光缆导致断网1天(大雾) 考虑直接枚举不在最小生成树上的边.但是边权可能与最小生成树上的边相等,这样删 ...

  4. 【洛谷】4180:【模板】严格次小生成树[BJWC2010]【链剖】【线段树维护最大、严格次大值】

    P4180 [模板]严格次小生成树[BJWC2010] 题目描述 小C最近学了很多最小生成树的算法,Prim算法.Kurskal算法.消圈算法等等.正当小C洋洋得意之时,小P又来泼小C冷水了.小P说, ...

  5. 「LuoguP4180」 【模板】严格次小生成树[BJWC2010](倍增 LCA Kruscal

    题目描述 小C最近学了很多最小生成树的算法,Prim算法.Kurskal算法.消圈算法等等.正当小C洋洋得意之时,小P又来泼小C冷水了.小P说,让小C求出一个无向图的次小生成树,而且这个次小生成树还得 ...

  6. 洛谷 P4180 【模板】严格次小生成树[BJWC2010]【次小生成树】

    严格次小生成树模板 算法流程: 先用克鲁斯卡尔求最小生成树,然后给这个最小生成树树剖一下,维护边权转点权,维护最大值和严格次大值. 然后枚举没有被选入最小生成树的边,在最小生成树上查一下这条边的两端点 ...

  7. 【洛谷 P4180】【模板】严格次小生成树[BJWC2010](倍增)

    题目链接 题意如题. 这题作为我们KS图论的T4,我直接打了个很暴力的暴力,骗了20分.. 当然,我们KS里的数据范围远不及这题. 这题我debug了整整一个晚上还没debug出来,第二天早上眼前一亮 ...

  8. 【【模板】严格次小生成树[BJWC2010]】

    树上的路径怎么能没有树剖 显然,次小生成树和最小生成树只在一条边上有差距,于是我们就可以枚举这一条边,将所有边加入最小生成树,之后再来从这些并不是那么小的生成树中找到那个最小的 我们往最小生成树里加入 ...

  9. POJ -1679(次小生成树)模板

    The Unique MST Time Limit: 1000MS   Memory Limit: 10000K Total Submissions:34617   Accepted: 12637 D ...

随机推荐

  1. poj 3104 dring 二分

    Drying Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 7684   Accepted: 1967 Descriptio ...

  2. Windows的图形设备接口与Windows绘图

    本次学习目标 理解DC, 映像模式, 坐标系统, 窗口和视口; 学习获取绘图工具(画笔/画刷)的句柄, 设置颜色, 能定义映像模式; 会使用常用的绘图函数. 编写程序: 在屏幕上出现一个圆心沿正弦曲线 ...

  3. python学习之老男孩python全栈第九期_day027知识点总结——反射、类的内置方法

    一. 反射 ''' # isinstance class A:pass class B(A):pass a = A() print(isinstance(a,A)) # 判断对象和类的关系 print ...

  4. CSS实现文字两端对齐

    最近的项目遇到了这样的需求:(要求标题部分不管文字多少,都必须两端对齐) 如下图: 当时也没有多想直接使用‘ ’进行代替,毕竟产品同学想快一点看到效果,不敢怠慢!不过到第二个页面就傻眼了. 如图: 这 ...

  5. Struts2详讲

    一 概述 1.什么是Struts2? Struts2是一个在WebWork框架基础上发展起来开源MVC框架. 2.StrutsPrepareAndExecuteFilter StrutsPrepare ...

  6. System.Media.Color与System.Drawing.Color转换方法

    public static System.Media.Color GetMediaColorFromDrawingColor(System.Drawing.Color color) {      re ...

  7. VC学习笔记---ATL MFC CLR三个库的区别

    MFC.ATL和CLR是VC2005内置的三大库,涵盖了Windows的各种开发方法和开发应用.当然关于C++开发的库不止这三个,不过这三个是微软推荐. 从编程所处层次而言,WIN32为最底层,其次是 ...

  8. Android 性能优化的方面方面都在这儿

    又到周六了,鸿洋的不定期的周六放送又来了~~这次来谈谈性能优化吧.大家在工作中或多或少都会拿自家的应用和竞品app做比对,不可避免的需要做一些app性能优化的活.很多时候可能是策略上的调整,不过还是有 ...

  9. 小程序-报错 xxx is not defined (已解决)

    小程序-报错 xxx is not defined (已解决) 问题情境: 这样一段代码,微信的小程序报错 is not defined 我 wxml 想这样调用 //wxml 代码 <view ...

  10. gulp 在 angular 项目中的使用

    gulp 在 angular 项目中的使用 keyword:gulp,angularjs,ng,ngAnnotate,jshint,gulpfile 最后附完整简洁的ng项目gulpfile.js 准 ...