BZOJ 1853

http://www.lydsy.com/JudgeOnline/problem.php?id=1853

岛娘在空间上发的题解就看了看果然被骗了。还以为是数位dp。

原来是容斥啊。好吧第一道正式的题目。

这题还不错，有几个小优化，像排序之类的啦很常见。

看了一下别人的代码，学习了一下姿势。

---

http://hi.baidu.com/greencloud/item/6f564ec036a3ba2ee90f2edb

这里直接把别人的代码贴过来。

 1 #include <iostream>
 2 #include <algorithm>
 3 #include <cstring>
 4  
 5 using namespace std ;
 6  
 7 #define rep( i , x , y ) for ( int i = x ; i <= y ; ++ i )
 8  
 9 typedef long long ll ;
  
 const int maxn =  ;
  
 ll a , b , num[ maxn ] ;
 int cnt =  ; 
  
 void make( ll now ) {
     if ( now > b ) return ; 
     num[ ++ cnt ] = now ;
     make( now *  +  ) , make( now *  +  ) ;
 }
  
 ll ans , n[ maxn ] ;
 int m =  ; 
  
 ll X ;
  
 ll gcd( ll x , ll y ) {
     if ( x < y ) swap( x , y ) ;
     ll k ;
     while ( y ) {
         k = y ;
         y = x % y ;
         x = k ;
     }
     return x ;
 }
  
 void dfs( int pos , int times , ll rec ) {
     if ( ! pos ) {
         if ( ! times ) return ;
         ll ret = ( b / rec ) - ( a / rec ) ;
         if ( times &  ) ans += ret ; else ans -= ret ;
         return ; 
     }
     ll temp = gcd( rec , n[ pos ] ) ;
     if ( rec / temp <= b / n[ pos ] ) dfs( pos -  , times +  , rec / temp * n[ pos ] ) ;
     dfs( pos -  , times , rec ) ;
 }
  
 int main(  ) {
     cin >> a >> b ; 
     make(  ) , make(  ) ;
     rep( i ,  , cnt ) if ( num[ i ] ) {
         rep( j ,  , cnt ) if ( i != j && ! ( num[ j ] % num[ i ] ) ) {
             num[ j ] =  ;
         }
     }
     memset( n ,  , sizeof( n ) ) ;
     rep( i ,  , cnt ) if ( num[ i ] ) {
         n[ ++ m ] = num[ i ] ;
     }
     sort( n +  , n + m +  ) ;
     -- a ;
     ans =  ;
     dfs( m ,  ,  ) ;
     cout << ans << endl ;
     return  ; 

68 }