记一次改题

最近考试题一直改不完，不过还是想写一篇博客

[NOIP模拟测试26（A）T2]

下午2：00 下定决心改过T2（跟他们说复习一下树剖板子）

下午3：00终于看懂原理，开始码

下午3：01突然发现考试代码改一改可以成标准的暴力

下午3：30发现暴力可以得到87分，于是动了邪念。。。。（决定用分块艹过）

期间请教死缠了skyh和remarkable 大佬，几次产生放弃念头（唉，我还是打正解吧）

8：30尝试更改块的大小，终于AC。

学到的东西：

1.分块万能。

自己YY的树上分块：

//预处理
void dfs(const int x,const int kp){
    dfn[x]=++tot;
    if(t&&d[x]%t==0)ed[x]=1;
    tp[x]=kp;
    for(int i=head[x],y;i;i=nxt[i])
    {
        y=to[i];
        d[y]=d[x]+1;
        if(ed[x])tp[y]=y,er[y].set(a[y]),dfs(y,y);
        else tp[y]=tp[x],er[y]=er[x],er[y].set(a[y]),dfs(y,kp);
    }
}
//询问
for(register int i=1;i<=c;++i)
{
    register int now=p[i];
    while(fa[tp[now]]&&d[fa[tp[now]]]>=d[lc])
    {
        s[i]|=er[now];
        now=fa[tp[now]];
    }
    while(1)
    {
        s[i].set(a[now]);
        if(now==lc)break;
        now=fa[now];
    }
}

2.手打bitset

 

struct bitset{
    unsigned long long bit[16];
    inline void reset(){
        for(register int i=0;i^16;i+=4) bit[i]=bit[i|1]=bit[i|2]=bit[i|3]=0;
    }
    inline friend void operator |=(bitset &a,const bitset &b){
        for(register int i=0;i^16;i+=4) a.bit[i]|=b.bit[i],a.bit[i|1]|=b.bit[i|1],a.bit[i|2]|=b.bit[i|2],a.bit[i|3]|=b.bit[i|3];
    }
    inline friend bitset operator |(const bitset &a,const bitset &b){
        bitset ans;
        for(register int i=0;i^16;i+=4)            ans.bit[i]=a.bit[i]|b.bit[i],ans.bit[i|1]=a.bit[i|1]|b.bit[i|1],ans.bit[i|2]=a.bit[i|2]|b.bit[i|2],ans.bit[i|3]=a.bit[i|3]|b.bit[i|3];
        return ans;
    }
    inline void set(const int x){
        bit[x/64]|=(unsigned long long)1<<(x&63);
    }
    inline int count(){
        register int cnt=0;
        for(int i=0;i<16;++i)
            cnt+=ct[bit[i]&65535]+ct[(bit[i]>>16)&65535]+ct[(bit[i]>>32)&65535]+ct[(bit[i]>>48)&65535];
        return cnt;
    }
}; 

 

 

 

3.适时调整块的大小

4.多想一些东西。

AC代码：

#include<cstring>
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int L=1<<20|1;
char buf[L],*S,*T;
#define getchar() ((S==T&&(T=(S=buf)+fread(buf,1,L,stdin),S==T))?EOF:*S++)
int dfn[300005],f[300005][20];
int t;
int head[300005],to[300005],nxt[300005],cnt,p[6],a[300005],nk[6],c,ans,n,m,q,tot,d[300005],tcnt[105],yk[105],maxd,tp[300005],ct[65600],fa[300005];
char ed[300005];
struct bitset{
    #define lowbit(x) (x&-x)
    unsigned long long bit[16];
    inline void reset(){
        for(register int i=0;i^16;i+=4) bit[i]=bit[i|1]=bit[i|2]=bit[i|3]=0;
    }
    inline friend void operator |=(bitset &a,const bitset &b){
        for(register int i=0;i^16;i+=4) a.bit[i]|=b.bit[i],a.bit[i|1]|=b.bit[i|1],a.bit[i|2]|=b.bit[i|2],a.bit[i|3]|=b.bit[i|3];
    }
    inline friend bitset operator |(const bitset &a,const bitset &b){
        bitset ans;
        for(register int i=0;i^16;i+=4) 
            ans.bit[i]=a.bit[i]|b.bit[i],ans.bit[i|1]=a.bit[i|1]|b.bit[i|1],ans.bit[i|2]=a.bit[i|2]|b.bit[i|2],ans.bit[i|3]=a.bit[i|3]|b.bit[i|3];
        return ans;
    }
    inline void set(const int x){
        bit[x/64]|=(unsigned long long)1<<(x&63);
    }
    inline int count(){
        register int cnt=0;
        for(int i=0;i<16;++i)
            cnt+=ct[bit[i]&65535]+ct[(bit[i]>>16)&65535]+ct[(bit[i]>>32)&65535]+ct[(bit[i]>>48)&65535];
        return cnt;
    }
}s[6],ak[32],er[300005];
inline void Add(const int u,const int v)
{
    to[++cnt]=v;
    nxt[cnt]=head[u];
    head[u]=cnt;
    return ;
}
void dfs(const int x,const int kp){
    dfn[x]=++tot;
    if(t&&d[x]%t==0)ed[x]=1;
    tp[x]=kp;
    for(int i=head[x],y;i;i=nxt[i])
    {
        y=to[i];
        d[y]=d[x]+1;
        if(ed[x])tp[y]=y,er[y].set(a[y]),dfs(y,y);
        else tp[y]=tp[x],er[y]=er[x],er[y].set(a[y]),dfs(y,kp);
    }
}
int Q[300005],l,r;
inline int lca(int a,int b)
{
    if(d[a]>d[b])a^=b^=a^=b;
    for(int t=18;t>=0;t--)
        if(f[b][t]&&d[a]<=d[f[b][t]])b=f[b][t];
    if(a==b)return a;
    for(int t=18;t>=0;t--)
        if(f[a][t]&&f[a][t]!=f[b][t])a=f[a][t],b=f[b][t];
    return f[a][0];
}
inline int read(register int x=0,register char ch=getchar()){
    while(ch<48||ch>57) ch=getchar();
    while(ch>=48&&ch<=57) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    return x;
}
char cp[10],top;
inline void write(int x){
    do{
        cp[top++]=x%10,x/=10;
    }while(x);
    do{
        putchar(cp[--top]+48);
    }while(top);
    putchar(10);
}
int main()
{
    n=read(); m=read(); q=read();
    for(int i=1;i<=65535;++i) ct[i]=ct[(i&-i)^i]+1;
    yk[1]=1,yk[2]=2,yk[4]=3,yk[8]=4,yk[16]=5;
    for(int i=1;i<=31;++i)
        for(int j=1;j<=5;++j)
            if(i&(1<<(j-1)))tcnt[i]++;
    for(int i=2;i<=n;Add(fa[i],i),++i) fa[i]=read();
    for(int i=1;i<=n;++i) a[i]=read();
    t=pow(n,0.555);
    dfs(1,1);
    for(int j=1;j<=n;++j){
        f[j][0]=fa[j];
        for(int i=0;i<18;++i)
            f[j][i+1]=f[f[j][i]][i];
    }
    while(q--)
    {
        c=read(); ans=0x7f7f7f7f;
        for(register int i=1;i<=c;++i) p[i]=read();
        register int now=-1,pa,pb;
        for(register int i=1;i<=c;++i)
            for(register int j=i+1;j<=c;++j)
                if(abs(dfn[p[i]]-dfn[p[j]])>now)pa=p[i],pb=p[j],now=abs(dfn[p[i]]-dfn[p[j]]);
        register int lc=lca(pa,pb);
        for(register int i=1;i<=c;++i)
        {
            register int now=p[i];
            while(fa[tp[now]]&&d[fa[tp[now]]]>=d[lc])
            {
                s[i]|=er[now];
                now=fa[tp[now]];
            }
            while(1)
            {
                s[i].set(a[now]);
                if(now==lc)break;
                now=fa[now];
            }
        }
        for(register int i=1;i<1<<c;++i)
        {
            ak[i]=ak[i^(i&-i)]|s[yk[i&-i]];
            ans=min(ans,(int)ak[i].count()/tcnt[i]);
        }
        for(register int i=1;i<=c;++i)s[i].reset();
        write(ans*c);
    }
    return 0;
}