Description:

给定 \(n\) 个整数 \(a_1, a_2, \dots, a_n, 0 \le a_i \le n\),以及 \(n\) 个整数 \(w_1, w_2, \dots, w_n\)。称 \(a_1, a_2, \dots, a_n\)的 一个排列 \(a_{p[1]}, a_{p[2]}, \dots, a_{p[n]}\)为 \(a_1, a_2, \dots, a_n\)的一个合法排列,当且仅当该排列满足:对于任意 的 \(k\) 和任意的 \(j\),如果 \(j \le k\),那么 \(a_{p[j]}\)不等于 \(p[k]\)。(换句话说就是:对于任意的 \(k\) 和任意的 \(j\),如果 \(p[k]\)等于 \(a_{ p[j] }\),那么 $k<j $ 。)定义这个合法排列的权值为 ​\(w_{p[1]} + 2w_{p[2]} + \dots + nw_{p[n]}\)。

Hint:

\(n \le 5*10^5\)

Solution:

估计很多人是因为看不懂题意才做不出来的。。。

。。。题面真的反人类

如果 \(p[k]\)等于 \(a_{ p[j] }\),那么 $ k<j $

换个说法:

如果 \(p[k]=a_{p[j]}\) 那么 \(k<j\)

即 \(y=a_x\) 且\(y\) 排在 \(x\)前面

又\(a_x=a_x\)

考虑拓扑排序,连边 \(a_x->x\)

这样每个点入度为1

连边后一定是一颗以0为根的树

为什么?因为\(a \in [1,n]\) ,而\(a_i \in [0,n]\)

若图中只有\(n\)个点,那么由于有\(n\)条边

一定成环

所有\(0\)必然为起点

如果有环,则无法满足限制说明合法排列不存在

接下来考虑在树上贪心

每次找到最小的权值\(w_i\),若此时\(fa[i]=0\),则可直接选\(i\)

否则需要把\(w_i​\)合并到\(fa[i]​\)上,等待\(fa[i]​\)选完再选

但有一个问题,这样操作后每个点变成了序列,如何比较优劣?

考虑 $ ab $ 和 $ ba $ 两种合并后的序列的答案(假设当前是第 $ i $ 位)

\[W_{ab}=\sum_{j=1}^{m_1}(i+j)w_{a_j}+\sum_{j=1}^{m_2}(i+j+m_1)w_{b_j}​$$$$W_{ba}=\sum_{j=1}^{m_2}(i+j)w_{b_j}+\sum_{j=1}^{m_1}(i+j+m_2)w_{a_j}​
\]

\[W_{ab}-W_{ba}=m_1W_b-m_2W_a​
\]

如果\(W_{ab}\gt W_{ba}\Rightarrow \frac{W_a}{m_1}\lt\frac{W_b}{m_2}​\)也就是平均权值小的放前面答案会更优

这里使用 稍加修改的堆+并查集 实现

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mxn=5e5+5;
int n,cnt,tot,f[mxn],s[mxn],hd[mxn],vis[mxn],anc[mxn];
ll ans,w[mxn]; struct ed {
int to,nxt;
}t[mxn<<1]; struct data {
int id,sz;
ll val;
friend bool operator < (data a,data b) {
return a.val*b.sz>b.val*a.sz; //按平均值排序
}
}tp; priority_queue<data > q; inline void add(int u,int v) {
t[++cnt]=(ed) {v,hd[u]}; hd[u]=cnt;
} int find(int x) {
return anc[x]==x?x:anc[x]=find(anc[x]);
} void dfs(int u,int fa) {
vis[u]=1; ++tot;
for(int i=hd[u];i;i=t[i].nxt) {
int v=t[i].to;
if(v==fa) continue ;
if(vis[v]) {
puts("-1");
exit(0);
}
dfs(v,u);
}
} int main()
{
scanf("%d",&n); s[0]=1; //别忘记赋初值
for(int i=1;i<=n;++i) {
scanf("%d",&f[i]);
add(f[i],i); anc[i]=i,s[i]=1;
}
dfs(0,-1); if(tot<n) {puts("-1"); return 0;}
for(int i=1;i<=n;++i) {
scanf("%lld",&w[i]);
q.push((data){i,1,w[i]});
}
while(!q.empty()) {
tp=q.top(); q.pop(); int u=find(tp.id);
if(s[u]!=tp.sz) continue ;/*已经有合并后新的节点,丢弃以前的*/ anc[u]=find(f[u]);
ans+=w[u]*s[anc[u]],s[anc[u]]+=s[u],w[anc[u]]+=w[u]; //每次先计算贡献,再合并
if(anc[u]) q.push((data){anc[u],s[anc[u]],w[anc[u]]});
}
printf("%lld",ans);
return 0;
}

[HNOI2018]排列的更多相关文章

  1. 【BZOJ5289】[HNOI2018]排列(贪心)

    [BZOJ5289][HNOI2018]排列(贪心) 题面 BZOJ 洛谷 题解 这个限制看起来不知道在干什么,其实就是找到所有排列\(p\)中,\(p_k=x\),那么\(k<j\),其中\( ...

  2. 5289: [Hnoi2018]排列

    5289: [Hnoi2018]排列 链接 分析: 首先将题意转化一下:每个点向a[i]连一条边,构成了一个以0为根节点的树,要求选一个拓扑序,点x是拓扑序中的第i个,那么价值是i*w[x].让价值最 ...

  3. bzoj 5289: [Hnoi2018]排列

    Description Solution 首先注意到实际上约束关系构成了一棵树 考虑这个排列 \(p\),编号为 \(a[i]\) 的出现了,\(i\) 才可以出现 那么如果连边 \((a[i],i) ...

  4. [HNOI2018]排列[堆]

    题意 给定一棵树,每个点有点权,第 \(i\) 个点被删除的代价为 \(w_{p[i]}\times i\) ,问最小代价是多少. 分析 与国王游戏一题类似. 容易发现权值最小的点在其父亲选择后就会立 ...

  5. BZOJ5289: [Hnoi2018]排列

    传送门 第一步转化,令 \(q[p[i]]=i\),那么题目变成: 有一些 \(q[a[i]]<q[i]\) 的限制,\(q\) 必须为排列,求 \(max(\sum_{i=1}^{n}w[i] ...

  6. loj2509 hnoi2018排列

    题意:对于a数组,求它的一个合法排列的最大权值.合法排列:对于任意j,k,如果a[p[j]]=p[k],那么k<j. 权值:sigma(a[p[i]]*i).n<=50W. 标程: #in ...

  7. BZOJ.5289.[AHOI/HNOI2018]排列(贪心 heap)

    BZOJ LOJ 洛谷 \(Kelin\)写的挺清楚的... 要求如果\(a_{p_j}=p_k\),\(k\lt j\),可以理解为\(k\)要在\(j\)之前选. 那么对于给定的\(a_j=k\) ...

  8. 【比赛】HNOI2018 排列

    这题原题... 这题题面七绕八绕,有点麻烦,反正最后转化就是一棵树,每个点有一个值,要把所有点选完,要求选择一个点必须是它的父亲和祖先已经全部被选了,贡献是这个点的权值乘上它被选择的排名 如果一个点是 ...

  9. [BZOJ5289][HNOI2018]排列(拓扑排序+pb_ds)

    首先确定将所有a[i]向i连边之后会形成一张图,图上每条有向边i->j表示i要在j之前选. 图上的每个拓扑序都对应一种方案(如果有环显然无解),经过一系列推导可以发现贪心策略与合并的块的大小和w ...

随机推荐

  1. .NET编码解码(HtmlEncode与HtmlDecode)

    编码代码: System.Web.HttpUtility.HtmlEncode("<a href=\"http://hovertree.com/\">何问起& ...

  2. vs问题解决:an operation is not legal in the current state

    debug时弹出提示框:内容有:an operation is not legal in the current state 解决方案: Go to Tools > Options > D ...

  3. C#的值传递与引用传递

    值传递:在使用值传递时,是把变量的值传给函数,函数中对此变量的任何修改都不影响该变量本身的值. 引用传递:使用引用传递时,在函数中对此变量的修改会影响变量的值. 说简单点,值传递,就是我把身份证复印件 ...

  4. Python_网络编程udp-飞秋自动攻击

    # 模拟一个接收数据import socketimport time def auto_hack(udp_socket, recv_msg, revc_ip, revc_port=2425): # 发 ...

  5. # 20165206 2017-2018-2 《Java程序设计》第4周学习总结

    20165206 2017-2018-2 <Java程序设计>第4周学习总结 教材学习内容总结 继承:继承是一种由已有的类创建新类的机制. 子类和父类:由继承得到的类称为子类,被继承的类称 ...

  6. tensorflow:验证码的识别(下)

    上两篇详细的说明了验证码的识别,不过我们采用的是方法二,下面采用方法一.注意和方法二的区别. 验证码识别方法一: 把标签转为向量,向量长度为40.(4位数字验证码) 验证码的生成和tf.record的 ...

  7. 使用python解决算法和数据结构--使用栈实现进制转换

    可以将10进制数据转换成2进制,8进制,16进制等. 晚上练练算法和数据结构哈. # coding = utf-8 class Stack: def __init__(self): self.item ...

  8. Node 杂技

    1.关于require 当文件夹a中含有index.js时,在b.js中如果有require("文件夹a的路径"),则将会自动执行index.js的语句

  9. kickstart-E

    A题 简答模拟题 #include <iostream> #include<stdio.h> #include <set> #include <algorit ...

  10. bzoj 5099: [POI2018]Pionek

    题解: 还是比较简单的一道题 考虑现在有一个向量,当且仅当下一个向量与它夹角<90度这个向量的模长才会增加 接下来怎么做呢 如果我们去枚举初始向量,向量方向会随着新增向量而变化 随着不断顺时针的 ...