Codeforces Round #614 选讲

http://codeforces.com/contest/1292/problem/C

注意到编号x的边对答案要有贡献，必须和0到x-1的边一起形成一条链，否则x及编号比x大的边都没有贡献。由此，对答案有贡献的边形成了一条链，且这条链的编号是个谷形，即中间编号小，往两边编号变大，编号最大的边在最外侧。由此可以进行dp，dp[u][v]表示如果上述链为点u到点v这条链的答案。令sz[u][v]为以u为根，子树v的大小；fa[u][v]为以u为根，点v的父亲，则有dp[u][v]=dp[v][u]=sz[u][v]*sz[v][u]+max(dp[u][fa[u][v]],dp[v][fa[v][u]])。（考虑放入编号最大的一条边，相当于某些路径的mex增加了1）

#include <bits/stdc++.h>
using namespace std;
#define LL long long

int n;
#define maxn 3011
struct Edge{int to,next;}edge[maxn<<1]; int first[maxn],le=2;
void in(int x,int y) {Edge &e=edge[le]; e.to=y; e.next=first[x]; first[x]=le++;}
void insert(int x,int y) {in(x,y); in(y,x);}

int sz[maxn][maxn],fa[maxn][maxn];
LL dp[maxn][maxn]; 
struct NODE{int x,y;};
vector<NODE> vec[maxn];
void dfs(int x,int f,int top,int dep)
{
    sz[top][x]=1; fa[top][x]=f;
    if (dep>0) vec[dep].push_back((NODE){top,x});
    for (int i=first[x];i;i=edge[i].next)
    {
        Edge &e=edge[i]; if (e.to==f) continue;
        dfs(e.to,x,top,dep+1); sz[top][x]+=sz[top][e.to];
    }
}

int main()
{
    scanf("%d",&n);
    for (int i=1,x,y;i<n;i++) {scanf("%d%d",&x,&y); insert(x,y);}
    for (int i=1;i<=n;i++) dfs(i,0,i,0);
    
    for (int len=1;len<n;len++)
        for (int i=0,to=vec[len].size();i<to;i++)
        {
            int a=vec[len][i].x,b=vec[len][i].y;
            dp[a][b]=dp[b][a]=sz[a][b]*sz[b][a]+max(dp[a][fa[a][b]],dp[b][fa[b][a]]);
        }
    
    LL ans=0;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            ans=max(ans,dp[i][j]);
    printf("%lld
",ans);
    return 0;
}

 http://codeforces.com/contest/1292/problem/D

回忆找重心的过程：从根节点开始，看点权和最大的子树，其点权和是否超过总的一半，是则在该子树中找重心，否则当前点为重心。以1为根时，在深度dep的节点x的子树中的点y必须满足：y的质因子中最大的dep个乘积即为x。由此，从1出发，对所有数进行分解后（相同的合并），每次取每个数的最大的质因子p，算出每个p的点权和，然后看重儿子的点权和是否超过n/2，若否则找到答案，若是则调整答案，并删去该子树中的最大质因子（将其次数-1），然后在该子树中循环进行该过程。

#include <bits/stdc++.h>
using namespace std;
#define LL long long

int n;
#define maxn 5011
#define MAXNUM 5000
int cnt[maxn],fra[maxn][maxn],poi[maxn],cntpri[maxn],pri[maxn],lp; bool vis[maxn];
void makeprime(int n)
{
    for (int i=2;i<=n;i++)
    {
        if (!vis[i]) pri[++lp]=i;
        for (int j=1;j<=lp && 1ll*i*pri[j]<=n;j++)
        {
            vis[i*pri[j]]=1;
            if (i%pri[j]==0) break;
        }
    }
}

int main()
{
    scanf("%d",&n);
    for (int i=1,x;i<=n;i++) scanf("%d",&x),cnt[x]++;
    makeprime(MAXNUM);
    
    LL ans=0;
    for (int x=0;x<=MAXNUM;x++)
        for (int i=1;i<=lp && pri[i]<=x;i++)
        {
            for (int j=pri[i];j<=x;j*=pri[i])
                fra[x][i]+=x/j;
            ans+=1ll*cnt[x]*fra[x][i];
        }
        //trick: The number of prime 2 in k! is k/2+k/4+k/8+...
    
    for (int x=0;x<=MAXNUM;x++) poi[x]=lp;
    while (1)
    {
        for (int i=1;i<=lp;i++) cntpri[i]=0;
        for (int x=0;x<=MAXNUM;x++)
        {
            while (poi[x]>0 && fra[x][poi[x]]==0) poi[x]--;
            if (poi[x]>0) cntpri[poi[x]]+=cnt[x];
        }
        bool findans=1;
        for (int i=1;i<=lp;i++) if (cntpri[i]*2>n)
        {
            findans=0;
            ans-=cntpri[i]-(n-cntpri[i]);
            for (int x=0;x<=MAXNUM;x++) if (poi[x]>0 && poi[x]==i) fra[x][poi[x]]--; else poi[x]=0;
            break;
        }
        if (findans) break;
    }
    printf("%lld
",ans);
    return 0;
}

。