Codeforces Round #507 (Div. 2, based on Olympiad of Metropolises) E. Network Safety(异或运算+SCC)

E. Network Safety

题意：

输入$n(5e5),m(min(\frac{n\times (n-1)}{2},5e5)),k(60)$
接下来一行输入$a_1,a_2,\dots ,a_n([0,2^k))$
接下来$m$行输入$u,v$,表示边，保证$a_u\ne a_v$。
问有多少个$(A,x)$,$A$是此图的点集(可以为空)，$x([0,2^k))$,使$x$和$A$中每个点点值异或后，每条边两端的点值依旧不同。

题解：

首先考虑对于每个值$x$分析，求有多少个子集符合要求。
$a_u\ne a_v所以a_u⨁x\ne a_v⨁x$
于是便可以分两种情况：

• $a_u⨁x\ne a_v$时,$a_v$有两种选择。
• $a_u⨁x=a_v$时，$a_v$只有一种选择。

因此可以把所有的$a_u⨁a_v=x$的$u,v$都拿出来，形成一种子图，求联通分量的个数$g$,把$2^g$加入答案即可。

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int mod=1e9+7;
const int N=5e5+9;
ll n,m,k,a[N];
map<ll,int>mp;
int cnt;
#define P pair<int,int>
P b[N];
vector<P>g[N];
ll ans;
int poww(int a,int b,int c){
    int ans=1,base=a;
    while(b){
        if(b&1)ans=1ll*ans*base%c;
        base=1ll*base*base%c;
        b>>=1;
    }
    return ans;
}
int solve1(){
    return 1ll*poww(2,n,mod)*(((1ll<<k)-cnt)%mod)%mod;
}
int f[N],siz[N];
int findf(int x){return f[x]==x?x:f[x]=findf(f[x]);}
int solve2(int x){
    set<int>s;
    for(auto i:g[x]){
        int fx=findf(i.first),fy=findf(i.second);
        if(fx!=fy){
            f[fx]=fy,siz[fy]+=siz[fx];
        }
    }
    int sum=0,sz=0;
    for(auto i:g[x]){
        if(f[i.first]==i.first&&!s.count(i.first))sum++,sz+=siz[i.first],s.insert(i.first);
        if(f[i.second]==i.second&&!s.count(i.second))sum++,sz+=siz[i.second],s.insert(i.second);
    }
    s.clear();
    for(auto i:g[x]){
        f[i.first]=i.first,f[i.second]=i.second;
        siz[i.first]=siz[i.second]=1;
    }
    return poww(2,n-sz+sum,mod);
}
int main(){
   // freopen("tt.in","r",stdin),freopen("tt.out","w",stdout);
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>n>>m>>k;
    for(int i=1;i<=n;i++)cin>>a[i];
    for(int i=1;i<=m;i++)cin>>b[i].first>>b[i].second;
    for(int i=1;i<=m;i++){
        ll t=a[b[i].first]^a[b[i].second];
        if(!mp.count(t))mp[t]=++cnt;
        g[mp[t]].push_back(b[i]);
    }
    ans=solve1();
    for(int i=1;i<=n;i++)f[i]=i,siz[i]=1;
    for(int i=1;i<=cnt;i++)ans=(ans+solve2(i))%mod;
    cout<<ans<<endl;
    return 0;
}