2017 多校2 I Curse Myself(第k小生成树)

题目: 给一张带权无向连通图,该图的任意一条边最多只会经过一个简单环,定义\(V(k)为第k小生成树的权值和\),求出\(\sum_{k=1}^{K}k \cdot V(k) mod 2^{32}\)

思路: 比赛的时候看了一眼,没有看清楚是仙人掌那句话,觉得好难啊

看完题解之后觉得就算看清了还是过不了嘛。

直接上题解

由于图是一个仙人掌,所以显然对于图上的每一个环都需要从环上取出一条边删掉。所以问题就变为有 M 个集合,每个集合里面都有一堆数字,要从每个集合中选择一个恰好一个数加起来。求所有的这样的和中,前 K 大的是哪些。这就是一个经典问题了。

对所有集合两个两个进行合并,设当前合并的集合是 A 和BB,合并的过程中用堆来求出当前 \(A_{i} + B_{j}\)

​ 这样的复杂度看起来是\(\mathcal{O}(MK \log K)\),但如果合并的时候保证堆内的元素个数是新集合里的元素个数,设每个集合的大小分别为 \(m_{0}, m_{1}, \cdots, m_{M-1}\)

​​ ,则复杂度为 \(\mathcal{O}(\sum{K \log{m_{i}}}) = \mathcal{O}(K \log{\prod{m_i}})\)。当 \(m_{i}\)

​​ 都相等时取得最大值 \(\mathcal{O}\left(MK \log{\frac{\sum{m_i}}{M}}\right)\),所以实际复杂度为 \(\mathcal{O}(MK)\)。

就照着题解敲了一遍,该优化的地方优化了,

找环那里用时间戳的方法一直没写对,所以换成直接用树再加边的方法去暴力找环了。

合并有序表的问题 刘汝佳的书上有讲

#include<bits/stdc++.h>
#define LL long long
#define P pair<int,int> using namespace std;
const LL mod = (1LL<<32);
const int N = 1200;
const int M = 4 * N;
int n, m, K;
int read(){
int x = 0;
char c = getchar();
while(c < '0' || c > '9') c = getchar();
while(c >= '0' && c <= '9') x = x * 10 + c - 48,c = getchar();
return x;
}
struct edge{
int v,w,nxt;
edge(){};
}e[M];
struct Edge{
int u, v, w;
Edge(int u,int v,int w):u(u),v(v),w(w){};
Edge(){};
};
vector<Edge> other;
int head[N],EN;
int pa[N];
int Find(int x){
return pa[x] == x?x:pa[x] = Find(pa[x]);
}
int A[100010],B[100010],ca,cb;
void init(){
ca = cb = EN = 0;
memset(head,-1,sizeof(head));
other.clear();
for(int i = 1;i <= n;i++) pa[i] = i;
}
void add(int u,int v,int w){
e[EN].v = v,e[EN].w = w,e[EN].nxt = head[u];
head[u] = EN++;
}
bool dfs(int u,int f,int ed){
if(u == ed) return true;
for(int i = head[u];~i;i = e[i].nxt){
if(i == f) continue;
if(dfs(e[i].v,i ^ 1,ed)){
B[cb++] = e[i].w;
return true;
}
}
return false;
}
bool cmp(int x,int y){
return x > y;
}
void Merge(){
if(ca > cb) swap(ca,cb),swap(A,B);
priority_queue<P> q;
for(int i = 0;i < ca;i++) q.push(P(A[i]+B[0],0));
int siz = min(K, ca * cb);
ca = 0;
for(int i = 0;i < siz;i++){
P cur = q.top();q.pop();
A[ca++] = cur.first;
if(cur.second + 1 < cb) q.push(P(cur.first - B[cur.second] + B[cur.second+1],cur.second + 1));
}
}
int main(){
int cas = 1, u, v, w;
while(scanf("%d%d",&n,&m)==2){
init();
int total = 0;
for(int i = 0;i < m;i++){
u = read(),v = read(),w = read();
if(Find(u) != Find(v)){
add(u, v, w);
add(v, u, w);
pa[Find(u)] = Find(v);
}else {
other.push_back(Edge(u,v,w));
}
total += w;
}
K = read();
int first = 1;
for(auto E:other){
dfs(E.u,-1,E.v);
B[cb++] = E.w;
sort(B, B + cb,cmp);
if(first) first = 0,swap(A,B),swap(ca,cb);
else Merge();
cb = 0;
}
LL ans = 0;
for(int i = 0;i < ca;i++) {
ans = (ans + 1LL * (i + 1) * (total - A[i]))%mod;
}
if(!ca) ans = total;
printf("Case #%d: %lld\n",cas++,ans);
}
return 0;
}

2017 多校1 I Curse Myself的更多相关文章

  1. 2017 多校5 hdu 6093 Rikka with Number

    2017 多校5 Rikka with Number(数学 + 数位dp) 题意: 统计\([L,R]\)内 有多少数字 满足在某个\(d(d>=2)\)进制下是\(d\)的全排列的 \(1 & ...

  2. 2017 多校5 Rikka with String

    2017 多校5 Rikka with String(ac自动机+dp) 题意: Yuta has \(n\) \(01\) strings \(s_i\), and he wants to know ...

  3. 2017 多校4 Wavel Sequence

    2017 多校4 Wavel Sequence 题意: Formally, he defines a sequence \(a_1,a_2,...,a_n\) as ''wavel'' if and ...

  4. 2017 多校4 Security Check

    2017 多校4 Security Check 题意: 有\(A_i\)和\(B_i\)两个长度为\(n\)的队列过安检,当\(|A_i-B_j|>K\)的时候, \(A_i和B_j\)是可以同 ...

  5. 2017 多校3 hdu 6061 RXD and functions

    2017 多校3 hdu 6061 RXD and functions(FFT) 题意: 给一个函数\(f(x)=\sum_{i=0}^{n}c_i \cdot x^{i}\) 求\(g(x) = f ...

  6. 2017 多校2 hdu 6053 TrickGCD

    2017 多校2 hdu 6053 TrickGCD 题目: You are given an array \(A\) , and Zhu wants to know there are how ma ...

  7. hdu6136[模拟+优先队列] 2017多校8

    有点麻烦.. /*hdu6136[模拟+优先队列] 2017多校8*/ #include <bits/stdc++.h> using namespace std; typedef long ...

  8. hdu6134[莫比乌斯反演] 2017多校8

    /*hdu6134[莫比乌斯反演] 2017多校8*/ #include <bits/stdc++.h> using namespace std; typedef long long LL ...

  9. hdu6098[RMQ+筛法] 2017多校6

    /*hdu6098[RMQ+筛法] 2017多校6*/ #include <bits/stdc++.h> using namespace std; ][], len[], a[]; voi ...

随机推荐

  1. poj_1284_Primitive root

    We say that integer x, 0 < x < p, is a primitive root modulo odd prime p if and only if the se ...

  2. Spring Cloud 升级最新 Finchley 版本,踩坑指南!

    https://blog.csdn.net/youanyyou/article/details/81530240 Spring Cloud 升级最新 Finchley 版本,踩了所有的坑! 2018年 ...

  3. JDK5 新特性

     JDK5新特性目录导航: 自动拆装箱 Foreach 静态导入 可变参数 Var args 枚举 格式化输出 泛型 ProcessBuilder 内省 线程并发库(JUC) 监控和管理虚拟机 元数据 ...

  4. Laravel系列之环境搭建 — VirtualBox+Vagrant+Homestead

    一.为啥需要搭建环境 为了解决环境不统一问题,所以要搭建这么个玩意儿 二.步骤  Laravel对环境有所要求(不使用Homestead情况下),具体参考官网 使用Homestead步骤 1. Hom ...

  5. 关于 PHP 程序员技术职业生涯规划

    原文地址:http://rango.swoole.com/archives/570 看到很多 PHP 程序员职业规划的文章,都是直接上来就提 Linux.PHP.MySQL.Nginx.Redis.M ...

  6. 数据库DDL

    自己对数据库的整理,也是对自己知识的梳理 SQL ( Structure query language ) 结构化查询语言 SQL语言分为4个部分 1.DDL(Data Definition Lang ...

  7. Python 编码格式的使用

    编码史 ASCII > Unicode > UTF-8 Unicode支持多语言,UTF-8自动转换长短细节节省空间 在计算机内存中,统一使用Unicode编码,当需要保存到硬盘或者需要传 ...

  8. Poweroj:2425-跳台阶(经典递推)

    题目链接:https://www.oj.swust.edu.cn/problem/show/2425 跳台阶 Edit Manage Data Rejudge Time Limit: 1000 MS ...

  9. 单源次短路径:poj:3255-Roadblocks

    Roadblocks Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 17521 Accepted: 6167 Descripti ...

  10. 安装python 第三方库遇到的安装问题 microsoft visual studio c++ 10.0 is required,Could not find function xmlCheckVersion in library libxml2. Is libxml2 installed?

    问题一: microsoft visual studio c++ 10.0 is required 安装scrapy时候出现需要vc c++ 10,有时安装其他也会有. 解决方法:安装vc 2010, ...