注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

ydc的博客

 
 
 

日志

 
 

bzoj 2005(noi 2010) 能量采集  

2013-06-04 20:28:17|  分类: bzoj |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
simga(gcd(i,j))
这种题做过很多次了
从开始学信息学开始,直到现在,时间复杂度一步步降低
这题n ln n就可以水过去了
我特地学习了jzp的O(n)预处理O(n^0.5)单次询问的方法
想写点推导的……但是不会打simga的符号
所以就算了吧
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxn 100010
using namespace std;
typedef long long LL;
LL n,m,phi[maxn],sum[maxn];
void Prepare(int n)
{
static int prime[maxn];
int tot=0;
phi[1]=1;
for(int i=2;i<=n;i++)
{
if(phi[i]==0)
{
phi[i]=i-1;
prime[++tot]=i;
}
for(int j=1;j<=tot&&prime[j]<=n/i;j++)
{
if(i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
else
phi[i*prime[j]]=phi[i]*phi[prime[j]];
}
}
for(int i=1;i<=n;i++)
sum[i]=sum[i-1]+phi[i];
}
LL work(int n,int m)
{
if(n>m)
swap(n,m);
LL ans=0;
for(int i=1,nexti;i<=n;i=nexti+1)
{
int ta=n/i,tb=m/i;
nexti=min(n/ta,m/tb);
ans=ans+(sum[nexti]-sum[i-1])*ta*tb;
}
return ans;
}
int main()
{
Prepare(100000);
scanf("%lld %lld",&n,&m);
printf("%lld\n",work(n,m)*2-n*m);
return 0;
}


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

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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