BZOJ2005 NOI2010 能量采集


Description

栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量。在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起。 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x的范围是1至n,表示是在第x列,y的范围是1至m,表示是在第x列的第y棵。 由于能量汇集机器较大,不便移动,栋栋将它放在了一个角上,坐标正好是(0, 0)。 能量汇集机器在汇集的过程中有一定的能量损失。如果一棵植物与能量汇集机器连接而成的线段上有k棵植物,则能量的损失为2k + 1。例如,当能量汇集机器收集坐标为(2, 4)的植物时,由于连接线段上存在一棵植物(1, 2),会产生3的能量损失。注意,如果一棵植物与能量汇集机器连接的线段上没有植物,则能量损失为1。现在要计算总的能量损失。 下面给出了一个能量采集的例子,其中n = 5,m = 4,一共有20棵植物,在每棵植物上标明了能量汇集机器收集它的能量时产生的能量损失。 在这个例子中,总共产生了36的能量损失。

Input

仅包含一行,为两个整数n和m。

Output

仅包含一个整数,表示总共产生的能量损失。

Sample Input

【样例输入1】
5 4
【样例输入2】
3 4

Sample Output

【样例输出1】
36
【样例输出2】
20

HINT

对于100%的数据:1 ≤ n, m ≤ 100,000。



#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define N 100010
LL n,m,tot=0,ans=0;
LL pri[N],mu[N],F[N];
bool mark[N]={0};
void init(){
mu[1]=1;
for(int i=2;i<N;i++){
if(!mark[i])pri[++tot]=i,mu[i]=-1;
for(int j=1;j<=tot&&pri[j]*i<N;j++){
mark[pri[j]*i]=1;
if(i%pri[j]==0){
mu[i*pri[j]]=0;
break;
}else mu[i*pri[j]]=-mu[i];
}
}
for(int i=1;i<N;i++)F[i]=F[i-1]+mu[i];
}
LL solve(int d){
LL res=0,n1=n/d,m1=m/d,up=min(n1,m1);
for(int i=1,j;i<=up;i=j+1){
j=min(n1/(n1/i),m1/(m1/i));
res+=(F[j]-F[i-1])*(n1/i)*(m1/i);
}
return res*d;
}
int main(){
init();
scanf("%lld%lld",&n,&m);
int up=min(n,m);
for(int i=1;i<=up;i++)ans+=solve(i);
ans*=2;
ans-=n*m;
printf("%lld",ans);
return 0;
}

但是我们发现这样做似乎不是最优秀的

然后定义k=d∗pk=d∗p
转换一下:
ans=−n∗m+2∗∑min(n,m)d=1d∑min(n,m)d|kμ(k/d)⌊nk⌋⌊mk⌋ans=−n∗m+2∗∑d=1min(n,m)d∑d|kmin(n,m)μ(k/d)⌊nk⌋⌊mk⌋

把k提到前面枚举:

ans=−n∗m+2∗∑min(n,m)k=1∑min(n,m)d|kdμ(k/d)⌊nk⌋⌊mk⌋ans=−n∗m+2∗∑k=1min(n,m)∑d|kmin(n,m)dμ(k/d)⌊nk⌋⌊mk⌋

ans=−n∗m+2∗∑min(n,m)k=1∑min(n,m)d|kkdμ(d)⌊nk⌋⌊mk⌋ans=−n∗m+2∗∑k=1min(n,m)∑d|kmin(n,m)kdμ(d)⌊nk⌋⌊mk⌋

ans=−n∗m+2∗∑min(n,m)k=1⌊nk⌋⌊mk⌋∑min(n,m)d|kkdμ(d)ans=−n∗m+2∗∑k=1min(n,m)⌊nk⌋⌊mk⌋∑d|kmin(n,m)kdμ(d)

然后我们发现∑min(n,m)d|kkdμ(d)∑d|kmin(n,m)kdμ(d)是可以线性筛的

时间复杂度O(n)预处理+O(sqrt(n))查询O(n)预处理+O(sqrt(n))查询

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define N 100010
LL n,m,tot=0,ans=0;
LL pri[N],mu[N],F[N];
bool mark[N]={0};
void init(){
mu[1]=1;
for(int i=2;i<N;i++){
if(!mark[i])pri[++tot]=i,mu[i]=-1;
for(int j=1;j<=tot&&pri[j]*i<N;j++){
mark[pri[j]*i]=1;
if(i%pri[j]==0){
mu[i*pri[j]]=0;
break;
}else mu[i*pri[j]]=-mu[i];
}
}
for(int i=1;i<N;i++)
for(int j=1;j*i<N;j++)
F[i*j]+=mu[i]*j;
for(int i=1;i<N;i++)F[i]+=F[i-1];
}
int main(){
init();
scanf("%lld%lld",&n,&m);
int up=min(n,m);
for(int i=1,j;i<=up;i=j+1){
j=min(n/(n/i),m/(m/i));
ans+=(F[j]-F[i-1])*(n/i)*(m/i);
}
ans*=2;
ans-=n*m;
printf("%lld",ans);
return 0;
}

然后我们发现可能代码2比代码1慢?为什么呢?
因为预处理的时候失去了线性的性质

但是当查询很多的时候代码优势就得以凸显

各取所需吧

BZOJ2005 NOI2010 能量采集 【莫比乌斯反演】的更多相关文章

  1. BZOJ2005: [Noi2010]能量采集 莫比乌斯反演的另一种方法——nlogn筛

    分析:http://www.cnblogs.com/huhuuu/archive/2011/11/25/2263803.html 注:从这个题收获了两点 1,第一象限(x,y)到(0,0)的线段上整点 ...

  2. BZOJ2005:[NOI2010]能量采集(莫比乌斯反演,欧拉函数)

    Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得 ...

  3. luogu1447 [NOI2010]能量采集 莫比乌斯反演

    link 冬令营考炸了,我这个菜鸡只好颓废数学题了 NOI2010能量采集 由题意可以写出式子: \(\sum_{i=1}^n\sum_{j=1}^m(2\gcd(i,j)-1)\) \(=2\sum ...

  4. BZOJ 2005: [Noi2010]能量采集 [莫比乌斯反演]

    题意:\((0,0)\)到\((x,y),\ x \le n, y \le m\)连线上的整点数\(*2-1\)的和 \((0,0)\)到\((a,b)\)的整点数就是\(gcd(a,b)\) 因为. ...

  5. bzoj2005: [Noi2010]能量采集

    lsj师兄的题解 一个点(x, y)的能量损失为 (gcd(x, y) - 1) * 2 + 1 = gcd(x, y) *  2 - 1. 设g(i)为 gcd(x, y) = i ( 1 < ...

  6. [BZOJ2005][Noi2010]能量采集 容斥+数论

    2005: [Noi2010]能量采集 Time Limit: 10 Sec  Memory Limit: 552 MBSubmit: 4394  Solved: 2624[Submit][Statu ...

  7. 【莫比乌斯反演】BZOJ2005 [NOI2010]能量采集

    Description 求sigma gcd(x,y)*2-1,1<=x<=n, 1<=y<=m.n, m<=1e5. Solution f(n)为gcd正好是n的(x, ...

  8. BZOJ2005: [Noi2010]能量采集(容斥原理 莫比乌斯反演)

    Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 4727  Solved: 2877[Submit][Status][Discuss] Descript ...

  9. [bzoj2005][Noi2010][能量采集] (容斥 or 欧拉函数)

    Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后, 栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种 ...

随机推荐

  1. WinCE数据通讯之Web Service分包传输篇

    前面写过<WinCE数据通讯之Web Service篇>那篇对于数据量不是很大的情况下单包传输是可以了,但是对于大数据量的情况下WinCE终端的内存往往会在解包或者接受数据时产生内存溢出. ...

  2. vue element new vue const

    new Vue{ el:"app", } ========= const app = new Vue({ router, data:{....} }).$mount('#app') ...

  3. [PostgreSql]PostgreSql创建函数及函数的执行

    1.准备好创建函数的脚本 -- FUNCTION: public.dynamic_placelist_select(integer, timestamp without time zone) -- D ...

  4. 交通部道路运输车辆卫星定位系统部标JTT808、809、796标准大全

    无论是开发GPS设备硬件还是开发应用软件,都要面临一个标准,这个标准就是国家交通部发布的道路运输车辆卫星定位系统部标认证标准,它涵盖了GPS硬件设备参数.功能标准,也包括了设备上传到应用平台的协议标准 ...

  5. [嵌入式培训 笔记]-----Vim编辑器使用简介

    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 第一讲小结 1. 光标在屏幕文本中的移动既 ...

  6. C++调用Python脚本中的函数

    1.环境配置 安装完python后,把python的include和lib拷贝到自己的工程目录下 然后在工程中包括进去 2.例子 先写一个python的测试脚本,如下 这个脚本里面定义了两个函数Hel ...

  7. (2) iOS开发之UI处理-UILabel篇

    我们经常要根据内容去动态计算控件的高度,比如一个UILabel控件,常常要显示多行内容,并且计算出总高度,如果每个UILabel要多行显示,都要写这么一段代码是非常痛苦的,看代码如下:     我想大 ...

  8. restframework api(基础1)

    最近项目忙成狗,都没时间好好看看开发的东西了,正好最近开始看rest api的东西,真是个好东西啊.可以前后端分离,但是在学习的过程中,遇到各种问题.还是基础不够扎实.本次rest api的会一边遇坑 ...

  9. 伪基站,卒于5G——本质上是基于网络和UE辅助的伪基站检测,就是将相邻基站的CI、信号强度等信息通过测量报告上报给网络,网络结合网络拓扑、配置信息等相关数据,对所有数据进行综合分析,确认在某个区域中是否存在伪基站

    伪基站,卒于5G from:https://www.huxiu.com/article/251252.html?h_s=h8 2018-07-05 21:58收藏27评论6社交通讯     本文来自微 ...

  10. MySql设计规范及SQL索引优化【呕心之作】

    数据库及表结构基本设计规范 1. 所有表必须使用Innodb存储引擎 没有特殊要求(即Innodb无法满足的功能如:列存储,存储空间数据等)的情况下,所有表必须使用Innodb存储引擎(mysql5. ...