hdu 5398 动态树LCT
GCD Tree
Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 415 Accepted Submission(s): 172
You need to find a subset of the edges that forms a tree that includes every vertex, where the total weight of all the edges in the tree is maximized. Print the total weight of these edges.
For each test case, there is only one line contains one number n(1≤n≤105).
2
3
4
5
1
2
4
5
/*
hdu 5398 动态树LCT problem:
给你n个点,让你构成一棵树,边的权值为端点的最大公约数. 求树的最大边权和 solve:
枚举所有点,然后贪心的思想.如果n-1个点构成了一个最大的生成树,那么加入第n个点时我们应尽可能把它往
自己的约数上连(大->小).这样的话就在连接第二个约数的时候就会构成环,所以需要将 当前点-->当前约数点的链上
断开最小的一条边.
所以可以枚举 点和其约数,然后利用LCT维护链上面的最小边. 就这个想了很久,不知道应该吧边权保存在哪个端点
比较好,总感觉会有矛盾 = =||. 结果发现完全可以弄个点来替代边,将点上的值初始化无限大就行了T_T,好2... hhh-2016-08-21 16:43:04
*/
#pragma comment(linker,"/STACK:124000000,124000000")
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <vector>
#include <map>
#define lson ch[0]
#define rson ch[1]
#define ll long long
#define clr(a,b) memset(a,b,sizeof(a))
#define key_val ch[ch[root][1]][0]
using namespace std;
const int maxn = 100100;
const int INF = 0x3f3f3f3f; struct Node* null;
struct Node
{
Node* ch[2] ;
Node* fa;
int Min;
Node* Minnode;
int Size,lpos,rpos ;
int rev;
void newnode(int v)
{
Min = v;
Size = 1 ;
Minnode = this;
lpos = rpos = 0;
fa = ch[0] = ch[1] = null ;
rev = 0;
}
void ini()
{
Min = INF;
Size = 1 ;
Minnode = this;
lpos = rpos = 0;
fa = ch[0] = ch[1] = null ;
rev = 0;
}
void update_rev()
{
if(this == null)
return ;
swap(ch[0],ch[1]);
rev ^= 1;
} void push_up ()
{
Minnode = this;
if(ch[0] && ch[0]->Minnode->Min < Minnode->Min)
Minnode = ch[0]->Minnode;
if(ch[1] && ch[1]->Minnode->Min < Minnode->Min)
Minnode = ch[1]->Minnode;
} void push_down()
{
if(rev)
{
ch[0]->update_rev();
ch[1]->update_rev();
rev = 0;
}
} void link_child ( Node* o , int d )
{
ch[d] = o ;
o->fa = this ;
} int isroot()
{
return fa == null || this != fa->ch[0] && this != fa->ch[1] ;
}
void sign_down ()
{
if ( !isroot () ) fa->sign_down () ;
push_down () ;
}
void Rotate ( int d )
{
Node* p = fa ;
Node* gp = fa->fa ;
p->link_child ( ch[d] , !d ) ;
if ( !p->isroot () )
{
if ( gp->ch[0] == p ) gp->link_child ( this , 0 ) ;
else gp->link_child ( this , 1 ) ;
}
else fa = gp ;
link_child ( p , d ) ;
p->push_up () ;
} void splay ()
{
sign_down () ;
while ( !isroot () )
{
if ( fa->isroot () )
{
this == fa->ch[0] ? Rotate ( 1 ) : Rotate ( 0 ) ;
}
else
{
if ( fa == fa->fa->ch[0] )
{
this == fa->ch[0] ? fa->Rotate ( 1 ) : Rotate ( 0 ) ;
Rotate ( 1 ) ;
}
else
{
this == fa->ch[1] ? fa->Rotate ( 0 ) : Rotate ( 1 ) ;
Rotate ( 0 ) ;
}
}
}
push_up () ;
} void access()
{
Node* o = this ;
Node* x = null ;
while ( o != null )
{
o->splay () ;
o->link_child ( x , 1 ) ;
o->push_up () ;
x = o ;
o = o->fa ;
}
splay () ;
} void make_root()
{
access();
update_rev();
} void cut()
{
access();
ch[0]->fa = null;
ch[0] = null;
push_up();
}
Node* find_root ()
{
access () ;
Node* o = this ;
while ( o->ch[0] != null )
{
o->push_down () ;
o = o->ch[0] ;
}
return o ;
}
void cut(Node* to)
{
to->make_root();
cut();
} void link(Node* to)
{
to->make_root();
to->fa = this;
}
Node* query(Node* to)
{
to->make_root();
access();
return Minnode;
}
};
Node memory_pool[maxn*2];
Node* now;
Node* node[maxn]; void Clear()
{
now = memory_pool;
now->newnode(INF);
null = now ++;
null->Size = 0;
}
int dp[maxn];
vector<int > vec[maxn]; void init()
{
Clear();
for(int i = 100000; i >= 1; i--)
{
now->newnode(INF);
node[i] = now++;
for(int j = i+i; j <= 100000; j+=i)
{
vec[j].push_back(i);
}
}
int tans = 0;
dp[1] = 0;
now->newnode(INF),node[1] = now++;
for(int i =2; i <= 100000; i++)
{
// cout << "number:" << i<<endl;
now->newnode(vec[i][0]);
Node *tp = now++;
tp->link(node[i]);
tp->link(node[vec[i][0]]);
tp->lpos = i,tp->rpos = vec[i][0];
tans += vec[i][0];
// cout << tp->lpos <<" " <<tp->rpos <<" " <<tans <<endl;
for(int j = 1;j < vec[i].size();j++)
{
int to = vec[i][j];
// cout <<"to "<<to <<endl; Node* tMin = node[to]->query(node[i]);
// cout <<"Min:" << tMin->lpos <<" " <<tMin->rpos <<" "<<tMin->Min <<endl;
if(tMin->Min < to)
{
tans -= tMin->Min;
tMin->cut(node[tMin->lpos]);
// cout <<"cut1 ";
tMin->cut(node[tMin->rpos]);
// cout << "cut2" <<endl;
tMin->ini();
tMin->Min = to,tMin->lpos = i,tMin->rpos = to;
// cout << i <<" " <<to<<endl; tMin->link(node[i]);
// cout <<"link1 ";
tMin->link(node[to]);
// cout << "link2" <<endl;
tans += to;
}
}
dp[i] = tans ;
}
} int main()
{ // freopen("in.txt","r",stdin);
init();
int n;
while( scanf("%d",&n) !=EOF)
{
printf("%d\n",dp[n]);
}
return 0;
}
hdu 5398 动态树LCT的更多相关文章
- hdu 5002 (动态树lct)
Tree Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submi ...
- hdu 5314 动态树
Happy King Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Tot ...
- HDU 4718 The LCIS on the Tree (动态树LCT)
The LCIS on the Tree Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Oth ...
- 动态树LCT小结
最开始看动态树不知道找了多少资料,总感觉不能完全理解.但其实理解了就是那么一回事...动态树在某种意思上来说跟树链剖分很相似,都是为了解决序列问题,树链剖分由于树的形态是不变的,所以可以通过预处理节点 ...
- bzoj2049-洞穴勘测(动态树lct模板题)
Description 辉辉热衷于洞穴勘测.某天,他按照地图来到了一片被标记为JSZX的洞穴群地区.经过初步勘测,辉辉发现这片区域由n个洞穴(分别编号为1到n)以及若干通道组成,并且每条通道连接了恰好 ...
- [模板] 动态树/LCT
简介 LCT是一种数据结构, 可以维护树的动态加边, 删边, 维护链上信息(满足结合律), 单次操作时间复杂度 \(O(\log n)\).(不会证) 思想类似树链剖分, 因为splay可以换根, 用 ...
- 动态树LCT(Link-cut-tree)总结+模板题+各种题目
一.理解LCT的工作原理 先看一道例题: 让你维护一棵给定的树,需要支持下面两种操作: Change x val: 令x点的点权变为val Query x y: 计算x,y之间的唯一的最短路径的点 ...
- SPOJ OTOCI 动态树 LCT
SPOJ OTOCI 裸的动态树问题. 回顾一下我们对树的认识. 最初,它是一个连通的无向的无环的图,然后我们发现由一个根出发进行BFS 会出现层次分明的树状图形. 然后根据树的递归和层次性质,我们得 ...
- BZOJ 2002: [Hnoi2010]Bounce 弹飞绵羊 (动态树LCT)
2002: [Hnoi2010]Bounce 弹飞绵羊 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 2843 Solved: 1519[Submi ...
随机推荐
- 敏捷冲刺报告--Day5
敏捷冲刺报告--Day5 情况简介 GUI框架重写, 添加功能 任务进度 赵坤: 后端爬虫bug修复 李世钰: GUI编写 黄亦薇:更新sprint backlog.编写每日报告 王成科:召集小组成员 ...
- 【iOS】Swift if let 和 if var
if let unwrappedOptional = postDict { print("The optional has a value! It's \(unwrappedOptional ...
- Flask 页面缓存逻辑,jinja2 过滤器,测试器
回调接入点-页面缓存逻辑 from flask import Flask,request,render_template from werkzeug.contrib.cache import Simp ...
- highcharts 具体参数详解
<script type="text/javascript" src="js/jquery.min.js"></script> < ...
- mysql5.5中datetime默认值不能为NOW或者CURRENT_TIMESTAMP,用触发器解决
mysql5.6及以上的版本datatime默认值可以为CURRENT_TIMESTAMP或者NOW 那我们要用的是mysql5.5及以下版本呢? 请看代码 delimiter // DROP TRI ...
- React 深入系列2:组件分类
文:徐超,<React进阶之路>作者 授权发布,转载请注明作者及出处 React 深入系列2:组件分类 React 深入系列,深入讲解了React中的重点概念.特性和模式等,旨在帮助大家加 ...
- Python内置函数(5)——pow
英文文档: pow(x, y[, z]) Return x to the power y; if z is present, return x to the power y, modulo z (co ...
- js new到底干了什么,new的意义是什么?
学过JS的都知道 创建对象可以这样 var obj=new Object(); var obj=new Function(); 用内置的函数对象来构造对象 还可以这样自定义函数 function te ...
- 新概念英语(1-29)Come in, Amy.
How must Amy clean the floor? A:Come in, Amy. Shut the door, please. This bedroom's very untidy. B:W ...
- notepad++运行Python
1.打开notepad++的菜单栏,点击run 2.输入cmd /k python "$(FULL_CURRENT_PATH)" & PAUSE & EXIT 3. ...