[BJOI2014]想法(随机算法,神奇思路,拓扑排序)

对于这种随机数据或者随机算法的题……

都是神仙题吧。

要求的就是对每个点前 \(m\) 个点中有多少个可以到达它。

由于评分方式这么奇怪,不妨考虑随机。

随机 127 次(可以选别的数,够多而且不 T 就行),每次给前 \(m\) 个数随机赋值,然后拓扑求出能到达每个点的最小值。

可能脸黑,所以多跑几次取平均数。最后每个点的平均最小值就可以看成真的期望最小值。

有一个结论:\([0,v]\) 中取 \(x\) 个数,最小值的期望值是 \(\frac{v}{x+1}\)

所以就能算出每个点的最小值由多少个数取得了。也就是答案。

#include
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);
}

转载于:https://www.cnblogs.com/1000Suns/p/11295973.html

你可能感兴趣的:([BJOI2014]想法(随机算法,神奇思路,拓扑排序))