HDU 4609——3-idiots
题意:
给定n个数,随机从这n个数中取3个数,问能组成三角形的概率是多少?
思路:
首先把统计这n个数出现的个数,那么会得到一个向量,这个向量的自我的乘积就是a[i]+a[j]的可能的方案数,这样,我么就很方便求出了两条边的和的方案数。
但是在加的过程中,我么多加了a[i]+a[i]的情况,也就是自己跟自己相加的情况,所以我们要减去这一段。而且,在加的过程a[i]+a[j]和a[j]+a[i]是同一种情况,我们加了两次,所以要减去一半。
sum是num的前缀和,假设a[i]是三角形内要取得最长的边,那么只需要在前面取两条边并且这两条边的和大于a[i];
而长度和大于a[i]的取两个的取法是sum[len]-sum[a[i]].
但是这里面有不符合的。
一个是包含了取一大一小的
cnt -= (long long)(n-1-i)*i;
一个是包含了取一个本身i,然后取其它的
cnt -= (n-1);
还有就是取两个都大于的了
cnt -= (long long)(n-1-i)*(n-i-2)/2;
这样把i从0~n-1累加,就答案了。
参考了kuangbin 的博客:
http://www.cnblogs.com/kuangbin/archive/2013/07/24/3210565.html
code:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;const double PI=acos(-1.0);
typedef long long ll;struct complex
{double l,r;complex(double ll=0.0,double rr=0.0){l=ll;r=rr;}complex operator +(const complex& B){return complex(l+B.l,r+B.r);}complex operator - (const complex& B){return complex(l-B.l,r-B.r);}complex operator *(const complex& B){return complex(l*B.l-r*B.r,l*B.r+B.l*r);}
};/** 进行FFT和IFFT前的反转变换。* 位置i和j(i二进制反转后位置)互换* len必须是2的幂*/
void change(complex y[],int len){int i,j,k;for (int i=1,j=len/2;i<len-1;i++){if (i<j) swap(y[i],y[j]);k=len/2;while (j>=k){j-=k;k>>=1;}if (j<k) j+=k;}
}
/** 做FFT* len必须为2^k形式,* on==1时是DFT,on==-1时是IDFT*/
void fft(complex y[],int len,int on){change(y,len);for (int h=2;h<=len;h<<=1){complex wn(cos(-on*2*PI/h),sin(-on*2*PI/h));for (int j=0;j<len;j+=h){complex w(1,0);for (int k=j;k<j+h/2;k++){complex u=y[k];complex t=w*y[k+h/2];y[k]=u+t;y[k+h/2]=u-t;w=w*wn;}}}if (on==-1){for (int i=0;i<len;i++){y[i].l/=len;}}
}
const int MAXN=400040;
complex x1[MAXN];
ll sum[MAXN],num[MAXN];
int a[MAXN/4];int main()
{int T,n;scanf("%d",&T);while (T--){scanf("%d",&n);memset (num,0,sizeof(num));for (int i=0;i<n;i++){scanf("%d",a+i);num[a[i]]++;}sort(a,a+n);int len1=a[n-1]+1,len=1;while (len<2*len1) len<<=1;for (int i=0;i<len1;i++){x1[i]=complex(num[i],0);}for (int i=len1;i<len;i++){x1[i]=complex(0,0);}fft(x1,len,1);for (int i=0;i<len;i++) x1[i]=x1[i]*x1[i];fft(x1,len,-1);for (int i=0;i<len;i++) {num[i]=(ll) (x1[i].l+0.5);}len=2*a[n-1];for (int i=0;i<n;i++) num[a[i]+a[i]]--;for (int i=1;i<=len;i++) num[i]/=2;sum[0]=0;for (int i=1;i<=len;i++) sum[i]=sum[i-1]+num[i];ll cnt=0;for (int i=0;i<n;i++){cnt+=sum[len]-sum[a[i]];cnt-=(ll)(n-i-1)*i;cnt-=(n-1);cnt-=(ll) (n-1-i)*(n-i-2)/2;}ll tot=(ll)n*(n-1)*(n-2)/6;printf("%.7f\n",1.0*cnt/tot);}
}
HDU 4609——3-idiots相关推荐
- HDU 4609 3-idiots
HDU - 4609 思路: 记录每个木棍长度出现的次数,然后就可以用用类似多项式的乘法(专业术语:卷积,因为是下标和为一特定值的积的和(x+y=k),相当于在笛卡尔坐标系中将这条直线卷起来,故得名卷 ...
- 快速傅里叶变换应用之二 hdu 4609 3-idiots
快速傅里叶变化有不同的应用场景,hdu4609就比较有意思.题目要求是给n个线段,随机从中选取三个,组成三角形的概率. 初始实在没发现这个怎么和FFT联系起来,后来看了下别人的题解才突然想起来:组合计 ...
- hdu 4609 3-idiots(FFT计数)
题目链接:hdu 4609 3-idiots 题意: 给你n条线段.问随机取三个,可以组成三角形的概率. 题解: FFT搞,具体可看kuangbin菊苣的详细题解:传送门 1 #include< ...
- HDU 4609 3-idiots FFT
3-idiots 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=4609 Description King OMeGa catched three m ...
- HDU 4609 3idiots
第一道FFT,研究了kuangbin神的版子,发现想不太通,不过基本的意思都懂了,就是用一个数组的权和另一个数组的权求和,计算个数.题目中有不少坑点. #include <cstdio> ...
- HDU 4609 FFT
题目大意 给定n条边的边值,求任意取三条边能组成三角形的概率 这里概率 P = valid/tot tot = (n-2)*(n-1)*n/6是没问题的 valid表示合法的方式 先考虑,任意两条边组 ...
- HDU 4609 3-idiots(FFT)
3-idiots 思路 多项式卷积的经典应用了: 看样例一: 4 1 3 3 4 我们用多项式系数表示得到一个与aia_iai有关的多项式:010210\ 1\ 0\ 2\ 10 1 0 2 1, ...
- 3-idiots HDU - 4609 FFT模板
给N个木棍问任选三个可组合出三角形的概率 以前写的,整理整理存个档 //#pragma comment(linker, "/STACK:1024000000,1024000000" ...
- [kuangbin]各种各样的题单
[kuangbin]各种各样的题单 专题1 简单搜索 POJ 1321 POJ 2251 POJ 3278 POJ 3279 POJ 1426 POJ 3126 POJ 3087 POJ 3414 F ...
最新文章
- JSP大作业数据库_本地MySQL【种种问题】
- angularjs loading, animate
- OPNET仿真zigbee
- 来吧!我教你画真正的流程图
- php 请求url 返回json,php解析请求url并返回json数据
- java 图片刷新页面_js修改img的src属性刷新图片时的图片缓存问题
- 浅谈一个缓存小实践,另外还有更多案例,欢迎关注~
- 深入理解计算机系统:计算机系统概述
- IntelliJ IDEA for Mac中Maven菜单用法详解
- 《华为交换机学习指南》学习笔记·二
- 下周把交接文档做完就交接实习结束离职了
- AI(adobe illustrator)怎么设置导出图片的像素尺寸
- SAMSUNG S3C2440的简易BootLoader ㈢
- Chrome插件:UserAgent Switcher 浏览器UA切换插件 平台与浏览器可独立选择
- 信息系统综合测试与管理
- MATLAB在线文档打不开,怎么解决
- hasLayout引发的CSS Bug表 及处理办法【zoom:1】
- 高斯滤波与双边滤波原理、MATLAB实现及结果对照
- 【python辅助excel】(3)
- 什么是MBR/UEFI/GPT?MBR与UEFI具体区别是什么?