muntasir007

Untitled

Feb 20th, 2016
121
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.06 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define MAX 20000000
  4. vector <int> primes;
  5. long long int GCD( long long int a, long long int b ) {
  6. if( b == 0 ) return a;
  7. if( a < b ) swap( a, b );
  8. int r = a % b;
  9. return GCD( b, r );
  10. }
  11. int prime[MAX];
  12. int setBit( int n, int position ) {
  13. n = n | ( 1 << position );
  14. return n;
  15. }
  16. bool checkBit(int n, int position) {
  17. return n & ( 1 << position );
  18. }
  19. void primeGenerator( int n ) {
  20. int x = sqrt( n );
  21. primes.push_back(2);
  22. prime[0] = setBit( prime[0], 0 );
  23. prime[0] = setBit( prime[0], 1 );
  24. for( int i = 4; i <= n; i += 2 )
  25. prime[ i >> 5 ] = setBit( prime[ i >> 5 ], i & 31 );
  26. int i;
  27. for( i = 3; i <= x; i += 2 ) {
  28. if( !checkBit( prime[ i >> 5 ], i & 31 ) ) {
  29. primes.push_back(i);
  30. for( int j = i+i; j <= n; j += i )
  31. prime[ j >> 5 ] = setBit( prime[ j >> 5 ], j & 31 );
  32. }
  33. }
  34. for (; i<20000000; i+=2) {
  35. if( !checkBit( prime[ i >> 5 ], i & 31 ) ) {
  36. primes.push_back(i);
  37. }
  38. }
  39. return;
  40. }
  41. int main() {
  42. primeGenerator(20000000);
  43.  
  44. int T;
  45.  
  46. scanf("%d", &T);
  47. int cs;
  48.  
  49. for (cs=1; cs<=T; cs++) {
  50. vector <int> ber;
  51. vector <int> aer;
  52. long long int a, b, bb;
  53. scanf("%lld %lld", &a, &b);
  54. bb = b;
  55. long long int i;
  56.  
  57. for (i=0; primes[i]*primes[i] <= b; i++) {
  58. if (b % primes[i] == 0) {
  59. while (b%primes[i] == 0) {
  60. ber.push_back(primes[i]);
  61. b/=primes[i];
  62. }
  63. }
  64. }
  65. if (b>=2) ber.push_back(b);
  66.  
  67. for (i=0; primes[i]*primes[i] <= a; i++) {
  68. if (a % primes[i] == 0) {
  69. while (a%primes[i] == 0) {
  70. aer.push_back(primes[i]);
  71. a/=primes[i];
  72. }
  73. }
  74. }
  75. if (a>=2) aer.push_back(a);
  76. long long int gcd = 0, cnt=0, j =0, k = 0, xx = aer[0];
  77. i=0;
  78. while (aer[i] == xx && i<aer.size()) {
  79. i++;
  80. gcd++;
  81. }
  82.  
  83. for (; i<aer.size(); i++) {
  84. xx = aer[i];
  85. cnt =0;
  86. while (aer[i] == xx && i<aer.size()) {
  87. i++;
  88. cnt++;
  89. }
  90. gcd = GCD(cnt, gcd);
  91. }
  92.  
  93. for (i=0; primes[i]*primes[i] <= gcd; i++) {
  94. if (gcd % primes[i] == 0) {
  95. while (gcd%primes[i] == 0) {
  96. ber.push_back(primes[i]);
  97. gcd/=primes[i];
  98. }
  99. }
  100. }
  101. if (gcd>=2) ber.push_back(gcd);
  102. long long int ans = 1;
  103.  
  104. sort(ber.begin(), ber.end());
  105.  
  106. for (i=0; i<ber.size(); ) {
  107. xx = ber[i];
  108. cnt =0;
  109. while (ber[i] == xx && i<ber.size()) {
  110. i++;
  111. cnt++;
  112. }
  113. ans*= (cnt+1);
  114. }
  115.  
  116. printf("Case %d: %lld\n",cs, ans-1);
  117. }
  118. return 0;
  119. }
Advertisement
Add Comment
Please, Sign In to add comment