记codevs第一次月赛

第一次参加这种有奖励的比赛（没错，我就是为猴子而去的

一年没怎么碰代码果然手生，还是用没写多久的C++，差点全跪了

T1数学奇才琪露诺：

首先定义一个函数F(x)，F(x)=x的各个数位上的数字和

然后在区间[l,r]求F(x)^k*p+q=x的所有x，按升序输出

T1题解：

枚举x肯定是不行的，F(x)的值只有0到81，我们枚举F(x)的值，然后算出F(x)^k*p+q，看数位和是不是符合

可能会爆int，所以要开long long

#include <cstdio>
#include <algorithm>
#include <vector>
#include <set>
using namespace std;
 
long long k,p,q,l,r,a[],tot;
 
long long s(long long x){
    long long ss=;
    while(x>){
        ss+=x%;
        x=x/;
    }
    return ss;
}
 
long long pow(long long x,long long y){
    long long ipow=;
    long long i;
    for(i=;i<=y;++i)ipow=ipow*x;
    return ipow;
}
 
int main(){
    scanf("%lld%lld%lld%lld%lld",&k,&p,&q,&l,&r);
    long long i;
    tot=;
    for(i=;i<=;++i){
        long long tmp;
        tmp=pow(i,k)*p+q;
        if((tmp>=)&&(s(tmp)==i)&&(tmp>=l)&&(tmp<=r))a[++tot]=tmp;
    }
    sort(a+,a++tot);
    printf("%lld\n",tot);
    for(i=;i<=tot;++i)printf("%lld ",a[i]);
    return ;
}

T1

T2完美拓印：

给你一条轮廓线和一个印章，问有多少种方法可以让印章的一条边缘（上下两边）与轮廓线重合，印章可以180°旋转

T2题解：

我们注意到只要相邻两格的高度差一样就行，所以我们把相邻两项的差拿出来做kmp就行了

注意要做四次kmp，下面平的也要考虑，还有就是特判n=1的情况

#include <cstdio>
#include <algorithm>
#include <vector>
#include <set>
using namespace std;
 
int n,m,a[],b[],nex[],ans;
 
int work(){
    int i,j;
    for(i=n-;i>;--i)a[i]=a[i-]-a[i];
    nex[]=;
    nex[]=;
    for(i=;i<n;++i){
        j=nex[i-];
        while((j>)&&(a[i]!=a[j+]))j=nex[j];
        if(a[i]==a[j+])nex[i]=j+;
        else nex[i]=;
    }
    i=;j=;
    while(j<m-){
        if(a[i+]==b[j+])++i,++j;
        else
        if(i==)++j;
        else i=nex[i];
        if(i==n-){
            ++ans;
            i=nex[i];
        }
    }
    for(i=;i<n;++i)a[i]=a[i-]-a[i];
    return ;
}
 
int swap(int &a,int &b){
    int t;
    t=a;a=b;b=t;
    return ;
}
 
int main(){
    int i;
    ans=;
    scanf("%d%d",&n,&m);
    for(i=;i<n;++i)scanf("%d",&a[i]);
    for(i=;i<m;++i)scanf("%d",&b[i]);
    if(n==){
        printf("%d\n",*m);
        return ;
    }
    for(i=m-;i>;--i)b[i]=b[i-]-b[i];
    work();
    for(i=;i<n/;++i)swap(a[i],a[n-i-]);
    for(i=;i<n;++i)a[i]=-a[i];
    work();
    for(i=;i<n;++i)a[i]=;
    work();
    work();
    printf("%d\n",ans);
    return ;
}

T2

T3幻影阁的难题：

给你两棵树，你可以分别在两棵树上找一个点，然后在这两个点之间连一条长度为t的边

1.求连边之后最长路的最小值

2.求最长路的期望长度

T3题解：

首先我们看第一问，我们只要树dp算出从某个点出发的最长链，然后取两棵树的最小值相加再加t，但是我们还要考虑不经过新边的情况，那么就是原来树上的最长链，从这两个中取最大值就行了

第二问要求期望，所以我们要把所有的情况都算出来，然后总长除以n*m，我们先把两棵树计算好的最长链排好序

然后我们要计算的就是∑max(最长链,a[i]+b[j])，因为我们排好了序，所以对于每一个i，取最长链的是一段连续的区间，我们只要记一下前缀和就可以快速计算了

#include <cstdio>
#include <algorithm>
#include <vector>
#include <set>
using namespace std;
 
long long n,m,t,tot;
long long nex[],las[],l[],aa[];
long long s1[],s2[],fa[],fir[],q[],len[];
bool f[];
long long ss2[],g[];
long long ans,ans1,ans2;
 
int insert(int x,int y,int z){
    ++tot;
    l[tot]=z;
    las[tot]=y;
    nex[tot]=fir[x];
    fir[x]=tot;
}
 
int work(){
    int i,x,y,z;
    tot=;
    for(i=;i<=n;++i)fir[i]=,f[i]=false,s1[i]=,s2[i]=;
    for(i=;i<n;++i){
        scanf("%d%d%d",&x,&y,&z);
        insert(x,y,z);
        insert(y,x,z);
    }
    f[]=true;
    int ll,rr;
    ll=;rr=;
    q[]=;
    for(;ll<=rr;++ll){
        int j;
        for(j=fir[q[ll]];j!=;j=nex[j])
        if(!f[las[j]]){
            f[las[j]]=true;
            q[++rr]=las[j];
            len[rr]=l[j];
            fa[rr]=ll;
        }
    }
    len[]=;g[]=;
    s1[]=;s2[]=;
    for(i=n;i>;--i){
        if(s1[i]+len[i]>s1[fa[i]])s2[fa[i]]=s1[fa[i]],s1[fa[i]]=s1[i]+len[i];
        else
        if(s1[i]+len[i]>s2[fa[i]])s2[fa[i]]=s1[i]+len[i];
    }
    long long tmp=;
    for(i=;i<=n;++i){
        if(s1[i]+len[i]!=s1[fa[i]])g[i]=s1[fa[i]]+len[i];
        else g[i]=s2[fa[i]]+len[i];
        g[i]=max(g[i],g[fa[i]]+len[i]);
        tmp=min(tmp,max(g[i],s1[i]));
        ans1=max(ans1,g[i]+s1[i]);
    }
    ans2+=tmp;
    return ;
}
 
long long gcd(long long a,long long b){
    if(b==)return a;
    return gcd(b,a%b);
}
 
int main(){
    scanf("%d%d%d",&n,&m,&t);
    ans=;ans1=;ans2=;
    work();
    int i,j;
    for(i=;i<=n;++i)aa[i]=max(g[i],s1[i]);
    swap(n,m);
    work();
    for(i=;i<=n;++i)g[i]=max(g[i],s1[i]);
    sort(aa+,aa++m);
    sort(g+,g++n);
    ans1-=t;
    ss2[]=aa[];
    for(i=;i<=m;++i)ss2[i]=ss2[i-]+aa[i];
    j=m;
    for(i=;i<=n;++i){
        while((j>)&&(g[i]+aa[j]>ans1))--j;
        ans+=ans1*j+g[i]*(m-j)+ss2[m]-ss2[j];
    }
    ans+=t*n*m;
    ans2=max(ans2+t,ans1+t);
    printf("%lld\n",ans2);
    long long tmp;
    tmp=gcd(ans,n*m);
    printf("%lld",ans/tmp);
    printf("/");
    printf("%lld",n*m/tmp);
    return ;
}

T3