题意:有\(n\)个正整数,求随机选取一个3组合,能构成三角形的概率。

Solution: 很容易想到构造权值序列,对其卷积得到任取两条边(可重复)总长度为某数时的方案数序列,我们希望将它转化为两条边不可重复,并去掉顺序。不妨设给定的 \(N\) 个正整数的集合为 \(S\),卷积后的权值序列为 \(\{c_i\}\) ,那么我们对每一个 \(x \in S\), 对 \(c_x\) 减去 \(1\) 即可。

不妨设选出的组合为 \((i,j,k)\),假设 \({x_i}\) 为已经排序的长度序列,那么我们不妨枚举最长边下标 \(k\),那么这一次的贡献为

\[\sum_{k \in S}(\sum_{i=k+1}^{\infty}{c_i})
\]

前缀和扫一遍即可。当然问题并没有这么简单,我们这样获得的答案无法满足下列约束条件

\[k \geq i, k \geq j, i \neq k, j \neq k
\]

如果所有数字都不同,那么去重将是非常容易的。我们对排序后序列枚举元素下标 \(i \in [0,n)\),那么很显然以 \(i\) 为最大边时非法的共有三类。不妨设\(a_i<a_j\)

满足\(a_i>a_k, a_j>a_k\)的共有

\[\frac{(n-i-1)(n-i-2)}{2}
\]

满足\(a_i<a_k, a_j>a_k\)的共有

\[(n-i-1)i
\]

满足\(a_i=a_k or a_j=a_k\)的共有

\[n-1
\]

于是我们暴力统计一遍就可以得到非法数量,进而得出答案。

原问题是可以出现相同大小的数字的,但容易发现它的去重情况与全部不相同并没有什么差异。因而直接按照上述做法即可得出答案。

我们也可以对这个转化的正确性稍作证明。我们将原始序列中的每一个元素 \(x_i\) 映射到一个新元素

\[y_i = x_i - 0.2 + 0.1 \cdot \frac{i}{n}
\]

容易证明\(x_i,x_j,x_k\)可以构成三角形当且仅当\(y_i,y_j,y_k\)可以构成三角形。因此我们通过构造说明了这个转化的正确性。

笔者菜甚,尝试用前缀和直接求解带重复情况失败,留坑在此。

const int N = 200005;
int T,n,m,x[N],a[N],s[N];
double c[N];
double t[N]; int main() {
ios::sync_with_stdio(false);
cin>>T;
while(T--) {
cin>>n;
memset(a,0,sizeof a);
memset(x,0,sizeof x);
memset(s,0,sizeof s);
memset(c,0,sizeof c);
memset(t,0,sizeof t);
for(int i=0;i<n;i++) cin>>x[i];
for(int i=0;i<n;i++) m=max(x[i],m);
for(int i=0;i<n;i++) a[x[i]]++;
poly p,q;
p.c.resize(m+1);
for(int i=0;i<=m;i++) p.c[i]=a[i];
q=p;
double ans = 0;
p*=q;
for(int i=0;i<=2*m;i++) c[i]=(p.c[i]);
for(int i=0;i<n;i++) c[x[i]*2]-=1;
for(int i=0;i<=2*m;i++) c[i]/=2.0;
t[0]=c[0];
for(int i=1;i<=2*m;i++) t[i]=t[i-1]+c[i];
for(int i=0;i<n;i++) ans+=(t[2*m]-t[x[i]]);
sort(x,x+n);
for(int i=0;i<n;i++) {
ans -= (long long)(n-i-1)*(long long)(n-i-2)/2 + (long long)(n-i-1)*(long long)i + (n-1);
}
printf("%.7f\n",6.0*(double)ans / ((double)n*(n-1.0)*(n-2.0)));
}
}

[HDU4609] 3-idiots - 多项式乘法,FFT的更多相关文章

  1. 多项式乘法(FFT)学习笔记

    ------------------------------------------本文只探讨多项式乘法(FFT)在信息学中的应用如有错误或不明欢迎指出或提问,在此不胜感激 多项式 1.系数表示法  ...

  2. 【learning】多项式乘法&fft

    [吐槽] 以前一直觉得这个东西十分高端完全不会qwq 但是向lyy.yxq.yww.dtz等dalao们学习之后发现这个东西的代码实现其实极其简洁 于是趁着还没有忘记赶紧来写一篇博 (说起来这篇东西的 ...

  3. 洛谷.3803.[模板]多项式乘法(FFT)

    题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...

  4. @总结 - 1@ 多项式乘法 —— FFT

    目录 @0 - 参考资料@ @1 - 一些概念@ @2 - 傅里叶正变换@ @3 - 傅里叶逆变换@ @4 - 迭代实现 FFT@ @5 - 参考代码实现@ @6 - 快速数论变换 NTT@ @7 - ...

  5. [uoj#34] [洛谷P3803] 多项式乘法(FFT)

    新技能--FFT. 可在 \(O(nlogn)\) 时间内完成多项式在系数表达与点值表达之间的转换. 其中最关键的一点便为单位复数根,有神奇的折半性质. 多项式乘法(即为卷积)的常见形式: \[ C_ ...

  6. UOJ 34 多项式乘法 FFT 模板

    这是一道模板题. 给你两个多项式,请输出乘起来后的多项式. 输入格式 第一行两个整数 nn 和 mm,分别表示两个多项式的次数. 第二行 n+1n+1 个整数,表示第一个多项式的 00 到 nn 次项 ...

  7. [HNOI2017] 礼物 - 多项式乘法FFT

    题意:给定两个 \(n\) 元环,环上每个点有权值,分别为 \(x_i, y_i\).定义两个环的差值为 \[\sum_{i=0}^{n-1}{(x_i-y_i)^2}\] 可以旋转其中的一个环,或者 ...

  8. 【Luogu3808】多项式乘法FFT(FFT)

    题目戳我 一道模板题 自己尝试证明了大部分... 剩下的还是没太证出来... 所以就是一个模板放在这里 以后再来补东西吧.... #include<iostream> #include&l ...

  9. 【模板】多项式乘法(FFT)

    题目描述 给定一个n次多项式F(x),和一个m次多项式G(x). 请求出F(x)和G(x)的卷积. 输入输出格式 输入格式: 第一行2个正整数n,m. 接下来一行n+1个数字,从低到高表示F(x)的系 ...

随机推荐

  1. springboot打成war包并携带第三方jar包

    1.修改打包方式为war     <packaging>war</packaging> 2.添加第三方依赖的jar到pom     我的第三方jar包在resoueces目录下 ...

  2. 剑指offer-拓展训练-字符的所有组合-全组合

    /* 题目: 给定不含重复字符字符串的全组合. */ /* 思路: 递归法. 例给定abc,输出的组合长度为1,2,3. 对于长度为2的组合,分选择a(ab,ac)和不选择a的情况(bc). 选择a, ...

  3. git commond 详解

    Git commit git commit 主要是将暂存区里的改动给提交到本地的版本库.每次使用git commit 命令我们都会在本地版本库生成一个40位的哈希值,这个哈希值也叫commit-id, ...

  4. [HNOI2016]网络 [树链剖分,可删除堆]

    考虑在 |不在| 这条链上的所有点上放上一个 \(x\),删除也是,然后用可删除堆就随便草掉了. // powered by c++11 // by Isaunoya #pragma GCC opti ...

  5. centos7虚拟机分配静态IP但是得不到IP、不能上网一种可能的原因和解决办法

    1.首先通过ifconfig查看网卡,发现网卡名称为ens33 2. 在/etc/sysconfig/network-scripts/目录下查看网络配置文件 3. 发现有ifcfg-eth0的配置文件 ...

  6. BFS和队列

    深度优先搜索(DFS)和广度优先搜索(BFS)是基本的暴力技术,常用于解决图.树的遍历问题. 首先考虑算法思路.以老鼠走迷宫为例: (1):一只老鼠走迷宫.它在每个路口都选择先走右边,直到碰壁无法继续 ...

  7. 《操作系统真象还原》ELF文件

    下面是第五章部分内容的收获. 用C语言编写内核 一直以来我们都是用汇编语言编写程序的,但接下来我们或许很少用汇编语言编写代码了,大多数都是使用C语言.为什么要这样呢?书上的解释我看的不是很懂,只能结合 ...

  8. CentOS 7 部署 Redis(单机版)

    一.部署环境说明 软件 版本 安装包 CentOS 7.2 CentOS 7.2 Redis 4.0.14 redis-4.0.14.tar.gz 二.开始部署 安装gcc依赖 [root@bmsof ...

  9. java中数据类型转换注意事项

    1.byte.short.char这三种类型互相做数学运算时都会先提升为int类型后再做运算 char a = 'A'; short b = 1; int num = a + b;//a和b在做运算前 ...

  10. 安装DHCP到CentOS(YUM)

    运行环境 系统版本:CentOS Linux release 7.3.1611 (Core) 软件版本:DHCP-x 硬件要求:无 安装过程 1.安装YUM源,由EPEL提供 [root@localh ...