标签：tmp Antihype int tt Kuroni CF1305G fi id se
一、题目
二、解法
你感觉这道题有点像生成树模型，但是因为边有方向所以麻烦。
可以巧妙地转化成无向生成树模型，我们把 \((i,j)\) 之间边的权值设置成 \(a_i+a_j\)，那么如果是 \(i\) 把 \(j\) 拉进连通块，我们多算了 \(a_j\)，如果是 \(j\) 把 \(i\) 拉进连通块我们多算了 \(a_i\)，所以每个点的点权会被多算一次，我们增加 \(a_{n+1}=0\) 连向所有点来表示直接进入组织的情况，那么用最大生成树减去点权和就是答案。
那么问题转化成了求最大生成树，可以考虑 \(\tt kruskal\) 算法，也就是从大到小枚举边权 \(w\)，接下来就看哪些边有这种权值，其实就是把 \(w\) 的二进制拆分成两个不相交的子集，那么就可以套子集枚举了，用并查集暴力连起来即可，时间复杂度 \(O(3^{18}\cdot \alpha)\)
还有一种更好的方法，考虑 \(\tt Borůvka\) 算法，大概内容是对于当前的每一个连通块寻找权值最大的一个出边，这一轮把找到的边都连起来，因为每次连通块个数减半，所以一共只会有 \(O(\log n)\) 轮，复杂度关键在于找边。
二进制的东西不好用数据结构维护的，干脆每一轮我们现场算。注意到边权是 \(a_i+a_j\)，我们先固定 \(a_i\)，然后对于 \(i\) 关于全集的补集 \(s\)，我们找到其中最大的 \(a_j\)，考虑它们的连边，但是要注意 \(a_i\) 和 \(a_j\) 需要在不同的连通块内。
这相当于一个子集内求最大值的问题，可以用 \(\tt fwt\) 正变换的方法来做，也就是我们从小到大枚举数位 \(i\)，然后对于包含数位 \(i\) 的集合 \(s\) 考虑 \(s-2^i\) 来更新 \(s\)，这本质上是一个分治统计的过程，时间复杂度 \(O(2^{18}\cdot 18)\)
因为不能在同一连通块内所以我们要算的是不在同一联通块内的最大值和次大值，时间复杂度 \(O(2^{18}\cdot18\cdot\log n)\)
最后小结一下吧，为什么 \(\tt Borůvka\) 算法能做到更优的复杂度？我觉得是这道题的边权是和点权高度相关的，而且我们不是很清楚具体的边，这样 \(\tt kruskal\) 算法可能不优，而 \(\tt prim\) 算法是单源扩展所以可能也不行，\(\tt Borůvka\) 算法是和点高度相关的。
#include <cstdio>
#include <iostream>
using namespace std;
#define pii pair<int,int>
#define make make_pair
#define fi first
#define se second
const int M = 300005;
int read()
{
int x=0,f=1;char c;
while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
return x*f;
}
int n,m,cnt,a[M],fa[M];pii tmp[M],b[M][2];long long ans;
int find(int x)
{
if(x!=fa[x]) fa[x]=find(fa[x]);
return fa[x];
}
void merge(int u,int v)
{
int x=find(u),y=find(v);fa[x]=y;
}
void solve()
{
for(int i=0;i<m;i++)
b[i][0]=b[i][1]=tmp[i]=make(-1,-1);
for(int i=1;i<=n;i++)
{
int id=find(i);
if(b[a[i]][0].se==-1)
b[a[i]][0]=make(a[i],id);
else if(b[a[i]][1].se==-1 && b[a[i]][0].se!=id)
b[a[i]][1]=make(a[i],id);
}
for(int i=0;i<18;i++)
for(int s=0;s<m;s++) if((s>>i)&1)
{
int t=s^(1<<i);
for(int j=0;j<2;j++)
{
if(b[t][j].fi>b[s][0].fi)
{
if(b[s][0].se!=b[t][j].se) b[s][1]=b[s][0];
b[s][0]=b[t][j];
}
else if(b[t][j].fi>b[s][1].fi && b[s][0].se!=b[t][j].se)
b[s][1]=b[t][j];
}
}
for(int i=1;i<=n;i++)
{
int s=(m-1)^a[i],id=find(i);
if(b[s][0].fi!=-1 && b[s][0].fi+a[i]>tmp[id].fi && b[s][0].se!=id)
tmp[id]=b[s][0],tmp[id].fi+=a[i];
if(b[s][1].fi!=-1 && b[s][1].fi+a[i]>tmp[id].fi && b[s][1].se!=id)
tmp[id]=b[s][1],tmp[id].fi+=a[i];
}
for(int i=1;i<=n;i++)
if(tmp[i].fi!=-1 && find(i)!=find(tmp[i].se))
merge(i,tmp[i].se),ans+=tmp[i].fi,cnt--;
}
signed main()
{
cnt=n=read()+1;m=1<<18;
for(int i=1;i<=n;i++)
fa[i]=i;
for(int i=1;i<n;i++)
a[i]=read(),ans-=a[i];
while(cnt>1) solve();
printf("%lld\n",ans);
}
标签：tmp,Antihype,int,tt,Kuroni,CF1305G,fi,id,se 来源： https://www.cnblogs.com/C202044zxy/p/14878009.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。