JZOJ 5643. 【NOI2018模拟4.10】最小代价

最新推荐文章于 2019-06-15 15:05:49 发布

Felix-Lee

最新推荐文章于 2019-06-15 15:05:49 发布

阅读量314

点赞数 1

分类专栏：单调性 LCT 图论贪心文章标签： LCT 动态树图论贪心

本文链接：https://blog.csdn.net/liyizhixl/article/details/79911857

版权

图论同时被 3 个专栏收录

84 篇文章 0 订阅

订阅专栏

贪心

30 篇文章 1 订阅

订阅专栏

单调性

13 篇文章 0 订阅

订阅专栏

Description

给定一张n个点m条边的无向图，点编号1到n，每个点x有两个权值ax和bx。给定k，选出图中一个大小为k的点集S，使得S中任意两个点之间存在仅经过这个点集中的点的路径。S也存在两个权值aS和bS：aS为S中所有点的ax的最大值；bS为S中所有点的bx的最大值。求最小化aS+bS并输出这个最小值。

Input

第一行三个数 n,m,k。
接下来n行，第i+1行两个数ai和bi。
接下来m行，每行两个数x,y表示x和y之间有一条无向边。保证x和y不相等，任意两个点之间最多一条边。

Output

输出要求的最小值，若无解，输出“no solution”（引号不输出）。

Sample Input

5 7 3
1 2
2 1
2 2
1 1
1 2
1 3
4 5
2 3
5 1
4 2
1 4
3 4

Sample Output

Data Constraint

Subtask1：1 ≤ n ≤ 20, 1 ≤ m ≤ 100，时间限制1s，分值7。
Subtask2：1 ≤ n, m ≤ 5000，时间限制1s，分值17。
Subtask3：2 ≤ n ≤ 300000, m = n − 1，保证图连通，时间限制2s，分值32。
Subtask4：1 ≤ n ≤ 300000, 1 ≤ m ≤ 500000，时间限制3s，分值44。
保证权值均为正整数且不超过1,000,000,000。

Solution

看到这种有两个限制的题，先考虑固定一种限制，再想办法同时处理另一个限制。
由于是选择点集，权值又在点上不好处理，于是我们可以将权值转移到边上：
一条边的边权 $(a,b)$ 就设为两端点的 (a,b) max起来： $(max(a1,a2),max(b1,b2))$
那么选了一条边就相当于选了其两个端点，这样就方便多了。
先将 $a$ 和 $b$ 的值离散化，从小到大排序。
设 $mxa$ 和 $mxb$ ，那么答案即为 $ans=min(ans,mxa+mxb)$ 。
我们使 $mxa$ 从小到大的依次增加（开始时使 $mxb$ 为最大的 $b$ ），
每次把当前 $mxa$ 对应的边加入图中（用边集数组连边即可），
我们发现可以像维护最小生成树那样，加入小边，删去之前较大的边（LCT维护）（问题①）。
即：若该边的 $b$ 小于树中边最大的 $b$ ，就用该边替换之，否则就不加入。
若某个时刻加边后连通块大小 $\ge k$ 就可以更新答案了。
这时我们就将 $mxb$ 减少，并删去树中对应的边（没有就不删），
这样继续更新答案，知道删到某一个 $mxb$ 中的边后连通块大小 <script type="math/tex" id="MathJax-Element-5385"> $mxb$ 了，
之后就继续增加 $mxa$ ，重复以上操作。
现在解决 问题①：为什么这样就是正确的呢？
因为删去 $b$ 值较大的边，（在 $mxb$ 减少时它会先被删），就尽可能保持了它之后的连通性。
因为 $mxa$ 和 $mxb$ 都是单调的，于是时间复杂度就是 $O(N\ log\ N)$ 。

Code

#include<cstdio>
#include<algorithm>
#include<map>
#include<cctype>
using namespace std;
const int N=3e5+5,M=N*3,inf=2e9+1;
struct data
{
    int x,y,a,b;
}c[M];
int n,m,k,na,nb,tot1,tot2,cnta,cntb,top,ans=inf;
int firsta[M],nexa[M],ena[M];
int firstb[M],nexb[M],enb[M];
int a[N],b[N],aa[M],bb[M];
int fa[M],size[M],sz[M],s[M][2],num[M],mx[M],st[M];
bool rev[M],bz[M];
map<int,int>mpa,mpb;
inline int read()
{
    int X=0,w=0; char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
inline int max(int x,int y)
{
    return x>y?x:y;
}
inline int min(int x,int y)
{
    return x<y?x:y;
}
inline bool pd(int x)
{
    return x==s[fa[x]][1];
}
inline bool isroot(int x)
{
    return x^s[fa[x]][0] && x^s[fa[x]][1];
}
inline void reverse(int x)
{
    if(x) swap(s[x][0],s[x][1]),rev[x]^=1;
}
inline bool update(int x)
{
    size[x]=size[s[x][0]]+size[s[x][1]]+1+sz[x];
    num[x]=mx[num[s[x][0]]]>mx[num[s[x][1]]]?num[s[x][0]]:num[s[x][1]];
    if(mx[x]>mx[num[x]]) num[x]=x;
}
inline void down(int x)
{
    if(rev[x])
    {
        reverse(s[x][0]),reverse(s[x][1]);
        rev[x]=false;
    }
}
inline void rotate(int x)
{
    int y=fa[x],w=pd(x);
    if((fa[x]=fa[y]) && !isroot(y)) s[fa[y]][pd(y)]=x;
    if(s[y][w]=s[x][w^1]) fa[s[y][w]]=y;
    s[fa[y]=x][w^1]=y;
    update(y);
}
inline void splay(int x)
{
    for(int y=st[top=1]=x;!isroot(y);y=fa[y]) st[++top]=fa[y];
    while(top) down(st[top--]);
    for(int y;!isroot(x);rotate(x))
        if(!isroot(y=fa[x])) rotate(pd(x)==pd(y)?y:x);
    update(x);
}
inline void access(int x)
{
    for(int y=0;x;x=fa[y=x])
    {
        splay(x);
        sz[x]+=size[s[x][1]];
        sz[x]-=size[s[x][1]=y];
        update(x);
    }
}
inline void mkroot(int x)
{
    access(x),splay(x),reverse(x);
}
inline void link(int x,int y)
{
    mkroot(x),mkroot(y);
    sz[fa[x]=y]+=size[x];
    update(y);
}
inline void cut(int x,int y)
{
    mkroot(x),access(y),splay(y);
    s[y][0]=fa[x]=0;
}
inline bool connect(int x,int y)
{
    mkroot(x),access(y),splay(y);
    int z=y;while(s[z][0]) z=s[z][0];
    return x==z;
}
inline void inserta(int x,int y)
{
    nexa[++tot1]=firsta[x];
    firsta[x]=tot1;
    ena[tot1]=y;
}
inline void insertb(int x,int y)
{
    nexb[++tot2]=firstb[x];
    firstb[x]=tot2;
    enb[tot2]=y;
}
inline void insa(int x,int pos)
{
    if(!mpa[x]) mpa[x]=++cnta;
    inserta(mpa[x],pos);
}
inline void insb(int x,int pos)
{
    if(!mpb[x]) mpb[x]=++cntb;
    insertb(mpb[x],pos);
}
int main()
{
    n=read(),m=read(),k=read();
    for(int i=1;i<=n;i++) a[i]=read(),b[i]=read(),size[i]=1;
    if(k==1)
    {
        int sum=inf;
        for(int i=1;i<=n;i++) sum=min(sum,a[i]+b[i]);
        return 0&printf("%d",sum);
    }
    for(int i=1;i<=m;i++)
    {
        c[i].x=read(),c[i].y=read();
        c[i].a=max(a[c[i].x],a[c[i].y]);
        c[i].b=max(b[c[i].x],b[c[i].y]);
        insa(c[i].a,i);
        insb(c[i].b,i);
        aa[i]=c[i].a,bb[i]=c[i].b;
    }
    sort(aa+1,aa+1+m);
    sort(bb+1,bb+1+m);
    na=unique(aa+1,aa+1+m)-(aa+1);
    nb=unique(bb+1,bb+1+m)-(bb+1);
    int mxb=nb;
    for(int i=1;i<=na;i++)
    {
        for(int j=firsta[mpa[aa[i]]];j;j=nexa[j])
        {
            int nn=ena[j],x=c[nn].x,y=c[nn].y,z=nn+n;
            if(bz[nn]) continue;
            mx[num[z]=z]=c[nn].b;
            if(!connect(x,y)) link(x,z),link(z,y); else
            {
                mkroot(x),access(y),splay(y);
                if(c[nn].b<mx[num[y]])
                {
                    int pos=num[y];
                    cut(c[pos-n].x,pos);
                    cut(pos,c[pos-n].y);
                    link(x,z),link(z,y);
                }
            }
            mkroot(x),access(y),splay(y);
            while(size[y]>=k+k-1)
            {
                ans=min(ans,aa[i]+bb[mxb]);
                for(int l=firstb[mpb[bb[mxb]]];l;l=nexb[l])
                {
                    int nn1=enb[l],xx=c[nn1].x,yy=c[nn1].y;
                    bz[nn1]=true;
                    if(connect(xx,yy))
                    {
                        mkroot(xx),access(yy),splay(yy);
                        int pos=num[yy];
                        cut(c[pos-n].x,pos);
                        cut(pos,c[pos-n].y);
                    }
                }
                if(!--mxb)
                {
                    if(ans==inf) puts("no solution"); else printf("%d",ans);
                    return 0;
                }
                mkroot(x),access(y),splay(y);
            }
        }
    }
    if(ans==inf) puts("no solution"); else printf("%d",ans);
    return 0;
}

Felix-Lee

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
JZOJ 5643. 【NOI2018模拟4.10】最小代价

JZOJ 5643. 【NOI2018模拟4.10】最小代价 Description给定一张n个点m条边的无向图，点编号1到n，每个点x有两个权值ax和bx。给定k，选出图中一个大小为k的点集S，使得S中任意两个点之间存在仅经过这个点集中的点的路径。S也存在两个权值aS和bS：aS为S中所有点的ax的最大值；bS为S中所有点的bx的最大值。求最小化aS+bS并输出这个最小值。I
复制链接

扫一扫

专栏目录