【BZOJ】【2005】【NOI2010】能量采集

欧拉函数



  玛雅,我应该先看看JZP的论文的……贾志鹏《线性筛法与积性函数》例题一

  这题的做法……仔细想下可以得到:$ans=2*\sum_{a=1}^n\sum_{b=1}^m gcd(a,b)-n*m$

  那么重点就在于算$\sum_{a=1}^n\sum_{b=1}^m gcd(a,b)$这个东西



copy一下JZP的推导过程:

$$ \sum_{a=1}^n \sum_{b=1}^m gcd(a,b)=&\sum_{a=1}^n \sum_{b=1}^m \sum_{d|gcd(a,b)} \varphi(d) $$

$$ = \sum_{a=1}^n \sum_{b=1}^m \sum_{d|a && d|b} \varphi(d) \\ = \sum \varphi(d) \sum_{1 \leq a \leq n && d|a} \sum{1 \leq b \leq m && d|b} 1 \\ = \sum \varphi(d) ( \sum_{1 \leq a \leq n && d|a} 1) * ( \sum_{1 \leq b \leq m && d|b} 1) \\ = \sum \varphi(d) \frac{n}{d} \frac{m}{d}  \]

 1 /**************************************************************
 2     Problem: 2005
 3     User: Tunix
 4     Language: C++
 5     Result: Accepted
 6     Time:40 ms
 7     Memory:2152 kb
 8 ****************************************************************/
 9
10 //BZOJ 2005
11 #include<cstdio>
12 #include<cstring>
13 #include<cstdlib>
14 #include<iostream>
15 #include<algorithm>
16 #define rep(i,n) for(int i=0;i<n;++i)
17 #define F(i,j,n) for(int i=j;i<=n;++i)
18 #define D(i,j,n) for(int i=j;i>=n;--i)
19 using namespace std;
20 typedef long long LL;
21 inline int getint(){
22     int r=1,v=0; char ch=getchar();
23     for(;!isdigit(ch);ch=getchar()) if(ch==‘-‘)r=-1;
24     for(; isdigit(ch);ch=getchar()) v=v*10+ch-‘0‘;
25     return r*v;
26 }
27 const int N=1e5+10,INF=~0u>>2;
28 /*****************template**********************/
29 int phi[N],prime[N],tot,n,m;
30 bool check[N];
31 void getphi(int n){
32     memset(check,0,sizeof check);
33     phi[1]=1;
34     int tot=0;
35     F(i,2,n){
36         if(!check[i]){
37             prime[++tot]=i;
38             phi[i]=i-1;
39         }
40         F(j,1,tot){
41             if(i*prime[j]>n) break;
42             check[i*prime[j]]=1;
43             if(i%prime[j]==0){
44                 phi[i*prime[j]]=phi[i]*prime[j];
45                 break;
46             }
47             else phi[i*prime[j]]=phi[i]*(prime[j]-1);
48         }
49     }
50 }
51 int main(){
52     n=getint(); m=getint();
53     if (n>m) swap(n,m);
54     getphi(N-1);
55     LL ans=0;
56     F(i,1,n)
57         ans+=(LL)phi[i]*(n/i)*(m/i);
58     printf("%lld\n",ans*2-(LL)n*m);
59     return 0;
60 }

时间: 2024-10-03 13:46:54

【BZOJ】【2005】【NOI2010】能量采集的相关文章

BZOJ 2005: [Noi2010]能量采集( 数论 + 容斥原理 )

一个点(x, y)的能量损失为 (gcd(x, y) - 1) * 2 + 1 = gcd(x, y) *  2 - 1. 设g(i)为 gcd(x, y) = i ( 1 <= x <= n, 1 <= y <= m ) 的数对(x, y)个数. 这个不好求, 考虑容斥, 设f(i) 为含有公因数 i 的数对(x, y)(1 <= x <= n, 1 <= y <= m)个数 , 显然f(i) = (n / i) * (m / i). 则 g(i) = f

bzoj 2005: [Noi2010]能量采集 筛法||欧拉||莫比乌斯

2005: [Noi2010]能量采集 Time Limit: 10 Sec  Memory Limit: 552 MB[Submit][Status][Discuss] Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后, 栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列 有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x的范

BZOJ 2005 [Noi2010]能量采集 (容斥)

[Noi2010]能量采集 Time Limit: 10 Sec  Memory Limit: 552 MB Submit: 2324  Solved: 1387 [Submit][Status][Discuss] Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以

bzoj 2005: [Noi2010]能量采集

1 #include<cstdio> 2 #include<iostream> 3 #define ll long long 4 using namespace std; 5 ll f[100009],n,m,ans; 6 int main() 7 { 8 scanf("%d%d",&n,&m); 9 if(n>m) 10 swap(n,m); 11 for(int i=n;i;i--) 12 { 13 f[i]=(n/i)*(m/i);

BZOJ 2005 NOI2010 能量采集 数论

题目大意:给定n和m,求Σ(1<=i<=n)Σ(1<=j<=m)GCD(i,j)*2-1 i和j的限制不同,传统的线性筛法失效了,这里我们考虑容斥原理 令f[x]为GCD(i,j)=x的数对(i,j)的个数,这个不是很好求 我们令g[x]为存在公因数=x的数对(i,j)的个数(注意不是最大公因数!),显然有g[x]=(n/x)*(m/x) 但是这些数对中有一些的最大公因数为2d,3d,4d,我们要把他们减掉 于是最终f[x]=(n/x)*(m/x)-Σ(2*x<=i*x<

2005: [Noi2010]能量采集 - BZOJ

Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x的范围是1至n,表示是在第x列,y的范围是1至m,表示是在第x列的第y棵. 由于能量汇集机器较大,不便移动,栋栋将它放在了一个角上,坐标正好是(0, 0). 能量汇集机器在汇集的过

2005: [Noi2010]能量采集

2005: [Noi2010]能量采集 Time Limit: 10 Sec  Memory Limit: 552 MBSubmit: 1831  Solved: 1086[Submit][Status] Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标

BZOJ 2005 [Noi2010]能量採集 (容斥)

[Noi2010]能量採集 Time Limit: 10 Sec  Memory Limit: 552 MB Submit: 2324  Solved: 1387 [Submit][Status][Discuss] Description 栋栋有一块长方形的地.他在地上种了一种能量植物,这样的植物能够採集太阳光的能量.在这些植物採集能量后,栋栋再使用一个能量汇集机器把这些植物採集到的能量汇集到一起. 栋栋的植物种得很整齐.一共同拥有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物.栋栋

[NOI2010]能量采集

469. [NOI2010]能量采集 ★★☆   输入文件:energy2010.in   输出文件:energy2010.out   简单对比时间限制:1 s   内存限制:512 MB [问题描述] 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x

【BZOJ2005】[Noi2010]能量采集 欧拉函数

[BZOJ2005][Noi2010]能量采集 Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x的范围是1至n,表示是在第x列,y的范围是1至m,表示是在第x列的第y棵. 由于能量汇集机器较大,不便移动,栋栋将它放在了一个角上,