请问这样记忆化建立点分树为啥是错的呢?

P3806 【模板】点分治 1

Richard_Whr @ 2023-11-30 17:28:52

#include<bits/stdc++.h>
using namespace std;
const int N=1e4+10,M=N*2,S=1e7+10;
int h[N],e[M],ne[M],w[M],idx;
bool st[N];
int root[N];
int p[N],q[N];
bool f[S];
int n,m,k;
bool ans;

void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}

int get_size(int u,int fa)
{
    if(st[u]) return 0;
    int sz=1;
    for(int i=h[u];~i;i=ne[i])
    {
        int v=e[i];
        if(v==fa) continue;
        sz+=get_size(v,u);
    }   
    return sz;
} 

int get_wc(int u,int fa,int tot,int &wc)
{
    if(st[u]) return 0;
    int sum=1,ms=0;
    for(int i=h[u];~i;i=ne[i])
    {
        int v=e[i];
        if(v==fa) continue;
        int t=get_wc(v,u,tot,wc);
        sum+=t;
        ms=max(ms,t);
    }
    ms=max(ms,tot-sum);
    if(ms<=tot/2) wc=u;
    return sum;
}

void get_dist(int u,int fa,int dist,int &ql)
{
    if(st[u]||dist>k) return;
    q[ql++]=dist;
    for(int i=h[u];~i;i=ne[i])
    {
        int v=e[i];
        if(v==fa) continue;
        get_dist(v,u,dist+w[i],ql);
    }
}

void calc(int u)
{
    if(st[u]) return;

    if(root[u]) u=root[u];
    else get_wc(u,-1,get_size(u,-1),root[u]),u=root[u];

    st[u]=true;

    int pl=0;
    for(int i=h[u];~i;i=ne[i])
    {
        int v=e[i],ql=0;
        get_dist(v,u,w[i],ql);
        for(int j=0;j<ql;j++)
        {
            p[pl++]=q[j];
            if(q[j]==k) ans=true;
            if(f[k-q[j]]) ans=true;
        }
        for(int j=0;j<ql;j++) f[q[j]]=true;
    }

    for(int i=0;i<pl;i++) f[p[i]]=false;

    for(int i=h[u];~i;i=ne[i]) calc(e[i]);
}

int main()
{
    ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);

    freopen("P3806_1.in","r",stdin);    
    cin>>n>>m;
    memset(h,-1,sizeof h);
    for(int i=1,a,b,c;i<n;i++)
    {
        cin>>a>>b>>c;
        add(a,b,c),add(b,a,c);
    }

    while(m--)
    {
        cin>>k;
        memset(st,0,sizeof st),ans=false;
        calc(1);
        if(ans) cout<<"AYE"<<"\n";
        else cout<<"NAY"<<"\n";
    }

    return 0;
}

by 冷月葬T魂 @ 2023-11-30 20:25:19

@Richard_Whr

for(int i=h[u];~i;i=ne[i])
{
    int v=e[i],ql=0;

是否应该加上 if(st[v]) continue;


by Richard_Whr @ 2023-11-30 20:37:38

@冷月葬T魂 应该不用吧,get_dist函数中遇到 st[u]=true的情况会自动跳过,不会修改


by zzzyyyyhhhhh @ 2024-03-30 19:34:08

%%%


|