loj2509 hnoi2018排列
题意:对于a数组,求它的一个合法排列的最大权值。合法排列:对于任意j,k,如果a[p[j]]=p[k],那么k<j。
权值:sigma(a[p[i]]*i)。n<=50W。
标程:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll read()
{
ll x=,f=;char ch=getchar();
while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<='') x=(x<<)+(x<<)+ch-'',ch=getchar();
return x*f;
}
const int N=;
vector<ll> vec[N];
ll ans,sum[N],he[N],cnt,head[N],n,a[N],fa[N],w[N],sz[N],cn,vis[N],tail[N],f[N];
int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
struct _node{ll sum,he,sz,id;_node(ll A,ll B,ll C,ll D){sum=A;he=B;sz=C;id=D;}};
struct cmp{
bool operator () (const _node &A,const _node &B)
{return A.he*B.sz>B.he*A.sz||A.he*B.sz==B.he*A.sz&&A.sz<B.sz;}
};
priority_queue<_node,vector<_node>,cmp> q;
int main()
{
n=read();
for (int i=;i<=n;i++) f[i]=i;
for (int i=;i<=n;i++)
{
a[i]=read();
if (<a[i]&&a[i]<=n) vec[a[i]].push_back(i),fa[i]=a[i]; else fa[i]=-;
}
for (int i=;i<=n;i++) w[i]=read(),q.push(_node(sum[i]=w[i],he[i]=w[i],sz[i]=,i));
while (!q.empty())
{
int x=q.top().id,fx=fa[find(x)];q.pop();
if (vis[x]) continue; vis[x]=; //dijkstra的思想,肯定先访问最后一次合并过的点,其他过去版本直接continue,这样就不用再记录一个del的堆。
if (fx==-)
{
ans+=sum[x]+cn*he[x];
for (int i=;i<vec[x].size();i++) fa[vec[x][i]]=-;
cn+=sz[x];
}
else {
if (vis[fx]) return puts("-1"),;
for (int i=;i<vec[x].size();i++) f[vec[x][i]]=find(x);
sum[fx]+=sum[x]+he[x]*sz[fx];he[fx]+=he[x];sz[fx]+=sz[x];
q.push(_node(sum[fx],he[fx],sz[fx],fx));
}
}
printf("%lld\n",ans);
return ;
}
易错点:1.居然碰到了yhx钦定的最难调错误没有之一,记!
return A.he*B.sz>B.he*A.sz||A.he*B.sz==B.he*A.sz&&A.sz<B.sz;
如果不判定相等的情况就不一定取到最后一个。
2.判断无解:vis表示已经被合并/删除的节点,重新连爸爸后爸爸应该是没有被删除的,如果vis[fa]=1,那么必然矛盾。
3.更改父亲的操作如果用vector暴力加,时间复杂度会到O(n^2)。用并查集保存同父亲的点是最快的做法。
题解:堆+并查集+建树+贪心
做法好神。如果a[p[j]]=p[k],那么k<j:也就是说比如a[1]=3,那么在排列中3一定在1前面。对于1<=a[i]<=n的点,连边a[i]->i,表示先取a[i],再取i。那么就形成了一棵树,如果有环必然无解。
这棵树肯定是每次取一个没有父亲的点作为p[i]。基于贪心,i越小,选越小的w[i]更优。
因此我们每次用堆/set维护权值最小的点,如果它没有父亲肯定直接取走,反之和其父亲合并,表示如果取走父亲后接下来肯定就取它。
合并之后,该点的儿子都连边向它父亲,也就是说fa[son[x]]=fa[x],可以用并查集维护。这样这个点的权值用sigma/size来代替。
可以证明:1.比较两个点的sigma/size就相当于比较它们sigma(i*w[p[i]])的权值。
2.对于同一个点,sigma/size随着合并不严格单调减。(设he1/sz1<he2/sz2,那么1向2合并,必然有he2/sz2>(he1+he2)/(sz1+sz2)。化简he2/sz2>(he1+he2)/(sz1+sz2),则he1*sz2<he2*sz1,同假设成立)
时间复杂度O(nlogn+na(n))。
loj2509 hnoi2018排列的更多相关文章
- 【BZOJ5289】[HNOI2018]排列(贪心)
[BZOJ5289][HNOI2018]排列(贪心) 题面 BZOJ 洛谷 题解 这个限制看起来不知道在干什么,其实就是找到所有排列\(p\)中,\(p_k=x\),那么\(k<j\),其中\( ...
- 5289: [Hnoi2018]排列
5289: [Hnoi2018]排列 链接 分析: 首先将题意转化一下:每个点向a[i]连一条边,构成了一个以0为根节点的树,要求选一个拓扑序,点x是拓扑序中的第i个,那么价值是i*w[x].让价值最 ...
- bzoj 5289: [Hnoi2018]排列
Description Solution 首先注意到实际上约束关系构成了一棵树 考虑这个排列 \(p\),编号为 \(a[i]\) 的出现了,\(i\) 才可以出现 那么如果连边 \((a[i],i) ...
- [HNOI2018]排列
Description: 给定 \(n\) 个整数 \(a_1, a_2, \dots, a_n, 0 \le a_i \le n\),以及 \(n\) 个整数 \(w_1, w_2, \dots, ...
- [HNOI2018]排列[堆]
题意 给定一棵树,每个点有点权,第 \(i\) 个点被删除的代价为 \(w_{p[i]}\times i\) ,问最小代价是多少. 分析 与国王游戏一题类似. 容易发现权值最小的点在其父亲选择后就会立 ...
- BZOJ5289: [Hnoi2018]排列
传送门 第一步转化,令 \(q[p[i]]=i\),那么题目变成: 有一些 \(q[a[i]]<q[i]\) 的限制,\(q\) 必须为排列,求 \(max(\sum_{i=1}^{n}w[i] ...
- BZOJ.5289.[AHOI/HNOI2018]排列(贪心 heap)
BZOJ LOJ 洛谷 \(Kelin\)写的挺清楚的... 要求如果\(a_{p_j}=p_k\),\(k\lt j\),可以理解为\(k\)要在\(j\)之前选. 那么对于给定的\(a_j=k\) ...
- 【比赛】HNOI2018 排列
这题原题... 这题题面七绕八绕,有点麻烦,反正最后转化就是一棵树,每个点有一个值,要把所有点选完,要求选择一个点必须是它的父亲和祖先已经全部被选了,贡献是这个点的权值乘上它被选择的排名 如果一个点是 ...
- [BZOJ5289][HNOI2018]排列(拓扑排序+pb_ds)
首先确定将所有a[i]向i连边之后会形成一张图,图上每条有向边i->j表示i要在j之前选. 图上的每个拓扑序都对应一种方案(如果有环显然无解),经过一系列推导可以发现贪心策略与合并的块的大小和w ...
随机推荐
- 扩展BSGS-传送门
很好的讲解:ZigZagK 好的讲解:mjtcn 某个模板:here 模板题: BSGS:ZigZagK的poj2417 exBSGS:ZigZagK的poj3243 — AC_Gibson 一般的板 ...
- LeetCode 24. Swap Nodes in Pairs (两两交换链表中的节点)
题目标签:Linked List 题目给了我们一组 linked list,让我们把每对nodes 互换位置. 新键一个dummy node,然后遍历list,每次建立 s1 和 s2 记录两个点,然 ...
- 中断控制及basepri 与 basepri_max
1.总开关 每个CPU有一个中断总开关.通过CPU中断控制寄存器实现.Cortex-M的中断控制寄存器包括:FAULTMASK.PRIMASK.BASEPRI.BASEPRI_MAX.总开关的本质是变 ...
- Rabbit MQ 客户端 API 开发
项目开始 第一步首先需要引入对应的 jar 包 <!-- https://mvnrepository.com/artifact/com.rabbitmq/amqp-client --> & ...
- 根据单个或多个字段对list对象去重
pojo 省略 在list 对象中,根据某一字段进行去重,重写Comparator /** * 去重 * * @param orderList * @return * @author ziggo * ...
- HIVE文件
注册表的本地实体文件, 察看位置,以及映射本地文件到注册表中的位置, HKEY_LOCAL_MACHINE\SYSTEM\CurrentControlSet\Control\hivelist 在这里写 ...
- JS函数 函数的作用,可以写一次代码,然后反复地重用这个代码。
什么是函数 函数的作用,可以写一次代码,然后反复地重用这个代码. 如:我们要完成多组数和的功能. var sum; sum = 3+2; alert(sum); sum=7+8 ; alert(sum ...
- PHP算法之无重复字符的最长子串
给定一个字符串,请你找出其中不含有重复字符的 最长子串 的长度. 示例 1: 输入: "abcabcbb"输出: 3 解释: 因为无重复字符的最长子串是 "abc&quo ...
- 阿里云 Aliplayer高级功能介绍(三):多字幕
基本介绍 国际化场景下面,播放器支持多字幕,可以有效解决视频的传播障碍难题,该功能适用于视频内容在全球范围内推广,阿里云的媒体处理服务提供接口可以生成多字幕,现在先看一下具体的效果: WebVTT格式 ...
- VS开发工具的常用插件
转 http://www.spersky.com/post/vsPlugins.html 我目前主要用的是Hide Main Page——公司配给的电脑屏幕分辨率好小,还是1366*768的,去掉头可 ...