题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4609

题意:给出n个正整数(数组A)。每次随机选出三个数。问这三个数能组成三角形的概率为多大?

思路:求出有多少种选择的方案,除以总选择方案即可。用num[i]表示长度为i的出现几次。

对于样例1 3 3 4,我们得到num={0,1,0,2,1},

对num求卷积,得到:num={0,0,1,0,4,2,4,4,1}。此时的num[i]表示选择两个数和为i的选择方案的种数。 但是这里有重复的:

(1)两次不能选同一个数,这种情况要删去;

(2)对于两个数x和y,先选x后选y和先选y后选x,所以多算了一倍。

减去上面两种重复之后,我们得到:num={0,0,0,0,2,1,1,2}。接着设sum[i]为num的前缀和。接着将给出的数组A排序,开始统计答案。对于A[i]我们认为三个数中它是最大的,因此我们要选出在i之前的两个数使得其和大于A[i],首先大于A[i]的情况有sum[Maxlen]-sum[A[i]],然后要减去重的:

(1)一个选了i后面的一个选了i前面的;

(2)一个选了i另一个随便;

(3)两个都是i后面的。

struct node
{
    double x,y;
    
    node(double _x=0.0,double _y=0.0)
    {
        x=_x;
        y=_y;
    }
    
    node operator+(node a)
    {
        return node(x+a.x,y+a.y);
    }
    
    node operator-(node a)
    {
        return node(x-a.x,y-a.y);
    }
    
    node operator*(node a)
    {
        return node(x*a.x-y*a.y,x*a.y+y*a.x);
    }
};

node A[N];
int L;

int reverse(int x)
{
    int ans=0,i;
    FOR0(i,L) if(x&(1<<i)) ans|=1<<(L-1-i);
    return ans;
}

void bitReverseCopy(node a[],int n)
{
    int i;
    FOR0(i,n) A[i]=a[i];
    FOR0(i,n) 
    {
        a[reverse(i)]=A[i];
    }
}

void fft(node a[],int n,int on)
{
    bitReverseCopy(a,n);
    int len,i,j,k;
    node x,y,u,t;
    for(len=2;len<=n;len<<=1)
    {
        x=node(cos(-on*2*PI/len),sin(-on*2*PI/len));
        for(j=0;j<n;j+=len)
        {
            y=node(1,0);
            for(k=j;k<j+len/2;k++)
            {
                u=a[k];
                t=y*a[k+len/2];
                a[k]=u+t;
                a[k+len/2]=u-t;
                y=y*x;
            }
        }
    }
    if(on==-1)
    {
        FOR0(i,n) a[i].x/=n;
    }
}

node a[N];
int d[N],n,m;
i64 num[N],sum[N];

int main()
{
    rush()
    {
        RD(n);
        clr(num,0);
        int i;
        FOR0(i,n) RD(d[i]),num[d[i]]++;
        sort(d,d+n);
        int len=d[n-1]+1;
        m=1; L=0;
        while(m<2*len) m<<=1,L++;
        FOR0(i,len) a[i]=node(num[i],0);
        for(i=len;i<m;i++) a[i]=node(0,0);
        fft(a,m,1);
        FOR0(i,m) a[i]=a[i]*a[i];
        fft(a,m,-1);
        FOR0(i,m) num[i]=(i64)(a[i].x+0.5);
        m=d[n-1]<<1;
        FOR0(i,n) num[d[i]<<1]--;
        FOR1(i,m) num[i]>>=1,sum[i]=sum[i-1]+num[i];
        i64 ans=0;
        FOR0(i,n) 
        {
            ans+=sum[m]-sum[d[i]];
            ans-=(i64)i*(n-i-1);
            ans-=n-1;
            ans-=(i64)(n-1-i)*(n-2-i)/2;
        }
        i64 Sum=(i64)n*(n-1)*(n-2)/6;
        PR(1.0*ans/Sum);
     }
}

HDU 4609 3-idiots(FFT)的更多相关文章

  1. HDU 5763 Another Meaning(FFT)

    [题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5763 [题目大意] 给出两个串S和T,可以将S串中出现的T替换为*,问S串有几种表达方式. [题解 ...

  2. 2019.01.02 bzoj3513: [MUTC2013]idiots(fft)

    传送门 fftfftfft经典题. 题意简述:给定nnn个长度分别为aia_iai​的木棒,问随机选择3个木棒能够拼成三角形的概率. 思路:考虑对于木棒构造出生成函数然后可以fftfftfft出两个木 ...

  3. A * B Problem Plus HDU - 1402 (FFT)

    A * B Problem Plus HDU - 1402 (FFT) Calculate A * B.  InputEach line will contain two integers A and ...

  4. 快速傅里叶(FFT)的快速深度思考

    关于按时间抽取快速傅里叶(FFT)的快速理论深度思考 对于FFT基本理论参考维基百科或百度百科. 首先谈谈FFT的快速何来?大家都知道FFT是对DFT的改进变换而来,那么它究竟怎样改进,它改进的思想在 ...

  5. HDU 5938 Four Operations(四则运算)

    p.MsoNormal { margin: 0pt; margin-bottom: .0001pt; text-align: justify; font-family: Calibri; font-s ...

  6. HDU 5775 Bubble Sort(冒泡排序)

    p.MsoNormal { margin: 0pt; margin-bottom: .0001pt; text-align: justify; font-family: Calibri; font-s ...

  7. HDU 1711 Number Sequence(数列)

    HDU 1711 Number Sequence(数列) Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Ja ...

  8. HDU 1005 Number Sequence(数列)

    HDU 1005 Number Sequence(数列) Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Jav ...

  9. 【BZOJ3527】力(FFT)

    [BZOJ3527]力(FFT) 题面 Description 给出n个数qi,给出Fj的定义如下: \[Fj=\sum_{i<j}\frac{q_i q_j}{(i-j)^2 }-\sum_{ ...

  10. 【BZOJ4827】【HNOI2017】礼物(FFT)

    [BZOJ4827][HNOI2017]礼物(FFT) 题面 Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一 个送给她.每 ...

随机推荐

  1. bzoj 2751 快速幂

    首先我们知道,对于所有种情况,我们可以将每一位可以放的 数的值加起来,所有位置的乘起来,等于的就是最后的答案,具体 为什么正确,可以根据乘法分配律来想一想. 那么对于所有不做要求的,快速幂直接算就行了 ...

  2. 【Dancing Link专题】解题报告

    DLX用于优化精确覆盖问题,由于普通的DFS暴力搜索会超时,DLX是一个很强有力的优化手段,其实DLX的原理很简单,就是利用十字链表的快速删除和恢复特点,在DFS时删除一些行和列以减小查找规模,使得搜 ...

  3. [工作积累] Google Play Game SDK details

    https://developers.google.com/games/services/cpp/api/structgpg_1_1AndroidSupport For apps which targ ...

  4. IE8下兼容rgba颜色的半透明背景

    在工作中做一个图片半透明遮罩时发现在IE8下不兼容 一查再知道IE8不支持rgba颜色,再搜搜兼容性方法,没想到这么快就解决了. 先说说rgba的含义: r代表red,g代表green,b代表blue ...

  5. 一道PK赛题

    Problem Description I think that you might have played the traditional Chinese ring game: The Chines ...

  6. HDU-2604_Queuing

    题目:Problem Description Queues and Priority Queues are data structures which are known to most comput ...

  7. $q -- AngularJS中的服务

    此 承诺/延迟(promise/deferred)实现 的灵感来自于 Kris Kowal's QCommonJS Promise建议文档 将承诺(promise) 作为和 异步执行操作(action ...

  8. (转)OpenCV 2.4.8 +VS2010的开发环境配置

    转自:  http://blog.csdn.net/poem_qianmo/article/details/19809337 自己可能需要再进行修改 本系列文章由zhmxy555(毛星云)编写,转载请 ...

  9. excel设置下拉菜单,并且不同值会显示不同颜色

    工作中常常要用的excel,每次都会有新的需求,然后不会,然后百度,然后过段时间可能就又忘了,于是就想说,自己记录下来~~~因为自己用的都是2010,其实哪个版本都差不多,都是应该可以找到相应的按钮滴 ...

  10. MongoDB (九) MongoDB 投影

    mongodb 投影意思是只选择必要的数据而不是选择一个文件的数据的整个.如果一个文档有5个字段,需要显示只有3个,然后选择其中只有3个字段. find() 方法 MongoDB 的find()方法, ...