muntasir007

Untitled

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