#988297 (C++23) No.2822 Lights Up! (Tree Edition)

• 問題: No.2822 Lights Up! (Tree Edition)
• ユーザー: FplusFplusF
• 提出日時: 2024-06-14 06:23:43
• 言語: C++23 (gcc 13.3.0 + boost 1.87.0)
• 結果: WA
• 実行時間: -
• コード長: 7,140 bytes
• コンパイル時間: 4,039 ms
• コンパイル使用メモリ: 281,704 KB
• 実行使用メモリ: 65,172 KB
• 最終ジャッジ日時: 2024-06-18 19:45:51
• 合計ジャッジ時間: 10,778 ms
• ジャッジサーバーID （参考情報）: judge2 / judge3

ファイルパターン	結果
sample	AC * 1
other	AC * 105 WA * 37

ソースコード

#include <bits/stdc++.h>
using namespace std;
using ll=long long;
using ull=unsigned long long;
using pll=pair<ll,ll>;
using tll=tuple<ll,ll,ll>;
using ld=long double;
const ll INF=(1ll<<60);
#define rep(i,n) for (ll i=0;i<(ll)(n);i++)
#define replr(i,l,r) for (ll i=(ll)(l);i<(ll)(r);i++)
#define all(v) v.begin(),v.end()
#define len(v) ((ll)v.size())
template<class T> inline bool chmin(T &a,T b){
    if(a>b){
        a=b;
        return true;
    }
    return false;
}
template<class T> inline bool chmax(T &a,T b){
    if(a<b){
        a=b;
        return true;
    }
    return false;
}
template<class T> struct fenwick_tree{
    vector<T> v;
    int n;
    fenwick_tree(int x){
        n=x+1;
        v.assign(n+1,0);
    }
    fenwick_tree(vector<T> &a){
        n=(int)a.size()+1;
        v.assign(n+1,0);
        if((int)a.size()==0) return;
        v[1]=a[0];
        for(int i=0;i<n-2;i++) v[i+2]=a[i+1]-a[i];
        for(int idx=1;idx<=n;idx++){
            if(n<idx+(idx&(-idx))) continue;
            v[idx+(idx&(-idx))]+=v[idx];
        }
    }
    void add_sub(int i,T x){
        i++;
        for(int idx=i;idx<=n;idx+=(idx&(-idx))){
            v[idx]+=x;
        }
    }
    void add(int l,int r,T x){  //[l,r)
        assert(0<=l&&l<=n);
        assert(0<=r&&r<=n);
        assert(l<=r);
        add_sub(l,x);
        add_sub(r,-x);
    }
    T get(int i){
        assert(0<=i&&i<n);
        i++;
        T ret=0;
        for(int idx=i;0<idx;idx-=(idx&(-idx))){
            ret+=v[idx];
        }
        return ret;
    }
    int lower_bound(T w){
        if(w<=0) return 0;
        int x=0;
        static int r=1;
        if(r==1){
            while((r<<1)<=n) r<<=1;
        }
        while((r<<1)<=n) r<<=1;
        for(int k=r;0<k;k>>=1){
            if(x+k<=n&&v[x+k]<w){
                w-=v[x+k];
                x+=k;
            }
        }
        x++;
        return x;
    }
};
struct tree{
    int n,max_k=1;
    vector<vector<int>> g;
    vector<vector<int>> parent;
    vector<int> depth;
    tree(int g_size){
        n=g_size;
        g.resize(n);
    }
    void add_edge(int u,int v){
        g[u].push_back(v);
        g[v].push_back(u);
    }
    void LCA(int root=0){
        while((1<<max_k)<n) max_k++;
        parent.assign(max_k,vector<int>(n,-1));
        depth.assign(n,-1);
        dfs(root,-1,0);
        for(int k=0;k+1<max_k;k++){
            for(int i=0;i<n;i++){
                if(parent[k][i]==-1){
                    parent[k+1][i]=-1;
                }else{
                    parent[k+1][i]=parent[k][parent[k][i]];
                }
            }
        }
    }
    void dfs(int v,int par,int d){
        parent[0][v]=par;
        depth[v]=d;
        for(auto &i:g[v]){
            if(i!=par){
                dfs(i,v,d+1);
            }
        }
    }
    int doubling(int x,int diff){
        for(int k=0;k<max_k;k++){
            if(diff&(1<<k)){
                x=parent[k][x];
            }
        }
        return x;
    }
    int get_lca(int u,int v){
        if(depth[u]<depth[v]) swap(u,v);
        int diff=depth[u]-depth[v];
        for(int k=0;k<max_k;k++){
            if(diff&(1<<k)){
                u=parent[k][u];
            }
        }
        if(u==v) return u;
        for(int k=max_k-1;0<=k;k--){
            if(parent[k][u]!=parent[k][v]){
                u=parent[k][u];
                v=parent[k][v];
            }
        }
        return parent[0][u];
    }
    int get_dist(int u,int v){
        return depth[u]+depth[v]-2*depth[get_lca(u,v)];
    }
    int on_path(int x,int u,int v){
        return get_dist(u,x)+get_dist(x,v)==get_dist(u,v);
    }
};
struct HLD{
    int n;
    vector<vector<int>> g;
    vector<int> par,sz,in,out,head;
    int in_ord=0,out_ord=0;
    HLD(int x){
        n=x;
        g.resize(n);
        par.resize(n);
        sz.assign(n,0);
        in.resize(n);
        out.resize(n);
        head.resize(n);
    }
    void add_edge(int u,int v){
        g[u].push_back(v);
        g[v].push_back(u);
    }
    void dfs_sz(int x,int p){
        sz[x]++;
        if(0<g[x].size()&&g[x].front()==p) swap(g[x].front(),g[x].back());
        for(auto &i:g[x]){
            if(i==p) continue;
            dfs_sz(i,x);
            sz[x]+=sz[i];
            if(sz[g[x].front()]<sz[i]) swap(g[x].front(),i);
        }
    }
    void dfs_hld(int x,int p){
        par[x]=p;
        in[x]=in_ord;
        in_ord++;
        for(auto &i:g[x]){
            if(i==p) continue;
            if(i==g[x].front()) head[i]=head[x];
            else head[i]=i;
            dfs_hld(i,x);
        }
        out[x]=out_ord;
        out_ord++;
    }
    void solve(int root){
        dfs_sz(root,-1);
        dfs_hld(root,-1);
    }
    vector<pair<int,int>> query(int u,int v,bool edge){
        vector<pair<int,int>> ret;
        while(true){
            if(in[v]<in[u]) swap(u,v);
            if(head[u]==head[v]) break;
            ret.emplace_back(in[head[v]],in[v]+1);
            v=par[head[v]];
        }
        if(in[u]+edge!=in[v]+1) ret.emplace_back(in[u]+edge,in[v]+1);
        return ret;
    }
};

int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    ll n;
    cin >> n;
    ll sz=n+n-1;
    vector<ll> p(n,-1);
    tree g(sz);
    HLD h(sz);
    replr(i,1,n){
        cin >> p[i];
        p[i]--;
        g.add_edge(i,n+i-1);
        g.add_edge(p[i],n+i-1);
        h.add_edge(i,n+i-1);
        h.add_edge(p[i],n+i-1);
    }
    string s(n,'x');
    replr(i,1,n) cin >> s[i];
    ll q;
    cin >> q;
    vector<ll> u(q),v(q);
    rep(i,q){
        cin >> u[i] >> v[i];
        u[i]--;
        v[i]--;
    }
    g.LCA();
    rep(i,q){
        ll x=g.get_lca(u[i],v[i]);
        ll nu=u[i],nv=v[i];
        if(u[i]==x){
            nu=g.doubling(v[i],g.depth[v[i]]-g.depth[u[i]]-1);
        }else{
            nu=g.parent[0][u[i]];
        }
        if(v[i]==x){
            nv=g.doubling(u[i],g.depth[u[i]]-g.depth[v[i]]-1);
        }else{
            nv=g.parent[0][v[i]];
        }
        u[i]=nu;
        v[i]=nv;
    }
    h.solve(0);
    vector<ll> out(sz),ord(sz);
    rep(i,sz) out[i]=h.out[i];
    rep(i,sz) ord[out[i]]=i;
    vector<set<ll>> lr(sz);
    rep(i,q){
        if(out[v[i]]<out[u[i]]) swap(u[i],v[i]);
        lr[u[i]].insert(v[i]);
    }
    fenwick_tree<ll> BIT(sz);

    auto add=[&](ll a,ll b){
        auto now=h.query(a,b,false);
        for(auto [l,r]:now){
            BIT.add(l,r,1);
        }
    };

    auto get=[&](ll x){
        auto now=h.query(x,x,false);
        return BIT.get(now[0].first);
    };

    replr(i,1,n){
        if(s[i]=='#') add(n+i-1,n+i-1);
    }

    for(auto i:ord){
        if(!(get(i)&1)) continue;
        if(lr[i].empty()) continue;
        ll j=*begin(lr[i]);
        add(i,j);
        lr[i].erase(j);
        if(len(lr[i])<len(lr[j])){
            for(auto k:lr[i]) lr[j].insert(k);
        }else{
            for(auto k:lr[j]) lr[i].insert(k);
            swap(lr[i],lr[j]);
        }
        lr[i].clear();
    }

    replr(i,1,n){
        if(get(n+i-1)&1){
            cout << "No\n";
            return 0;
        }
    }
    cout << "Yes\n";
}