[BJOI2014]想法(随机算法,神奇思路,拓扑排序)
对于这种随机数据或者随机算法的题……
都是神仙题吧。
要求的就是对每个点前 \(m\) 个点中有多少个可以到达它。
由于评分方式这么奇怪,不妨考虑随机。
随机 127 次(可以选别的数,够多而且不 T 就行),每次给前 \(m\) 个数随机赋值,然后拓扑求出能到达每个点的最小值。
可能脸黑,所以多跑几次取平均数。最后每个点的平均最小值就可以看成真的期望最小值。
有一个结论:\([0,v]\) 中取 \(x\) 个数,最小值的期望值是 \(\frac{v}{x+1}\)。
所以就能算出每个点的最小值由多少个数取得了。也就是答案。
#include<bits/stdc++.h>
using namespace std;
const int maxn=2000200;
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
int x=0,f=0;char ch=getchar();
while(ch<'0' || ch>'9') f|=ch=='-',ch=getchar();
while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
return f?-x:x;
}
int n,m,mn[maxn],pre[maxn][2];
double sum[maxn];
int main(){
srand(20050818);
n=read();m=read();
FOR(i,1,n-m){
int x=read(),y=read();
pre[m+i][0]=x;
pre[m+i][1]=y;
}
FOR(i,1,127){
FOR(j,1,m) mn[j]=rand();
FOR(j,m+1,n) mn[j]=min(mn[pre[j][0]],mn[pre[j][1]]);
FOR(j,1,n) sum[j]+=mn[j]/127.0;
}
FOR(j,m+1,n) printf("%.0lf\n",RAND_MAX/sum[j]-1);
}