Advertisement
a_pramanik

GONE spoj ( digit dp)

Apr 1st, 2017
99
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 2.38 KB | None | 0 0
  1. ///:-)
  2. #include <bits/stdc++.h>
  3. #define pb push_back
  4. #define ll long long int
  5. #define inf 2000000000
  6. #define infLL 9000000000000000000
  7. #define infULL 18446744073709551615
  8. #define Aktaruzzaman using
  9. #define pi (2.0*acos(0.0))
  10. #define Pramanik namespace std;
  11. #define vsort(v)   sort(v.begin(),v.end())
  12. #define ull unsigned long long int
  13. #define mem(a, b) memset(a, b, sizeof a)
  14. #define cf 100009
  15. #define MOD 1000000007
  16. #define pii pair<int, int>
  17.  
  18. //int dx[]={0,1,0,-1};
  19. //int dy[]={1,0,-1,0};
  20. //int dx[]={1,1,0,-1,-1,-1,0,1};
  21. //int dy[]={0,1,1,1,0,-1,-1,-1};
  22.  
  23. template <class T> inline T gcd(T a,T b){if(b==0)return a;return gcd(b,a%b);}
  24. template< class T > inline T lcm(T a,T b) {a=abs(a);b=abs(b); return (a/gcd(a,b))*b;}
  25. template <class T> inline T is_prime(T a){if(a<2|a==4)return false;if(a==2||a==3||a==5)return true; T g=(T)sqrt(a)+1; for(T i=2; i<=g; i++){if(a%i==0)return false;}return false;}
  26. template <class T> inline T bigmod(T n,T p,T m){if(p==0)return 1; if(p%2)return ((n%m)*bigmod(n,p-1,m))%m; else {T c=bigmod(n,p/2,m);return ((c%m)*(c%m))%m;}}
  27.  
  28. Aktaruzzaman Pramanik
  29.  
  30. ll vis[11][85][3];
  31. ll a[10];
  32. ll prm[101];
  33. ll dp(ll pos, ll sum, ll choto)
  34. {
  35.     if(pos==0){
  36.         if(prm[sum]==0)return 1;
  37.         else return 0;
  38.     }
  39.     if(vis[pos][sum][choto]!=-1)return vis[pos][sum][choto];
  40.  
  41.     ll xx=0, hi;
  42.  
  43.     if(choto==1)hi=a[pos];
  44.     else hi=9;
  45.  
  46.     for(ll i=0; i<=hi; i++){
  47.         if(i==hi && choto==1)xx+=dp(pos-1, sum+i, 1);
  48.         else xx+=dp(pos-1, sum+i, 0);
  49.     }
  50.     return vis[pos][sum][choto]=xx;
  51. }
  52.  
  53. ll sv()
  54. {
  55.     prm[0]=1;prm[1]=1;
  56.     for(ll i=2; i<=100; i++){
  57.         if(prm[i]==0){
  58.             for(ll j = i*2; j<=100; j+=i){
  59.                 prm[j]=1;
  60.             }
  61.         }
  62.     }
  63. }
  64.  
  65.  
  66.  
  67. int main()
  68.  
  69. {
  70.     ios_base:: sync_with_stdio(0);
  71.     cin.tie(NULL);
  72.     cout.tie(NULL);
  73.  
  74.     ll t, i, j, k, l, r, sz, x, y;
  75.  
  76.     cin>>t;
  77.     sv();
  78.     while(t--){
  79.         cin>>l>>r;
  80.  
  81.         sz = 0;
  82.         k = l-1;
  83.  
  84.         while(k>0){
  85.             a[++sz]=k%10ll;
  86.             k/=10ll;
  87.         }
  88.  
  89.         mem(vis, -1);
  90.         x = dp(sz, 0, 1);
  91.  
  92.         sz = 0;
  93.         k = r;
  94.  
  95.         while(k>0){
  96.             a[++sz]=k%10ll;
  97.             k/=10ll;
  98.         }
  99. //        for(i=sz; i>=1; i--)cout<<a[i];
  100. //        cout<<endl;
  101.  
  102.         mem(vis, -1);
  103.  
  104.         y = dp(sz, 0, 1);
  105.  
  106.         cout<<y-x<<endl;
  107.  
  108.     }
  109.  
  110.  
  111.  
  112. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement