2020CCPC威海C. Rencontre，​ Clock Master​

2020CCPC威海C. Rencontre

对于三个点i,j,k和另外一个点p,dis(i,o)+dis(j,o)+dis(k,o)=1/2*(dis(i,j)+dis(i,k)+dis(j,k),还需要知道期望的性质：

这样期望就是E(f(i,j,k))=E(1/2*(dis(i,j)+dis(i,k)+dis(j,k)))=1/2*E(dis(i,j)+dis(i,k)+dis(j,k)),然后问题就转化成了求一个这样的式子

这个问题就是u1集合的各个点到u2集合的各个点的距离和;

之前有个题是计算树上任意两点的距离和，这个就是多了个分组的条件，那只要给siz多开一维表示组别就可以了

2020ccpc威海C Rencontre(树上各块中任意两点距离之和)_小菜鸡加油的博客-CSDN博客_c. rencontre

图片和思路来自2020CCPC(威海) - Rencontre(树形dp)_Frozen_Guardian的博客-CSDN博客

#include<bits/stdc++.h>
//#pragma-GCC-optimize("-Ofast");
#define int long long
//#define ll long long
#define lowbit(x) ((x)&(-x))
#define endl '\n'
#define pause system("pause")
using namespace std;
const int N=5e5+5;
const int mod=1e9+7;
const int inf=1e18;
const double pi=acos(-1);
int qpow(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
int getinv(int a){return qpow(a,mod-2);}
int head[N],cnt;
struct Edge
{
    int next,to,w;
}e[N];
void addedge(int from,int to,int w)
{
    e[++cnt].next=head[from];
    e[cnt].to=to;
    e[cnt].w=w;
    head[from]=cnt;
}
int n,siz[4][N],ct[4];
double ans=0;
void dfs1(int u,int fa)
{
    for(int i=head[u];i;i=e[i].next)
    {
        int j=e[i].to;
        if(j==fa) continue;
        dfs1(j,u);
        for(int k=1;k<=3;k++)
            siz[k][u]+=siz[k][j];
    }
}
void dfs2(int u,int fa)
{
    for(int i=head[u];i;i=e[i].next)
    {
        int j=e[i].to;
        if(j==fa) continue;
        dfs2(j,u);
        for(int k=1;k<=3;k++)
            for(int p=1;p<=3;p++)
        {
            if(k==p) continue;
            ans+=(siz[p][1]-siz[p][j])*(siz[k][j])*1.0*e[i].w/ct[p]/ct[k]/2.0;
        }
    }
}
signed main()
{
    //ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    cin>>n;
    for(int i=1;i<n;i++)
    {
        int u,v,w;
        cin>>u>>v>>w;
        addedge(u,v,w);
        addedge(v,u,w);
    }
    for(int i=1;i<=3;i++)
    {
        cin>>ct[i];
        for(int j=1;j<=ct[i];j++)
        {
            int x;cin>>x;
            siz[i][x]++;
        }
    }
    dfs1(1,0);
    dfs2(1,0);
    cout<<fixed<<setprecision(12)<<ans<<endl;
    return 0;
}

 Clock Master

 不同的质数的lcm是最大的，所以尽量选质数或者是质数的幂次，而且一个质数只要选一次就可以了，也就是在幂次里面选一个就可以，这就是分组背包了，预处理出log，跑一遍背包就可以

2020 CCPC 威海站 L Clock Master （数论+分组背包）_康宇的博客-CSDN博客

#include<bits/stdc++.h>
//#pragma-GCC-optimize("-Ofast");
#define int long long
//#define ll long long
#define lowbit(x) ((x)&(-x))
#define endl '\n'
#define pause system("pause")
using namespace std;
const int N=2e6+5;
const int mod=1e9+7;
const int inf=1e18;
const int up=1e6;
const double pi=acos(-1);
int qpow(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
int getinv(int a){return qpow(a,mod-2);}
int f[N],pr[N],cnt,ct[30005],a[30005][66];
double lg[30003],dp[300005];
bool isp[N];
void prime(int n)
{
    memset(isp,1,sizeof(isp));
    isp[0]=isp[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(isp[i]) pr[++cnt]=i;
        for(int j=1;j<=cnt&&i*pr[j]<=n;j++)
        {
            isp[i*pr[j]]=0;
            if(i%pr[j]==0) break;
        }
    }
}
signed main()
{
    //ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    prime(30000);
    for(int i=1;i<=cnt;i++)
    {
        int x=pr[i],y=pr[i];
        while(x<=30000)
        {
            a[y][++ct[y]]=x;
            x*=y;
        }
    }
    for(int i=1;i<=30000;i++) lg[i]=log(i);
    for(int i=1;i<=cnt;i++)
    {
        int x=pr[i];
        for(int j=30000;j>=0;j--)
            for(int k=1;k<=ct[x];k++)
                if(j>=a[x][k]) dp[j]=max(dp[j],dp[j-a[x][k]]+lg[a[x][k]]);
    }
    int T;
    cin>>T;
    while(T--)
    {
        int b;cin>>b;
        cout<<fixed<<setprecision(12)<<dp[b]<<endl;
    }
    return 0;
}