【JZOJ 3765】【BJOI2014】想法

Description

小强和阿米巴是好朋友。

小强要出一套题目。他的题目以涉及面广（偏）、考察深入（怪）、思维强度大（难）著称。他为了出题，一共攒了M个本质不同的想法，每个想法形成了一个题目。不过，他觉得拿这些题目去考察选手会把比赛搞的太过变态，所以，想请阿米巴来帮忙调整一下他的题目。

阿米巴指出，为了让一场考试的题目的考察点尽量全面，有一个通用的做法叫做“组合”。如果把两个题目A和B组合在一起，那么组合而成的题目涉及到的想法的集合就是A涉及到的想法的集合和B涉及到的想法的集合的并。

并且，题目是可以反复组合的。

例如，小强现在有三个想法1,2,3，分别对应了题目P1,P2,P3。

现在，小强把P1和P2组合得到P4。P4涉及的想法的集合是{1,2}。

之后，小强把P2和P3组合得到P5。P5涉及的想法的集合是{2,3}。

最后，小强把P4和P5组合得到P6。P6涉及的想法的集合是{1,2,3}。

现在，小强告诉你每个题目都是如何组合而来的。你要回答的就是，每个题目涉及的想法的集合有多大。

不过，这个问题是很难的。于是，你只需要能够以比较高的概率回答的比较准确即可。

【评分方法】

对于每个输出文件，如果其中你有95%以上的行的答案和正确答案的误差不超过25%，那么你就可以得到分数。所谓误差不超过25%，即，如果正确答案是X，那么你的答案在[0.8X,1.25X]这个闭区间内。

Solution

这题用随机化
我们先把初始的元素每个给一个随机值，
对于每个操作，我们只记录随机出的数值前T小的数(T一般是30)，
那么只要归并排序即可，
而对于每个的答案则是T∗(RAND_MAX)/(max{a1...T})

Code

#include <iostream>#include <cstdio>#include <cstdlib>#include <ctime>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fod(i,a,b) for(int i=a;i>=b;i--)#define min(q,w) ((q)>(w)?(w):(q))#define random(m) (rand()%m*rand()%m)#define max(q,w) ((q)<(w)?(w):(q))using namespace std;const int N=1100500,M=2e9,T=35,CS=8;int read(int &n){    char ch=' ';int q=0,w=1;    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());    if(ch=='-')w=-1,ch=getchar();    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;}int m,n;int sc[N][2];int a[N],b[N][T+1];double ANS[N];int c[T*2];void solve(){    int q,w,ans,mx=0;    fo(i,1,n)mx=max(mx,a[i]=rand()+1);    fo(i,1,n)b[i][1]=i,b[i][0]=1;    fo(I,n+1,m+n)    {        q=sc[I-n][0];w=sc[I-n][1];        int i=1,j=1;        c[0]=0;        while(i<=b[q][0]&&j<=b[w][0])        {            if(a[b[q][i]]<a[b[w][j]])c[++c[0]]=b[q][i],i++;            else if(a[b[q][i]]>a[b[w][j]])c[++c[0]]=b[w][j],j++;            else c[++c[0]]=b[w][j],j++,i++;            if(c[0]>T+10)break;        }        for(;i<=b[q][0]&&c[0]<=T;i++)c[++c[0]]=b[q][i];        for(;j<=b[w][0]&&c[0]<=T;j++)c[++c[0]]=b[w][j];        c[0]=min(c[0],T);        b[I][0]=min(c[0],T);        fo(i,1,b[I][0])b[I][i]=c[i];        if(c[0]>=T)ANS[I-n]+=c[0]*1.0*mx/a[c[c[0]]];        else ANS[I-n]+=c[0];    }}int main(){    srand(time(0));    int q,w;    read(m),read(n);    m=m-n;    fo(i,1,m)read(sc[i][0]),read(sc[i][1]);    fo(i,1,CS)solve();    fo(i,1,m)printf("%d\n",(int)(ANS[i]/CS));    return 0;}