[归并][随机算法] JZOJ P3765 想法

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}。

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

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

Input

第一行两个整数N，M，依次表示小强的题目数量和想法的数量

接下来N-M行，每行两个整数，依次表示小强组合出来的题目都是由哪两个题组合而成的。M个想法对应的题目依次编号为1~M。之后，小强组合出来的第一个题编号为M+1，组合出来的第二个题编号为M+2，依次类推。

Output

输出N-M行，每行一个整数表示小强组合出来的每个题都涉及了几个想法。

Sample Input

6 3

1 2

2 3

4 5

Sample Output

2

2

3

Data Constraint

对于30%的数据，M≤1000,N≤10000

对于60%的数据，M≤10000,N≤100000

对于100%的数据，M≤100000,N≤1000000

Hint

【评分方法】

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

题解

STM出题人一上来就说，这题是随机数算法，是一个允许有误差的SPJ

随机数大法好！！！（可惜子豪不在）

这题其实就是我们随机生产m个数，将每个想法用这些值赋值
然后我们再进行归并

那么怎么归并

我们可以每次加入！！可以形成它的思想的序号之前归并的随机数！！
每次从小到大加入，因为之前合并过的都是按升序的，所以可以将后面一串一样的bank掉
这样可以保证每次加进来的就是按升序的没有重复的序列

怎么求值

如果每次加入的数小于30,则直接累计
如果大于30，则加m*3*T/f[i][T]
（虽然我不知到怎么推到的）

引用讲题人一句话：

A：那这公式是怎么推出来的？
B：我也不知道
A：那你怎么A的
B：（手指这公式）这里不是有结论吗/
A：……

最后跑五次，求平均值输出

代码

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,T,mx,x[1000003],y[1000003],f[1000003][33];
double ans[1000003];
void gb(int *a,int *b,int *c)
{
    int i=1,j=1;
    c[0]=0;
    for (int z=1;z<=T&&(i<=a[0]||j<=b[0]);z++)
    {
        c[0]++;
        if ((a[i]<b[j]&&i<=a[0])||(j>b[0])) c[z]=a[i++]; else c[z]=b[j++];
        while (i<=a[0]&&c[z]==a[i]) i++;
        while (j<=b[0]&&c[z]==b[j]) j++;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=m+1;i<=n;i++) scanf("%d%d",&x[i],&y[i]);
    T=30;
    for (int t=1;t<=5;t++)
    {
        mx=3*m;
        for (int i=1;i<=m;i++) f[i][1]=rand()%mx+1,f[i][0]=1;
        for (int i=m+1;i<=n;i++)
        {
            gb(f[x[i]],f[y[i]],f[i]);
            if (f[i][0]<T) ans[i]+=f[i][0]; else ans[i]+=mx*T/f[i][T];
        }
    }
    for (int i=m+1;i<=n;i++) printf("%.0lf\n",ans[i]/5);
    return 0;
}