注册 登录  
 加关注

网易博客网站关停、迁移的公告:

将从2018年11月30日00:00起正式停止网易博客运营
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

GreenCloudS

 
 
 

日志

 
 

BZOJ-2005: [Noi2010]能量采集(莫比乌斯反演)  

2014-02-24 15:32:00|  分类: oi,bzoj,数论,莫 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2005


BZOJ-2005: [Noi2010]能量采集(莫比乌斯反演) - cjjlsdy - cjjlsdy的博客


代码:

  • #include <cstdio>
  • #include <algorithm>
  • #include <cstring>
  •  
  • using namespace std ;
  •  
  • #define ll long long
  • #define Sum( l , r ) ( sum[ r ] - sum[ l - 1 ] )
  • #define maxn 100100
  •  
  • bool f[ maxn ] ;
  • int u[ maxn ] , n , m , sum[ maxn ] ;
  • ll ans = 0 ;
  •  
  • void Init(  ) {
  •     memset( f , true , sizeof( f ) ) ;
  •     f[ 0 ] = f[ 1 ] = false , u[ 1 ] = 1 ; 
  •     for ( int i = 1 ; i ++ < min( n , m ) ; ) if ( f[ i ] ) {
  •         u[ i ] = - 1 ; 
  •         for ( int j = i << 1 ; j <= min( n , m ) ; j += i ) {
  •             f[ j ] = false ;
  •             if ( ! ( ( j / i ) % i ) ) u[ j ] = 0 ; else u[ j ] = u[ j / i ] * - 1 ;
  •         }
  •     }
  •     sum[ 0 ] = 0 ;
  •     for ( int i = 0 ; i ++ < min( n , m ) ; ) sum[ i ] = sum[ i - 1 ] + u[ i ] ;
  • }
  •  
  • ll query( int a , int b ) {
  •     ll ret = 0 ;
  •     for ( int i = 1 ; i <= min( a , b ) ; ) {
  •         int pos = min( a / ( a / i ) , b / ( b / i ) ) ;
  •         ret += ( ll )( Sum( i , pos ) ) * ( ll )( a / i ) * ( ll )( b / i ) ;
  •         i = pos + 1 ;
  •     }
  •     return ret ;
  • }
  •  
  • int main(  ) {
  •     scanf( "%d%d" , &n , &m ) ;
  •     Init(  ) ;
  •     for ( int i = 0 ; i ++ < min( n , m ) ; ) {
  •         ans += ( ll )( i ) * query( n / i , m / i ) * 2 ;
  •     }
  •     ans -= ( ll )( n ) * ( ll )( m ) ;
  •     printf( "%lld\n" , ans ) ;
  •     return 0 ;
  • }


  评论这张
 
阅读(3)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018