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
3
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 ,就用该边替换之,否则就不加入。若某个时刻加边后连通块大小 ≥k 就可以更新答案了。
这时我们就将 mxb 减少,并删去树中对应的边(没有就不删),
这样继续更新答案,知道删到某一个 mxb 中的边后连通块大小 <k <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;
}