【bzoj1853】: [Scoi2010]幸运数字 数论-容斥原理

【bzoj1853】: [Scoi2010]幸运数字

预处理出所有幸运数字然后容斥原理

但是幸运数字是2^logn个数的 直接搞会炸

所以把成倍数的处理掉

然后发现还是会T 所以数字要从大到小处理会快很多 我也不知道为什么

还有就是剪枝的时候要开double 不然爆longlong

 /* http://www.cnblogs.com/karl07/ */
 #include <cstdlib>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <map>
 #include <algorithm>
 using namespace std;
 #define ll long long

 const ll N=;
 ll num[];
 ll a,b,cnt,ans,cc;

 ll gcd(ll a,ll b){ return b== ? a : gcd(b,a%b);}
 ll lcm(ll a,ll b){ return a/gcd(a,b)*b;}

 void dfs(ll x,ll n){
     if (x>n) return;
     num[++cnt]=x;
     dfs(x*+,n);
     dfs(x*+,n);
 }

 void rc(ll c,ll x,ll n,ll nw){
     if (x==){
         if (c==) return;
         if (c&) ans+=n/nw; else ans-=n/nw;
         return;
     }
     if ((double)nw/gcd(num[x],nw)*num[x] <= n) rc(c+,x-,n,lcm(nw,num[x]));
     rc(c,x-,n,nw);
 }

 ll work(ll n){
     if (n<) return ;
     cnt=;ans=;
     dfs(,n);dfs(,n);cc=cnt;
     num[]=;
     sort(num+,num+cnt+);
     for (int i=;i<=cnt;i++){
         for (int j=i+;j<=cnt && num[i]!=N;j++){
             if (num[j]%num[i]==) num[j]=N,cc--;
         }
     }
     sort(num+,num+cnt+);
     rc(,cc-,n,num[cc]);
     rc(,cc-,n,);
     return ans;
 }

 int main(){
     scanf("%lld%lld",&a,&b);
     printf("%lld\n",work(b)-work(a-));
     return ;
 }