Advertisement
muntasir007

Untitled

Feb 20th, 2016
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.91 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define MAX 10000000
  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(10000000);
  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. for (i=0; primes[i]*primes[i] <= b; i++) {
  50. if (b % primes[i] == 0) {
  51. while (b%primes[i] == 0) {
  52. ber.push_back(primes[i]);
  53. b/=primes[i];
  54. }
  55. }
  56. }
  57. if (b>=2) ber.push_back(b);
  58. for (i=0; primes[i]*primes[i] <= a; i++) {
  59. if (a % primes[i] == 0) {
  60. while (a%primes[i] == 0) {
  61. aer.push_back(primes[i]);
  62. a/=primes[i];
  63. }
  64. }
  65. }
  66. if (a>=2) aer.push_back(a);
  67. long long int gcd = 0, cnt=0, j =0, k = 0, xx = aer[0];
  68. i=0;
  69. while (aer[i] == xx && i<aer.size()) {
  70. i++;
  71. gcd++;
  72. }
  73. for (; i<aer.size(); i++) {
  74. xx = aer[i];
  75. cnt =0;
  76. while (aer[i] == xx && i<aer.size()) {
  77. i++;
  78. cnt++;
  79. }
  80. gcd = GCD(cnt, gcd);
  81. }
  82.  
  83. for (i=0; primes[i]*primes[i] <= gcd; i++) {
  84. if (gcd % primes[i] == 0) {
  85. while (gcd%primes[i] == 0) {
  86. ber.push_back(primes[i]);
  87. gcd/=primes[i];
  88. }
  89. }
  90. }
  91. if (gcd>=2) ber.push_back(gcd);
  92. long long int ans = 1;
  93.  
  94. sort(ber.begin(), ber.end());
  95.  
  96. for (i=0; i<ber.size(); ) {
  97. xx = ber[i];
  98. cnt =0;
  99. while (ber[i] == xx && i<ber.size()) {
  100. i++;
  101. cnt++;
  102. }
  103. ans*= (cnt+1);
  104. }
  105.  
  106. printf("Case %d: %lld\n",cs, ans-1);
  107. }
  108. return 0;
  109. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement