NOIP2017金秋冲刺训练营杯联赛模拟大奖赛第二轮Day2题解

　　肝了两题...

　　T1一眼题，分解质因数，找出2的个数和5的个数取min输出

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
int n,m,T;
int fac2[maxn],fac5[maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
int main()
{
    freopen("zero.in","r",stdin);
    freopen("zero.out","w",stdout);
    read(T);
    for(int i=;i<=;i++)
    {
        int x=i;fac2[i]=fac2[i-];fac5[i]=fac5[i-];
        for(;(x&)==;x>>=)fac2[i]++;
        for(;x%==;x/=)fac5[i]++;
    }
    while(T--)
    {
        read(m);read(n);
        int cnt2=fac2[m]-fac2[n]-fac2[m-n];
        int cnt5=fac5[m]-fac5[n]-fac5[m-n];
        printf("%d\n",min(cnt2,cnt5));
    }
    return ;
}

　　T2是个环套树题

　　分三类：①两个点都在同个树上，求直径，两次bfs或者树形DP即可

　　　　　　②一个点在环上一个点在树上，答案为最大深度dep+len/2

　　　　　　③两个点在不同的树上，即求max(depi+depj+dis(i,j))，而dis(i,j)显然是<=len/2且递增的，用单调队列维护max(depj+dis(i,j))，队头与i距离超过len/2就出队，走一圈之后就可以算出这种情况的答案

　　三类取max

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
struct poi{int too,pre;}e[maxn];
int n,m,x,y,z,tot,ans,st,ed,cnt,len;
int mx1[maxn],mx2[maxn],d[maxn],last[maxn],isrt[maxn],rt[maxn],q[maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
inline void add(int x,int y){e[++tot].too=y;e[tot].pre=last[x];last[x]=tot;}
void dfs(int x,int fa)
{
    d[x]=d[fa]+;mx1[x]=d[x];mx2[x]=;
    for(int i=last[x];i;i=e[i].pre)
    if(e[i].too!=fa&&!isrt[e[i].too])
    {
        dfs(e[i].too,x);
        if(mx1[e[i].too]>mx1[x])mx2[x]=mx1[x],mx1[x]=mx1[e[i].too];
        else if(mx1[e[i].too]>mx2[x])mx2[x]=mx1[e[i].too];
    }
    ans=max(ans,mx1[x]+mx2[x]-(d[x]<<));
}
int dfs2(int x,int fa)
{
    if(x==ed)return len++,;
    for(int i=last[x];i;i=e[i].pre)
    if(e[i].too!=fa)if(dfs2(e[i].too,x))
    return rt[++cnt]=x,isrt[x]=,len++,;
    return ;
}
inline int dis(int x,int y){if(x>=y)return x-y;return len-(y-x);}
inline void work()
{
    int l=,r=;
    for(int i=(len>>);i;q[++r]=i,i--)
    while(l<=r&&mx1[rt[q[r]]]+dis(q[r],cnt)<=mx1[rt[i]]+dis(i,cnt))r--;
    for(int i=cnt;i;i--)
    {
        while(l<=r&&(dis(q[l],i)>(len>>)||q[l]==i))l++;
        ans=max(ans,mx1[rt[i]]+mx1[rt[q[l]]]+dis(q[l],i));
        while(l<=r&&mx1[rt[q[r]]]+dis(q[r],i)<=mx1[rt[i]])r--;
        q[++r]=i;
    }
}
int main()
{
    freopen("road.in","r",stdin);
    freopen("road.out","w",stdout);
    read(n);read(m);
    for(int i=;i<n;i++)read(x),read(y),add(x,y),add(y,x);
    for(int i=;i<=m;i++)
    {
        read(st);read(ed);ans=len=cnt=;
        memset(isrt,,sizeof(isrt));d[]=-;
        rt[++cnt]=ed;isrt[ed]=;dfs2(st,);
        for(int j=;j<=cnt;j++)dfs(rt[j],),ans=max(ans,mx1[rt[j]]+(len>>));
        work();printf("%d\n",ans);
    }
    return ;
}