2017.10.23一试

不鬼扯了……
今天的题我不知道从哪扯起……
好歹此次的题已经从个人层面上升到国家层面了嘛，对不对？
N ational O lympiad in I nformatics for P enguin
看看，人家比我还会扯。

——所以跳过这欢乐的环节。

T1

*签到题。
斐波那契数列递增地很快，所以在10^9前只有40来个数。
于是乎此题便在半打表的状态下AC。
代码如下：

#include<bits/stdc++.h>
#define fo(i,a,b) for(register int i=a;i<=b;i++)
using namespace std;

int fib[47]={0,1,1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,1597,2584,4181,6765,10946,17711,28657,46368,75025,121393,196418,317811,514229,832040,1346269,2178309,3524578,5702887,9227465,14930352,24157817,39088169,63245986,102334155,165580141,267914296,433494437,701408733,1134903170};

inline int read()
{
    int X=0,w=1; char ch=0;
    while(ch<'0' || ch>'9') {if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0' && ch<='9') X=(X<<3)+(X<<1)+ch-'0',ch=getchar();
    return X*w;
}

inline void write(int x)
{
     if(x<0) putchar('-'),x=-x;
     if(x>9) write(x/10);
     putchar(x%10+'0');
}

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

//  freopen("fib.in","r",stdin);
//  freopen("fib.out","w",stdout);

    int T;
    long long a;

    T=read();
    while(T--)
    {
        a=read();
        bool t=false;

        fo(i,0,47)
        {
            fo(j,0,47)
            {
                if(a==fib[i]*fib[j])
                {
                    cout<<"Yes"<<endl;
                    t=true;
                    break;
                }

                if(a<fib[i]*fib[j])
                break;
            }

            if(t)
            break;
        }

        if(!t)
        cout<<"No"<<endl;
    }

    return 0;
}

T2

低阶斯特拉算法全部爆炸。
所以还是LCA的好。
zz……
30分暴力都没有……

还有linux下的30分到OJ上交就全部AC我也是心悦诚服。
代码如下：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
inline int  read(){
    int  i=0,f=1;
    char ch;
    for(ch=getchar();!isdigit(ch);ch=getchar())
        if(ch=='-') f=-1;
    for(;isdigit(ch);ch=getchar())
        i=(i<<3)+(i<<1)+(ch^48);
    return i*f;
}
struct node
{
    int w,to,next;

}bian[200050];
int first[100050];
int tot;
int n;
int m;
int size[100050];
bool bz[100050];
int dep[100050];
int fa[100050][21];
int dis[100050];
bool bzz[100050];
int a,b;
void add(int x,int y,int z)
{
    tot++;
    bian[tot].next=first[x];
    first[x]=tot;
    bian[tot].to=y;
    bian[tot].w=z;
}
void dfs(int r)
{
    bz[r]=1;
    size[r]=1;
    for(int u=first[r];u;u=bian[u].next)
    {
        if(bz[bian[u].to]==0)
        {
            dfs(bian[u].to);
            size[r]+=size[bian[u].to];
        }
    }
}
inline void ddfs(int x)
{
    bzz[x]=true;
    for(int i=first[x];i;i=bian[i].next)
    {
        int v=bian[i].to;
        if(!bzz[v])
        {
            dep[v]=dep[x]+1;
            dis[v]=dis[x]+bian[i].w;
            fa[v][0]=x;ddfs(v);
        }
    }
}
int lca(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    int len=dep[x]-dep[y];
    for(int i=20;i>=0;--i)
    {
        if(len>=(1<<i))
        {
            len-=(1<<i);
            x=fa[x][i];
        }
    }
    if(x==y) return x;
    for(int i=20;i>=0;--i)
    {
        if(fa[x][i]!=fa[y][i])
        {
            x=fa[x][i];
            y=fa[y][i];
        }
    }
    return fa[x][0];
}



int lcaa(int x,int y)
{
    for(int i=20;i>=0;--i)
    {
        if(fa[x][i]!=fa[y][i])
        {
            x=fa[x][i];
            y=fa[y][i];
        }
    }
    return x;
}




int main()
{
    //freopen("equal.in","r",stdin);
    //freopen("equal.out","w",stdout);
    n=read();
    for(int i=1;i<n;i++)
    {
        int x;int y;
        x=read();
        y=read();
        add(x,y,1);
        add(y,x,1);
    }
    dfs(1);
    ddfs(1);
    for(int i=1;i<=20;++i)
        for(int j=1;j<=n;++j)
            if(fa[j][i-1]) fa[j][i]=fa[fa[j][i-1]][i-1];

    m=read();
    for(int i=1;i<=m;i++)
    {
        a=read();
        b=read();
        if(a==b)
        {
            cout<<n<<endl;
            continue;
        }
        int kkk=dis[a]+dis[b]-2*dis[lca(a,b)];

        if(kkk%2==1)
        {
            cout<<"0"<<endl;
            continue;
        }
        else
        {
            if(dep[a]==dep[b])
            {
                cout<<(n-size[lcaa(a,b)]-size[lcaa(b,a)])<<endl;
            }
            else
            {
                kkk=kkk/2;
                if(dep[a]<dep[b])swap(a,b);

                if(kkk==1)
                {
                    cout<<(size[fa[a][0]]-size[a])<<endl;
                }
                else
                {
                    kkk=kkk-1;
                    for(int i=20;i>=0;--i)
                    {
                        if(kkk>=(1<<i))
                        {
                            kkk-=(1<<i);
                            a=fa[a][i];
                        }
                    }
                    cout<<(size[fa[a][0]]-size[a])<<endl;
                }
            }
        }

    }

    return 0;
}

T3

大暴力都只有30分诶……
而判断一下奇偶然后输出n/2或（n+1）/2就有20分了。
虎虎虎啊哈哈呵呵呵。
我呀~听到了~嘞个~流动的~测试数据~
汹涌~啊嘞个~澎湃

而此题是一道DP。
一条边可以用两只企鹅站，这样的一条点对，越多越好。
如果是 ans对点，ans*2≥k，那么只需要 (k+1)/2 条边。
否则，需要 ans + (k-ans*2) 条边。
现在问题就转为求这样的点对有多少。
dp[i][0]表示以 i 为根的子树中能够两两配对的最大点数，不包含节点 i。
dp[i][1]表示以 i 为根的子树中能够两两配对的最大点数，包含节点 i。
转移方程：
dp[u][0]=Σv是 u的儿子 dp[v][1];
dp[u][1]=max(dp[u][1],dp[u][0]-dp[v][1]+dp[v][0]+2)。
最后 max(dp[1][0],dp[1][1])就是ans了。

*以上为题解。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
inline int  read()
{
    int  i=0,f=1;
    char ch;
    for(ch=getchar();!isdigit(ch);ch=getchar())
        if(ch=='-') f=-1;
    for(;isdigit(ch);ch=getchar())
        i=(i<<3)+(i<<1)+(ch^48);
    return i*f;
}

int T;
int n,k;
struct node
{
    int to,next;

}bian[200050];
int first[100050];
int tot;
int z[100050];
bool bz[100050];
int ans;
int d[100050];
void add(int x,int y)
{
    tot++;
    bian[tot].next=first[x];
    first[x]=tot;
    bian[tot].to=y;
}
void dfs(int r)
{
    z[r]=1;
    bz[r]=1;
    for(int u=first[r];u;u=bian[u].next )
    {
        if(bz[bian[u].to]==0)
        {
            dfs(bian[u].to);
            z[r]+=z[bian[u].to];
        }
    }
    if(z[r]>=2)
    {
        d[2]++;

        z[r]=0;
    }
}

int main()
{
//  freopen("tree.in","r",stdin);
//  freopen("tree.out","w",stdout);
    T=read();
    while(T--)
    {
        for(int i=1;i<=tot;i++)
        {
            bian[i].to=0;
            bian[i].next=0;
        }
        tot=0;
        memset(d,0,sizeof(d));
        memset(first,0,sizeof(first));
        memset(bz,0,sizeof(bz));
        memset(z,0,sizeof(z));
        ans=0;
        n=read();
        k=read();
        int x;
        for(int i=1;i<n;i++)
        {
            x=read();
            add(x,i+1);
            add(i+1,x);
        }
        dfs(1);
        while(d[2]!=0)
        {
            if(k>0)
            {
                ans++;
            }
            k=k-2;
            d[2]--;
        }
        if(k<=0)
        cout<<ans<<endl;
        if(k>0)
        {
            cout<<ans+k<<endl;
        }

    }
    return 0;
}

最后の总结（感觉意思不太对）：

看准T1
认真暴力
培养包容坚毅、德嘉学优、健康快乐、国际领先的中国人。（啊吔）

以上！