POJ 3904 Sky Code (容斥原理)
来源:互联网 发布:藏族歌曲下载软件 编辑:程序博客网 时间:2024/05/27 19:25
B - Sky Code
Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64uDescription
Stancu likes space travels but he is a poor software developer and will never be able to buy his own spacecraft. That is why he is preparing to steal the spacecraft of Petru. There is only one problem – Petru has locked the spacecraft with a sophisticated cryptosystem based on the ID numbers of the stars from the Milky Way Galaxy. For breaking the system Stancu has to check each subset of four stars such that the only common divisor of their numbers is 1. Nasty, isn’t it? Fortunately, Stancu has succeeded to limit the number of the interesting stars to N but, any way, the possible subsets of four stars can be too many. Help him to find their number and to decide if there is a chance to break the system.
Input
In the input file several test cases are given. For each test case on the first line the number N of interesting stars is given (1 ≤ N ≤ 10000). The second line of the test case contains the list of ID numbers of the interesting stars, separated by spaces. Each ID is a positive integer which is no greater than 10000. The input data terminate with the end of file.
Output
For each test case the program should print one line with the number of subsets with the asked property.
Sample Input
42 3 4 5 42 4 6 8 72 3 4 5 7 6 8
Sample Output
1 0 34
题意:给定n个数,从n个数找出四个数,使这四个数的最大公约数为1,找出有多少对这样的组合。
找最大公约数不是1的有多少对。
题解:四个数的公约数为1,并不代表四个数两两互质。比如(2,3,4,5)公约数为1,但是
2和4并不互质。从反面考虑,先求出四个数公约数不为1的情况个数,用总的方案个数
减去四个数公约数不为1的情况个数就是所求。
求四个数公约数不为1的情况个数,需要将N个数每个数质因数分解,纪录下所有不同
的素因子所能组成的因子(就是4个数的公约数),并统计构成每种因子的素因子个数,
和因子总数。然后再计算组合数。比如说因子2的个数为a,则四个数公约数为2的个数
为C(a,4),因子3的个数为b,则四个数公约数为3的个数为C(b,4),因子6(2*3)的个
数为c,则四个数公约数的个数为C(c,4)。
但是公约数为2的情况中或者公约数为3的情况中可能包括公约数为6的情况,相当于几
个集合求并集,这就需要容斥定理来做。
容斥原理应用,以2为因子的数有a个,3为因子 的数有b个,6为因子的数有c个,
n个数不互质的四元组个数为C(a,4)+C(b,4)-C(c,4) (含奇数个素因子的加,偶数个素因子的减),
下面就是统计出2,3,5这些因子的倍数的个数,对C(a,4)容斥!
代码:弄清思路以后就很好做了,一环扣一环,用二进制进行枚举,很棒!
#include <iostream>#include <stdio.h>#include <algorithm>#include <string.h>#include <math.h>#include <stdlib.h>using namespace std;typedef long long ll;ll c(ll n){ return n*(n-1)*(n-2)*(n-3)/24;}ll f[10010],cnt,ccount[10010],num[10010];void div(ll n){ cnt=0; for(int i=2;i<=sqrt(n);i++) { if(n%i==0) f[cnt++]=i; while(n%i==0) n/=i; } if(n!=1) f[cnt++]=n;}void solve(ll n){ div(n); for(int i=1;i<(1<<cnt);i++) { ll tmp=1,flag=0; for(int j=0;j<cnt;j++) { if(i&(1<<j)) { flag++; tmp*=f[j]; if(tmp>n) break; } } ccount[tmp]++; //count是关键字 不能用 num[tmp]=flag; }}int main(){ ll n,m; while(scanf("%lld",&n)!=EOF) { memset(ccount,0,sizeof(ccount)); for(int i=0;i<n;i++) { scanf("%lld",&m); solve(m); } ll ans=c(n); for(int i=0;i<=10000;i++) { if(ccount[i]) { if(num[i]&1) ans-=c(ccount[i]); else ans+=c(ccount[i]); } } printf("%lld\n",ans); } return 0;}
0 0
- POJ 3904 Sky Code 容斥原理
- poj 3904 sky code 容斥原理
- POJ 3904 Sky Code 容斥原理
- poj 3904 Sky Code【容斥原理】
- 容斥原理:poj 3904 Sky Code
- Sky Code - POJ 3904 容斥原理
- Sky Code - POJ 3904 容斥原理
- POJ 3904 Sky Code (容斥原理)
- POJ 3904 Sky Code (容斥原理)
- POJ 3904 Sky Code (容斥原理)
- poj 3904 Sky Code(容斥原理)
- [容斥原理] poj 3094 Sky Code
- POJ 3904:Sky Code _容斥原理
- POJ 3904 Sky Code 解题报告(容斥原理)
- POJ 3904 Sky Code (容斥原理)
- POJ 3904 Sky Code(容斥)
- POJ 3904 Sky Code(容斥)
- poj 3904 Sky Code 莫比乌斯反演 或 容斥原理
- 解析Json 出现——java.lang.RuntimeException: java.lang.reflect.InvocationTargetException
- HDU 2841 Visible Trees 数论+容斥原理
- UVA 10325 The Lottery( 容斥原理)
- iOS之一个超赞的视频直播、第三方库,直播看这个就够了,支持RTMP推流,美颜直播
- HDU 4334 Trouble (暴力)
- POJ 3904 Sky Code (容斥原理)
- Kubernetes& swarm对比
- Shell命令
- HDU 4059 容斥原理+快速幂+逆元
- HDU5072 容斥原理
- 8种机器学习算法
- 如何创建方框型的EditText
- CF GYM 100548 容斥原理+快速幂+逆元
- 【JavaScript】正则表达式