Advertisement
erfanul007

UVa 10140

Dec 26th, 2019
200
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 3.95 KB | None | 0 0
  1. #include<bits/stdc++.h>
  2. using namespace std;
  3.  
  4. // #pragma GCC optimize("Ofast,no-stack-protector")
  5. // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
  6. // #pragma GCC optimize("unroll-loops")
  7.  
  8. typedef unsigned long long ull;
  9. typedef long long          ll;
  10. typedef long double        lld;
  11. typedef vector<int>        vi;
  12. typedef vector<ll>         vll;
  13. typedef vector<vi>         vvi;
  14. typedef pair<int,int >     pii;
  15. typedef vector< pii >      vpii;
  16. typedef set<int>           sti;
  17.  
  18. #define sc                 scanf
  19. #define pf                 printf
  20. #define sci(n)             scanf("%d",&n)
  21. #define scii(n,m)          scanf("%d %d",&n,&m)
  22. #define scl(n)             scanf("%lld",&n)
  23. #define scd(n)             scanf("%lf",&n)
  24. #define scs(s)             scanf("%s",s)
  25. #define pfi(n)             printf("%d",n)
  26. #define pfl(n)             printf("%lld",n)
  27. #define pff(n)             cout<<n
  28. #define spc                printf(" ")
  29. #define line               printf("\n")
  30. #define loop(i,x,y)        for(int i=int(x); i<=int(y); i++)
  31. #define rloop(i,y,x)       for(int i=int(y); i>=int(x); i--)
  32. #define cspf(i)            printf("Case %d: ", i)
  33. #define vout(v)            for(int w=0;w<v.size();w++){if(w) cout<<' '; cout<<v[w];}
  34. #define clr(a,x)           memset(a,x,sizeof(a))
  35. #define tp(v,j)            get<j>(v)
  36. #define pb                 push_back
  37. #define mp                 make_pair
  38. #define mt                 make_tuple
  39. #define ff                 first
  40. #define ss                 second
  41. #define all(v)             v.begin(),v.end()
  42. #define rall(v)            v.rbegin(),v.rend()
  43. #define read()             freopen("input.txt", "r", stdin)
  44. #define write()            freopen("output.txt", "w", stdout)
  45. #define fastIO()           ios_base::sync_with_stdio(false); cin.tie(NULL);
  46. /// Constants
  47. #define eps                1e-6
  48. #define PI                 acos(-1.0)  // 3.1415926535897932
  49. #define MAX                (ll)2147483647
  50. #define MAXn               (ll)50000
  51.  
  52. int GCD(int a, int b) { return b == 0 ? a : GCD(b , a % b); }
  53. ll LCM(ll a, ll b) { return a * (b/GCD(a, b)); }
  54.  
  55. bool isprime[1000001];
  56. vll seive()
  57. {
  58.     clr(isprime,true);
  59.     isprime[0]=false;
  60.     isprime[1]=false;
  61.     for(int i=4;i<MAXn;i+=2) isprime[i]=false;
  62.     for(int i=3;i*i<MAXn;i+=2){
  63.         if(!isprime[i]) continue;
  64.         for(int j=i*i;j<MAXn;j+=i) isprime[j]=false;
  65.     }
  66.     vll primes;
  67.     primes.pb(2LL);
  68.     for(int i=3;i<MAXn;i+=2){
  69.         if(isprime[i]) primes.pb((ll)i);
  70.     }
  71.     return primes;
  72. }
  73.  
  74. vll seg_seive(ll l,ll r,vll sprime)
  75. {
  76.     clr(isprime,true);
  77.     for(int i=0; sprime[i]*sprime[i]<=r; i++){
  78.         ll base = (l/sprime[i])*sprime[i];
  79.         if(base<l) base+=sprime[i];
  80.         for(ll j=base; j<=r; j+=sprime[i]) isprime[j-l]=false;
  81.         if(base==sprime[i]) isprime[base-l]=true;
  82.     }
  83.     vll primes;
  84.     for(int i=0; (ll)i<=r-l; i++){
  85.         if(isprime[i]){
  86.             if(i+l!=1LL) primes.pb((ll)(i+l));
  87.         }
  88.     }
  89.     return primes;
  90. }
  91.  
  92. int main()
  93. {
  94.     #ifdef VAMP
  95.         clock_t tStart = clock();
  96.         freopen("input.txt", "r", stdin);
  97.         freopen("output.txt", "w", stdout);
  98.     #endif
  99.  
  100.     vll sprime = seive();
  101.     ll l,r;
  102.     while(sc("%lld %lld",&l,&r)!=EOF){
  103.         vll primes = seg_seive(l,r,sprime);
  104.         ll mn=MAX,mx=0LL,a,b,c,d;
  105.         for(int i=1;i<primes.size();i++){
  106.             ll x = primes[i]-primes[i-1];
  107.             if(x<mn){
  108.                 mn = x;
  109.                 a=primes[i-1];
  110.                 b=primes[i];
  111.             }
  112.             if(x>mx){
  113.                 mx = x;
  114.                 c=primes[i-1];
  115.                 d=primes[i];
  116.             }
  117.         }
  118.         if(mx) pf("%lld,%lld are closest, %lld,%lld are most distant.\n",a,b,c,d);
  119.         else pf("There are no adjacent primes.\n");
  120.     }
  121.    
  122.  
  123.     #ifdef VAMP
  124.         fprintf(stderr, "\n>> Runtime: %.10fs\n", (double) (clock() - tStart) / CLOCKS_PER_SEC);
  125.     #endif
  126.  
  127.     return 0;
  128. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement