bzoj1968: [Ahoi2005]COMMON 约数研究解题报告

Michael_Bryant

于 2018-01-25 22:53:53 发布

阅读量338

点赞数 1

分类专栏：数论

本文链接：https://blog.csdn.net/Michael_Byrant/article/details/79167633

版权

数论专栏收录该内容

6 篇文章 0 订阅

订阅专栏

题目链接

传送门

一句话题意

输入n 输出从1到n每个数约数的个数和

想说的话

有的时候麻烦跟简单
没差在复杂度
就差在想法上面了
这里写图片描述
下面给出这三种做法的题解

题解I

我先说最慢但是最好想的吧
欧拉筛
筛一遍
筛的过程中维护素数不用我多说
跑了300多ms有点慢（相对于另外两种方法）
相对来说代码长度也长

代码I

#include <cmath>
#include <queue>
#include <cstdio>
#include <iomanip>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 1000001
#define ll long long
using namespace std;
ll prime[100001];
bool isprime[N];
ll t[N],d[N],n;
void su()
{
  memset(isprime,true,sizeof isprime);
  isprime[1]=false;
  d[1]=1,t[1]=1;
  for(int i=2;i<=n;i++)
  {
    if(isprime[i])
    {
      prime[++prime[0]]=i;
      d[i]=2,t[i]=1;
    }
    for(int j=1;j<=prime[0]&&i*prime[j]<=n;j++)
    {
      isprime[i*prime[j]]=false;
      if(i%prime[j]==0)
      {
        t[i*prime[j]]=t[i]+1;
        d[i*prime[j]]=d[i]/(t[i]+1)*(t[i*prime[j]]+1);
        break;
      }
      else
      {
        d[i*prime[j]]=d[i]*2;
        t[i*prime[j]]=1;
      }
    }
  }
}
inline ll read(){
    char ch=getchar();ll sum=0,w=1;
    while(!(ch>='0'&&ch<='9')) {if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=getchar();
    return sum*w;
}
int main()
{

  n=read();
  su();
  ll ans=0;
  for(ll i=1;i<=n;i++) ans+=d[i];
  cout<<ans<<endl;
}

题解II

考虑到我们要找约数一定要筛不好写
我们换个思路
对于一个数n 我们要找到从1-n所有数约数的个数和
那么我们是不是可以找一下从1-n每个数的倍数出现了多少个
每出现一个证明这个数作为了一次其他数的约数
我们统计这个值相对来讲更简单
代码长度巨短时间也快了好多

代码II

#include<cstdio>
int n,ans;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)ans+=n/i;
    printf("%d\n",ans);
    return 0;
}

没错就是这么干净利落

题解III

看了上面的题解你们一定在想剩下0ms那个做法是何方神圣
其实没有什么特别特别高大上的东西
就是把我们的算法II进行了一步优化
时间复杂度从O(n)变成O(sqrt(n))
这个优化是什么样的呢
其实我是看了zyf学姐的博客才知道这种方法
当时我想了很长时间为什么可以采用这种优化
这个优化的方式是这样的
如果把所有的⌊n/i⌋都列出来形成一个数列，你会发现有一些区间内的数都是一样的，而如果你把所有的⌊n/n/i⌋求出来的话，就得到了一个与上一个序列“相反”的序列（其实求这个序列的目的只是为了方便代码实现）。
然后每一次只需要求一下区间和就可以了，时间复杂度降低
其实这个很难想到
我觉得第二种方法已经是考场的极限了

代码III

重新写了一遍内存和长度又优化了
2543184 Michael_Bryant 1968 Accepted 820 kb 0 ms C++/Edit 242 B 2018-01-25 22:44:24

#include<cstdio>
using namespace std;
#define LL long long
LL n,ans;
int main(){
    scanf("%lld",&n);
    LL i,j;
    for (i=1,j=0;i<=n;i=j+1){
        j=n/(n/i);
        ans+=(j-i+1)*(n/i);
    }
    printf("%lld\n",ans);
}