题目大意:

求最小方差生成树.N<=100,M<=2000,Ci<=100

题解:


首先我们知道这么一个东西:

  • 一些数和另一个数的差的平方之和的最小值在这个数是这些数的平均值时取得

所以我们可以枚举这个平均数,然后计算所有边与该值的差的平方

然后扔下去跑一个最小生成树

然后我们通过枚举这个平均数发现这个平均数和答案的对应函数的图像是一个波形函数

所以我们可以直接在这个波形图像上找函数最低点:

相应的就有

  • 爬山算法
  • 模拟退火

两种算法


所以我们可以先在全局用模拟退火然后在局部用爬山算法。

然而还是每三组数据就Wa一次

然后发现这样的话极限数据只需要0.8s,还有1.2s可以用

所以可以在全局再三分找函数最低点.

随机化左右端点然后再三分.

随机化22次端点极限数据可以跑到1.3s

#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
inline void read(int &x){
x=0;char ch;bool flag = false;
while(ch=getchar(),ch<'!');if(ch == '-') ch=getchar(),flag = true;
while(x=10*x+ch-'0',ch=getchar(),ch>'!');if(flag) x=-x;
}
const int maxn = 128;
const int maxm = 2048;
const double eps = 1e-10;
const double det = 0.99;
struct Node{
int u,v,bed;
double d;
bool friend operator < (const Node &a,const Node &b){
return a.d < b.d;
}
}G[maxm];
int fa[maxn];
inline int find(int x){return x == fa[x] ? x : fa[x] = find(fa[x]);}
int n,m;
inline double work(double mid){
for(int i=1;i<=n;++i) fa[i] = i;
for(int i=1;i<=m;++i) G[i].d = (G[i].bed - mid)*(G[i].bed - mid);
sort(G+1,G+m+1);
int cnt = 0;double sum = 0;
for(int i=1;i<=m;++i){
int x = find(G[i].u);
int y = find(G[i].v);
if(x != y){
sum += G[i].d;
fa[x] = y;
if(++cnt == n-1) break;
}
}
return sqrt(sum/(n-1));
}
inline double ran(){
return (1.0*rand())/1000.0;
}
double ans = 1e10,ans_p;
inline double f(double mid){
double x = work(mid);
if(ans > x){
ans = x;
ans_p = mid;
}return x;
}
int main(){
srand(2333);
read(n);read(m);
int minn = 0x7f7f7f7f,maxx = -0x7f7f7f7f;
for(int i=1;i<=m;++i){
read(G[i].u);
read(G[i].v);
read(G[i].bed);
minn = min(minn,G[i].bed);
maxx = max(maxx,G[i].bed);
}
double l = minn,r = maxx,nx,t,val_nx;
double T = 50.0,nw = (l+r)/2,val_nw;
while(T > eps){
nx = ( rand() % 2 == 0 ? -1 : 1)*ran()*T;
t = val_nw - (val_nx = f(nx));
if(t > 0 || exp(t/T) >= ran() ){
nw = nx;
val_nw = val_nx;
}
T *= det;
}
while(r-l > eps){
double midx = (l+l+r)/3;
double midy = (l+r+r)/3;
double x = f(midx);
double y = f(midy);
if(x < y) r = midy;
else l = midx;
}
ans = min(ans,f(l));
int num = 22;
while(num--){
double l = minn + rand()*ran()*0.1;
double r = maxx - rand()*ran()*0.1;
if(l > r) swap(l,r);
while(r-l > eps){
double midx = (l+l+r)/3;
double midy = (l+r+r)/3;
double x = f(midx);
double y = f(midy);
if(x < y) r = midy;
else l = midx;
}ans = min(ans,f(l));
}
printf("%.4lf\n",ans);
getchar();getchar();
return 0;
}

bzoj 3754: Tree之最小方差树 模拟退火+随机三分的更多相关文章

  1. BZOJ 3754 Tree之最小方差树 MST

    Description Wayne 在玩儿一个很有趣的游戏.在游戏中,Wayne 建造了N 个城市,现在他想在这些城市间修一些公路,当然并不是任意两个城市间都能修,为了道路系统的美观,一共只有M 对城 ...

  2. BZOJ 3754 Tree之最小方差树

    枚举平均数. mdzz编译器. #include<iostream> #include<cstdio> #include<cstring> #include< ...

  3. [BZOJ3754]Tree之最小方差树

    3754: Tree之最小方差树 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 402  Solved: 152[Submit][Status][Di ...

  4. [BZOJ3080]Minimum Variance Spanning Tree/[BZOJ3754]Tree之最小方差树

    [BZOJ3080]Minimum Variance Spanning Tree/[BZOJ3754]Tree之最小方差树 题目大意: 给定一个\(n(n\le50)\)个点,\(m(m\le1000 ...

  5. 【bzoj3754】Tree之最小方差树 最小生成树

    题目描述 给出一张无向图,求它的一棵生成树,使得选出的所有边的方差最小.输出这个最小方差. 输入 第一行两个正整数N,M 接下来M行,每行三个正整数Ui,Vi,Ci N<=100,M<=2 ...

  6. 【BZOJ 3754】Tree之最小方差树

    http://www.lydsy.com/JudgeOnline/problem.php?id=3754 核心思想:暴力枚举所有可能的平均数,对每个平均数排序后Kruskal. 正确的答案一定是最小的 ...

  7. 【BZOJ 3754】: Tree之最小方差树

    题目链接: TP 题解: 都是骗子233,我还以为是什么神奇的算法. 由于边权的范围很小,最小生成树和最大生成树之间的总和差不会太大,所以可以枚举边权和,再直接根据方差建最小生成树,每次更新答案即可. ...

  8. 【枚举】【最小生成树】【kruscal】bzoj3754 Tree之最小方差树

    发现,若使方差最小,则使Σ(wi-平均数)2最小即可. 因为权值的范围很小,所以我们可以枚举这个平均数,每次把边权赋成(wi-平均数)2,做kruscal. 但是,我们怎么知道枚举出来的平均数是不是恰 ...

  9. bzoj3754 Tree之最小方差树 最小生成树+推性质

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3754 题解 感觉这个思路挺神仙的. 后悔没有好好观察题目的数据范围,一直把 \(n\) 和 \ ...

随机推荐

  1. Mark指针的指针(**)和链表使用(*&)

    利用二级指针删除单向链表 彻底理解链表中为何使用指针的指针或者指针的引用 详解C++指针的指针和指针的引用

  2. Myeclipse 选中高亮

    打开显示功能 选择Windows->Preferences->Java-> Editor-> Mark Occurrences ,勾选选项.这时,当你单击一个元素的时候,代码中 ...

  3. 记录-MySQL中的事件调度Event Scheduler

    下面是自己的实例 /*查询event是否开启(查询结果Off为关闭 On为开启)*/show variables like '%sche%'; /*开启/关闭命令(1开启--0关闭)*/set glo ...

  4. java拾遗3----XML解析(三) StAX PULL解析

    使用PULL方式解析XML: Pull是STAX的一个实现 StAX是The Streaming API for XML的缩写,一种利用拉模式解析(pull-parsing)XML文档的API StA ...

  5. Android5.0以后版本把应用移动到SD或者TF卡的方法

    由于手机内存较小,才8G,用的时间一久,内部存储就满了,天天删垃圾,WIFI还老断线,终于忍无可忍了,决定把应用移动到SD卡,实践下来,只有少部分App默认支持移动到SD卡,大部分程序不支持只能装在内 ...

  6. java基础入门1到100的奇数求和

    /* Name:1-100所有奇数求和的程序 Power by Stuart Date:2015-4-23 */ public class DateTest01{ public static void ...

  7. Java并发之ArrayBlockingQueue

    ArrayBlockingQueue是一个由数组支持的有界阻塞队列.此队列按 FIFO(先进先出)原则对元素进行排序.队列的头部是在队列中存在时间最长的元素.队列的尾部是在队列中存在时间最短的元素.新 ...

  8. tomcat异常处理经验汇总

    1.Https: Feb 21, 2018 5:22:02 PM org.apache.coyote.AbstractProtocol initSEVERE: Failed to initialize ...

  9. Python 3 mysql 库操作

    Python 3 mysql 库操作 一.基础相关知识 MySQL数据库基本操作知识储备 数据库服务器:一台计算机(对内存要求比较高) 数据库管理系统:如mysql,是一个软件 数据库:oldboy_ ...

  10. vim设置tab为4空格

    vim的最基础设置 vim的设置需要编辑~/.vimrc文件,更改已有设置或者在后面添加相应的设置. 设置tab为4字符 # ts: tabstop set ts=4 将tab展开为空格 # expa ...