• 【ARC069F】Flags(2-SAT,Tarjan,线段树优化建图)


    注意:本题的点数可以相比题解优化一半。


    首先先二分答案。

    然后判断能否使得两两旗子之间的距离都大于 m i d mid mid

    然后发现这是一个 2-SAT 问题。

    2-SAT 问题:通俗地说,有 n n n 个 bool 变量 a i a_i ai,并给出一些形如 a i ⊕ a j = 0 / 1 a_ioplus a_j=0/1 aiaj=0/1 的条件(其中 ⊕ oplus 可以是 and ⁡ operatorname{and} and or ⁡ operatorname{or} or xor ⁡ operatorname{xor} xor),然后询问满足这组方程的一组解 a i a_i ai

    这题同样也给出了一些条件,比如说对于所有的 i i i x i x_i xi y i y_i yi 中必须选一个但是不能同时选。同时在 m i d mid mid 的约束下,也有一些数不能同时选。

    那我们可以设 bool 变量 a i a_i ai 表示取不取 x i x_i xi,bool 变量 b i b_i bi 表示取不取 y i y_i yi。那么 a i a_i ai b i b_i bi 不能同时等于 true ext{true} true,且两者中肯定有一者为 true ext{true} true,这就可以用 a i xor ⁡ b i = 1 a_ioperatorname{xor}b_i=1 aixorbi=1 来表示。

    那么这就变成了一个 2-SAT 问题了。

    但是这道题只用判断是否有解,也就是可行性。这时可以用 Tarjan 求解:

    首先把每一个 bool 变量 x x x 拆成两个命题:命题 x 0 x_0 x0 表示 x = false x= ext{false} x=false,命题 x 1 x_1 x1 表示 x = true x= ext{true} x=true。设命题 y y y 的反面是 y ′ y' y,显然 x 0 ′ x_0' x0 就是 x 1 x_1 x1

    那么显然当 x 0 x_0 x0 成立时, x 1 x_1 x1 不成立;当 x 1 x_1 x1 成立时, x 0 x_0 x0 不成立。而且 x 0 x_0 x0 x 1 x_1 x1 之间肯定有一者成立。

    那么 a i xor ⁡ b i = 1 a_ioperatorname{xor}b_i=1 aixorbi=1 就等价于:

    a i , 0 a_{i,0} ai,0 成立时, b i , 1 b_{i,1} bi,1 成立;当 b i , 1 b_{i,1} bi,1 成立时, a i , 0 a_{i,0} ai,0 成立。

    a i , 1 a_{i,1} ai,1 成立时, b i , 0 b_{i,0} bi,0 成立;当 b i , 0 b_{i,0} bi,0 成立时, a i , 1 a_{i,1} ai,1 成立。

    考虑用图论的方式来表达这种关系。

    设单向边 ( u , v ) (u,v) (u,v) 表示当命题 u u u 成立时,命题 v v v 也必定成立。

    那么上述关系就可以表示成单向边 ( a i , 0 , b i , 1 ) (a_{i,0},b_{i,1}) (ai,0,bi,1) ( b i , 1 , a i , 0 ) (b_{i,1},a_{i,0}) (bi,1,ai,0) ( a i , 1 , b i , 0 ) (a_{i,1},b_{i,0}) (ai,1,bi,0) ( b i , 0 , a i , 1 ) (b_{i,0},a_{i,1}) (bi,0,ai,1)

    然后题目中还有一些条件:某两个数不能同时取,即某两个 bool 变量 a a a b b b 不能同时为 1 1 1,即 a and ⁡ b = 0 aoperatorname{and} b=0 aandb=0,考虑也用图论来表示这个。

    发现可以用有向边 ( a 1 , b 0 ) (a_1,b_0) (a1,b0) ( b 1 , a 0 ) (b_1,a_0) (b1,a0) 来表示。

    那么我们就能把所有的条件都用图来表示了。

    至于如何判断解的可行性:

    我们可以先对这个图做一遍 Tarjan,然后看是否存在两个反命题在一个环中(即出现 “当命题 x x x 成立时,可得命题 x ′ x' x 成立,当 x ′ x' x 成立时,也可得 x x x 成立,但 x x x x ′ x' x 为相反的命题” 这种情况)。如果存在,那么原方程无解,否则有解。

    然后这道题需要用线段树优化建图才能过。

    代码如下:

    //1~n 取x[i]
    //n+1~2n 取y[i]
    //2n+1~3n 不取x[i]
    //3n=1~4n 不取y[i] 
    #include<bits/stdc++.h>
    
    #define N 20010
    
    using namespace std;
    
    struct data
    {
    	int val,id;
    	data(){};
    	data(int a,int b){val=a,id=b;} 
    }a[N<<1];
    
    int n;
    int node,id[N<<3];
    int idx,dfn[N*12],low[N*12];
    int top,sta[N*12];
    int tot,num[N*12];
    int cnt,head[N*12],to[N*44],nxt[N*44];
    bool ins[N*12];
    
    bool cmp(data a,data b)
    {
    	return a.val<b.val;
    }
    
    void adde(int u,int v)
    {
    	to[++cnt]=v;
    	nxt[cnt]=head[u];
    	head[u]=cnt;
    }
    
    void build(int k,int l,int r)
    {
    	if(l==r)
    	{
    		id[k]=(n<<1)+a[l].id;
    		return;
    	}
    	id[k]=++node;
    	int mid=(l+r)>>1;
    	build(k<<1,l,mid);
    	build(k<<1|1,mid+1,r);
    	adde(id[k],id[k<<1]);
    	adde(id[k],id[k<<1|1]);
    }
    
    void link(int k,int l,int r,int ql,int qr,int rt)
    {
    	if(ql<=l&&r<=qr)
    	{
    		adde(rt,id[k]);
    		return;
    	}
    	int mid=(l+r)>>1;
    	if(ql<=mid) link(k<<1,l,mid,ql,qr,rt);
    	if(qr>mid) link(k<<1|1,mid+1,r,ql,qr,rt);
    }
    
    void Tarjan(int u)
    {
    	dfn[u]=low[u]=++idx;
    	sta[++top]=u;
    	ins[u]=true;
    	for(int i=head[u];i;i=nxt[i])
    	{
    		int v=to[i];
    		if(!dfn[v])
    		{
    			Tarjan(v);
    			low[u]=min(low[u],low[v]);
    		}
    		else if(!num[v])
    			low[u]=min(low[u],dfn[v]);
    	}
    	if(low[u]==dfn[u])
    	{
    		tot++;
    		int v;
    		do
    		{
    			v=sta[top];
    			top--;
    			ins[v]=false;
    			num[v]=tot;
    		}while(u!=v);
    	}
    }
    
    bool check(int mid)
    {
    	cnt=idx=tot=0,node=n<<2;
    	memset(head,0,sizeof(head));
    	memset(dfn,0,sizeof(dfn));
    	memset(low,0,sizeof(low));
    	memset(num,0,sizeof(num));
    	build(1,1,n<<1);
    	for(int i=1;i<=n;i++)
    	{
    		adde(i,n+(n<<1)+i);//取x[i]则不取y[i] 
    		adde(n+(n<<1)+i,i);//不取y[i]则取x[i] 
    		adde(n+i,(n<<1)+i);//取y[i]则不取x[i]
    		adde((n<<1)+i,n+i); //不取x[i]则取y[i]
    	}
    	int nowl=1,nowr=1;
    	for(int i=1;i<=(n<<1);i++)
    	{
    		while(nowl<i&&a[i].val-a[nowl].val>=mid) nowl++;
    		while(nowr<(n<<1)&&a[nowr+1].val-a[i].val<mid) nowr++;
    		if(nowl<i) link(1,1,n<<1,nowl,i-1,a[i].id);
    		if(nowr>i) link(1,1,n<<1,i+1,nowr,a[i].id);
    	}
    	for(int i=1;i<=(n<<1);i++)
    		if(!dfn[i])
    			Tarjan(i);
    	for(int i=1;i<=(n<<1);i++)
    		if(num[i]==num[i+(n<<1)])
    			return false;
    	return true;
    }
    
    int main()
    {
    	scanf("%d",&n);
    	for(int i=1;i<=n;i++)
    	{
    		scanf("%d%d",&a[i].val,&a[n+i].val);
    		a[i].id=i,a[n+i].id=n+i;
    	}
    	sort(a+1,a+(n<<1)+1,cmp);
    	int l=0,r=1e9,ans;
    	while(l<=r)
    	{
    		int mid=(l+r)>>1;
    		if(check(mid)) ans=mid,l=mid+1;
    		else r=mid-1;
    	}
    	printf("%d
    ",ans);
    	return 0;
    }
    
  • 相关阅读:
    xgboost中XGBClassifier()参数详解
    xgboost使用经验总结
    特征选择之Chi卡方检验
    ldd 查看程序依赖库
    paddle——docker实践
    paddle实践
    java常用设计模式
    Java内存模型
    什么是线程?
    MySQL数据库提供了那些函数?
  • 原文地址:https://www.cnblogs.com/ez-lcw/p/14448648.html
Copyright © 2020-2023  润新知