标签：HDU include return gcd int ll 4497 LCM 31
题目链接：https://vjudge.net/problem/HDU-4497
题意：求有多少组（x,y,z)满足gcd(x,y,z)=a,lcm(x,y,z)=b。
思路：对于x,y,z都可以写成x = p1^a1+p2^a2+p3^a3....pn^an;y = p1^b1+p2^b2+p3^b3....pn^bn;z=p1^c1+p2^c2+p3^c3....pn^cn;
x,y,z的最大公约数可以写成gcd(x,y,z) = p1^min(a1,b1,c1)*p2^min(a2,b2,c2)......pn^min(an,bn,cn);
最小公倍数则是lcm(x,y,z) = p1^max(a1,b1,c1)*p2^max(a2,b2,c2)....*pn^max(an,bn,cn)
对于gcd(x,y,z) = g,lcm(x,y,z) = l的组合数等gcd(x,y,z) = 1,lcm(x,y,z) = l/g;
设l/g = k;
k = p1^d1+p2^d2+.....+pn^dn;
这样对于p1来说a1,b1,c1中一定有一个为0，因为gcd(x,y,z) = 1,同理也一定有一个为d1;
分为三种情况分别是2个0一个d1,2个d1一个0,1个d1一个0另外一个是[1，(d1-1)]中的一个数，则情况数便是3+3+(d1-1)*6 = 6*d1;
最终答案便是
ans = 1;
for(int i=1;i<=n;i++)
ans*=di*6;
1 #include<time.h>
2 #include <set>
3 #include <map>
4 #include <stack>
5 #include <cmath>
6 #include <queue>
7 #include <cstdio>
8 #include <string>
9 #include <vector>
10 #include <cstring>
11 #include <utility>
12 #include <cstring>
13 #include <iostream>
14 #include <algorithm>
15 #include <list>
16 using namespace std;
17 #define eps 1e-10
18 #define PI acos(-1.0)
19 #define lowbit(x) ((x)&(-x))
20 #define zero(x) (((x)>0?(x):-(x))<eps)
21 #define mem(s,n) memset(s,n,sizeof s);
22 #define ios {ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);}
23 typedef long long ll;
24 typedef unsigned long long ull;
25 const int maxn=1e5+5;
26 const int Inf=0x7f7f7f7f;
27 const ll Mod=999911659;
28 //const int N=3e3+5;
29 bool isPowerOfTwo(int n) { return n > 0 && (n & (n - 1)) == 0; }//判断一个数是不是 2 的正整数次幂
30 int modPowerOfTwo(int x, int mod) { return x & (mod - 1); }//对 2 的非负整数次幂取模
31 int getBit(int a, int b) { return (a >> b) & 1; }// 获取 a 的第 b 位，最低位编号为 0
32 int Max(int a, int b) { return b & ((a - b) >> 31) | a & (~(a - b) >> 31); }// 如果 a>=b,(a-b)>>31 为 0，否则为 -1
33 int Min(int a, int b) { return a & ((a - b) >> 31) | b & (~(a - b) >> 31); }
34 ll gcd(ll a, ll b) {return b ? gcd(b, a % b) : a;}
35 ll lcm(ll a, ll b) {return a / gcd(a, b) * b;}
36 int Abs(int n) {
37 return (n ^ (n >> 31)) - (n >> 31);
38 /* n>>31 取得 n 的符号，若 n 为正数，n>>31 等于 0，若 n 为负数，n>>31 等于 -1
39 若 n 为正数 n^0=n, 数不变，若 n 为负数有 n^(-1)
40 需要计算 n 和 -1 的补码，然后进行异或运算，
41 结果 n 变号并且为 n 的绝对值减 1，再减去 -1 就是绝对值 */
42 }
43 ll binpow(ll a, ll b,ll c) {
44 ll res = 1;
45 while (b > 0) {
46 if (b & 1) res = res * a%c;
47 a = a * a%c;
48 b >>= 1;
49 }
50 return res%c;
51 }
52 void extend_gcd(ll a,ll b,ll &x,ll &y)
53 {
54 if(b==0) {
55 x=1,y=0;
56 return;
57 }
58 extend_gcd(b,a%b,x,y);
59 ll tmp=x;
60 x=y;
61 y=tmp-(a/b)*y;
62 }
63 ll mod_inverse(ll a,ll m)
64 {
65 ll x,y;
66 extend_gcd(a,m,x,y);
67 return (m+x%m)%m;
68 }
69 ll eulor(ll x)
70 {
71 ll cnt=x;
72 ll ma=sqrt(x);
73 for(int i=2;i<=ma;i++)
74 {
75 if(x%i==0) cnt=cnt/i*(i-1);
76 while(x%i==0) x/=i;
77 }
78 if(x>1) cnt=cnt/x*(x-1);
79 return cnt;
80 }
81 ll p[maxn],cnt;
82 ll a,b,n;
83 ll num[maxn];
84 void getp(ll x)
85 {
86 cnt=0;
87 p[cnt++]=x;
88 for(ll i=2;i*i<=x;i++)
89 {
90 if(x%i==0)
91 {
92 p[cnt++]=i;
93 if(i!=x/i)
94 {
95 p[cnt++]=x/i;
96 }
97 }
98 }
99 sort(p,p+cnt);
100 return;
101 }
102 ll ans[maxn];
103 int main()
104 {
105 int t;
106 scanf("%d",&t);
107 while(t--)
108 {
109 mem(ans,0);
110 ll a,b;
111 scanf("%lld%lld",&a,&b);
112 if(b%a)
113 {
114 puts("0");
115 continue;
116 }
117 ll sum=1;
118 ll c=b/a;
119 int cnt=0;
120 for(int i=2;i<=c;i++)
121 {
122 if(c==0) break;
123 if(c%i==0)
124 {
125 while(c%i==0)
126 {
127 ans[cnt]++;
128 c/=i;
129 }
130 cnt++;
131 }
132 }
133 for(int i=0;i<cnt;i++)
134 {
135 if(ans[i])
136 {
137 //cout<<ans[i]<<endl;
138 sum*=ans[i]*6;
139 }
140 }
141 cout<<sum<<endl;
142 }
143 return 0;
144 }
View Code
标签：HDU,include,return,gcd,int,ll,4497,LCM,31 来源： https://www.cnblogs.com/zpj61/p/13551132.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。