开启掘金成长之旅！这是我参与「掘金日新计划 · 12 月更文挑战」的第1天，点击查看活动详情
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;
}
