题目链接:E-Xor Tree
题目大意: 给出一棵树,每个节点有一个权值,现在给出每个点权值的范围 [ l i , r i ] [l_i,r_i] [li,ri]和相邻两点权值异或的值,问全部节点所有可能的取值所构成的方案数。
题目分析: 考虑如果某一个节点的值已知,那么剩下所有点的值都可以求出来。假设
1
1
1号点的值为
0
0
0,求出其他点的值,然后发现,如果令第一个点的值
w
1
=
a
w_1=a
w1=a那么其他点的值都会
⊕
a
\oplus{a}
⊕a,因此问题变为了求合法的
a
a
a的数量,限制有
n
n
n个不等式,形式为
l
i
≤
w
i
⊕
a
≤
r
i
l_i \leq w_i \oplus a\leq r_i
li≤wi⊕a≤ri也就是
a
∈
∩
i
=
1
n
S
i
a \in \cap_{i=1}^{n} S_i
a∈∩i=1nSi
S
i
=
[
l
i
,
r
i
]
⊕
w
i
S_i=[l_i,r_i] \oplus w_i
Si=[li,ri]⊕wi因为如果
a
′
∈
S
a'\in S
a′∈S那么
a
′
⊕
w
i
∈
[
l
i
,
r
i
]
a'\oplus w_i \in [l_i,r_i]
a′⊕wi∈[li,ri],所以现在的问题是如何对每个点所在的区间求出
S
i
S_i
Si。我们考虑如果将区间
[
l
i
,
r
i
]
[l_i,r_i]
[li,ri]拆分成
l
o
g
2
w
i
log_2w_i
log2wi个区间,保证每个区间的左端点二进制拆分为
x
.
.
.
x
x
0...00
x...xx0...00
x...xx0...00而右端点二进制拆分为
x
.
.
x
x
1...11
x..xx1...11
x..xx1...11那么在这个区间中,异或
w
i
w_i
wi所形成的区间的值一定也是连续的,其中可以通过分类讨论求出异或后所得区间的端点。
具体可以想线段树一样,第一层表示的是
[
0
,
2
30
−
1
]
[0,2^{30}-1]
[0,230−1]依次类推,由于这棵线段树并不需要实际操作,因此没有空间复杂度,只需要每次以
O
(
l
o
g
N
)
O(logN)
O(logN)的代价求出能够被处理的区间。
这样再去维护端点所构成的若干区间的并集然后加入可行解,最后把
n
n
n个集合求交集,并求出方案数。
完整代码如下。
题目代码:
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<utility>
#define mxn 100005
using namespace std;
vector<pair<int,int> >Ns;
int EHead[mxn],ENext[mxn<<1],ETo[mxn<<1],ECost[mxn<<1],Lt;
int x[mxn],y[mxn],w[mxn],cnt;
void AddEdge(int fr,int to,int cc){
ETo[Lt]=to;
ECost[Lt]=cc;
ENext[Lt]=EHead[fr];
EHead[fr]=Lt++;
}
void DFS(int s,int fa){
for(int i=EHead[s];i+1;i=ENext[i])
if(ETo[i]^fa){
w[ETo[i]]=w[s]^ECost[i];
DFS(ETo[i],s);
}
}
void pushup(int L,int R,int d,int v){
int _d=0,_v=v,l,r,_L=L;
while(_v)++_d,_v>>=1;
if(_d<d)l=L,r=l+(1<<d)-1;
if(_d==d)l=(L>>d)<<d,r=L+(1<<d)-1;
if(_d>d)l=L^v,l=(l>>d)<<d,r=l+(1<<d)-1;
//printf("%d %d\n",l,r);
Ns.push_back(make_pair(l,r));
}
void work(int L,int R,int l,int r,int d,int v){
if(l<=L&&r>=R){
//printf("%d %d ",L,R);
pushup(L,R,d,v);
return;
}
int mid=(L+R)>>1;
if(r<=mid)work(L,mid,l,r,d-1,v);
else if(l>mid)work(mid+1,R,l,r,d-1,v);
else work(L,mid,l,mid,d-1,v),work(mid+1,R,mid+1,r,d-1,v);
}
int main()
{
int n;
memset(EHead,-1,sizeof(EHead));
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d %d",&x[i],&y[i]);
for(int i=1;i<=n-1;i++){
int fr,to,cc;
scanf("%d %d %d",&fr,&to,&cc);
AddEdge(fr,to,cc);
AddEdge(to,fr,cc);
}
DFS(1,0);
for(int i=1;i<=n;i++){
work(0,(1<<30)-1,x[i],y[i],30,w[i]);
// work(0,(1<<30)-1,5,20,30,6);
}
vector<pair<int,int> >res;
for(auto x:Ns){
res.push_back(make_pair(x.first,1));
res.push_back(make_pair(x.second+1,-1));
}
sort(res.begin(),res.end());
int dep=0,ans=0;
for(int i=0;i<=res.size()-1;i++){
dep+=res[i].second;
if(dep==n)if(i+1<res.size())
ans+=res[i+1].first-res[i].first;
}
printf("%d\n",ans);
return 0;
}