唯一分解定理+容斥-HDU4497
来源:互联网 发布:淘宝客用什么软件 编辑:程序博客网 时间:2024/06/08 19:27
https://vj.xtuacm.cf/contest/view.action?cid=59#problem/L
这个题目的关键在于map的用法,通过map的灵活运用减少代码量和思维难度
题意:
求同时满足gcd(x,y,z)==g&&lcm(x,y,z)==L的(x,y,z)组合数.
题解:
1.很直观的:若l%g!=0||g>l时,答案为0;
2.若果满足l%g==0,用唯一分解定理把个g,l分解,用map保存g,l相同因子个数的差(关键)
3.对于x,y,z,其对应的质因子pi应该属于[0,s],s为map[pi],pi的个数,而且x,y,z中必须有一个为0,一个为map[pi],这是为了保证lcm=L,gcd=g;
每个pi对应的种数累乘就是答案了。
#include <iostream>#include <cstdio>#include <cmath>#include <cstring>#include <algorithm>#include <map>#define ll long longusing namespace std;map<int,int>mp;map<int,int>::iterator it;void fen(int n,int d){ for(int i=2;i*i<=n;i++) { if(n%i==0) while(n%i==0) { mp[i]+=d; n=n/i; } } if(n>1) mp[n]+=d;}int main(){ int t; scanf("%d",&t); while(t--) { int g,l; mp.clear(); scanf("%d%d",&g,&l); if(g>l||l%g) { printf("0\n"); continue; } fen(g,-1); fen(l,1); ll ans=1; for(it=mp.begin();it!=mp.end();it++) { if(it->second) { int x=it->second; //ans*=((x+1)*(x+1)*(x+1)-2*x*x*x+(x-1)*(x-1)*(x-1)); //容斥原理 ans=ans*6*x;//为什么是6*x在后面说明 } } cout<<ans<<endl; } return 0;}
1.qi的值不等于ei或者wi,那么qi是(ei,wi),一共有wi-ei-1种选择,而且题目中说了相同数字顺序不同算不同情况,所以对于三个互不相同的数字,一共有A(3,3)中情况
也就是说这里ans=6*(wi-ei-1)
2.qi的值等于ei或者wi,那么qi一共有两种选择,且对于给定的三个数字里有两个相同数字的排列方式,一共有A(3,1)种
也就是说这里ans=3*2
两种情况加起来ans=6*(wi-ei) 再根据乘法原理,把所有质因数的6*(wi-ei)相乘即可得到结果
阅读全文
0 0
- 唯一分解定理+容斥-HDU4497
- count_prime----容斥+唯一分解定理
- poj2773 —— 二分 + 容斥原理 + 唯一分解定理
- hdu4497(唯一分解+排列组合算组数)
- poj3904 Sky Code —— 唯一分解定理 + 容斥原理 + 组合
- 唯一分解定理
- 唯一分解定理
- 数论-唯一分解定理
- 唯一分解定理
- [数论]唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 唯一分解定理
- 整数唯一分解定理
- 对四元数解算姿态的理解
- spring boot的ActiveMQ使用
- FZU
- 关于创建对象数组的引用数据类型空指针异常问题
- 【Java集合】HashMap和Hashtable
- 唯一分解定理+容斥-HDU4497
- HDU6069 Counting Divisors【数学】
- hdu 6069 类筛法
- explicit关键字
- XML配置文件头部理解
- all-alert
- Android自动化测试探索
- 2017.08.03
- C++基础之面向对象(二)