刷题记录（NC17450 因数个数和、NC218398 小G的约数、NC204916 浅尝辄止）(数论分块)

NC17450 因数个数和

题目链接

关键点：

1、题目要求求因数的个数，考虑到数据范围，可以将所求的x分为两部分，一部分<=sqrt(x),一部分>sqrt(x);

2、对于<=sqrt(x),就直接求该范围内的数做的贡献，即1-x中有多少个数字有该范围的因子，那么对于范围内的i，i的贡献为x/i，然后利用该i取求>sqrt(x)范围内的贡献，即找i*(刚刚好大于sqrt(x))还在x范围内的数字的个数，那么该数的范围为(sqrt(x)+1 -> x/i)

比方说求1-100的因数个数和

首先是看1的贡献，为100/1，共有100个数有1这个因子

然后是找1*(刚刚好大于sqrt(x))还在x范围内的数字数，满足要求的数字有1*11， 1*12...1*100

即11有因子11，12有因子12....个数为（100-11+1）这就是利用<=sqrt(x)的因子来计算>sqrt(x)的因子个数

再比方说2的贡献，为100/2，共有50个数有2这个因子

然后有2*11,2*12,...2*50，即22有11因子，24有12因子。。。，再加上(50-11+1)

以此类推，推到sqrt(100) = 10

这样就将所有因子数计算，为分块思想

3、注意个数要开long long

完整代码：

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll x, q;
int main()
{
    scanf("%lld", &q);
    while (q--)
    {
        scanf("%lld", &x);
        ll tmp = sqrt(x);
        ll t = tmp+1, total = 0;
        for (int i=1; i<=tmp; i++)
        {
            total += x/i;
            total += (x/i-t+1);
        }
        cout<<total<<endl;
    }
    
    
    
    return 0;
}

NC218398 小G的约数

题目链接

关键点：

1、该题与上一题相似，只不过要求求的为所有约数的和

2、一样的想法，分为两个部分来求，对于<=sqrt(n)，直接求该数在1-n的贡献，求出所有有该因数的个数，然后乘上该数，比方说求1的贡献(假设1<=sqrt(n))，即n/1*1即可，其他类似

3、对于大于sqrt(n)的数，还是一样利用<=sqrt(n)的数来间接的求，只不过这次求的次数改为求的和，即求一个等差序列的和，比如n = 100, sqrt(n) = 10

对于1来说，1*大于sqrt(n)的数中还在100范围内的有

11，12，13，14....100 ，个数为100/1 - (10+1) +1,首项为10+1,末项为100/1

那么对于任意的i<=sqrt(n)

其计算i的贡献和与i有关的>sqrt(n)的贡献的公式为：设sqrt(n) + 1 = t

n/i*i(i的贡献) + (t+n/i)*(n/i - t+1)/2(与i有关的大于sqrt(n)的贡献和，即等差数列求和)

完整代码：

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int n;
ll ans, rans;
int main()
{
    scanf("%d", &n);
    int s = sqrt(n);
    ll t = s+1;
    for (int i=1; i<=s; i++)
    {
        ll a = (n/i)*i;
        ll b = (t+n/i)*(n/i-t+1)/2;
        ans += a+b;
    }
    
    ll s1 = sqrt(ans);
    ll t1 = s1+1;
    for (int i=1; i<=s1; i++)
    {
        ll a = (ans/i)*i;
        ll b = (t1+ans/i)*(ans/i-t1+1)/2;
        rans += a+b;
    }
    cout<<rans<<endl;
    
    return 0;
}

NC204916 浅尝辄止

题目链接

关键点：

1、该题想法与上述类似

2、对于数n，其除以小于n的数即n/i，得数为1有很多种，得数为2，有很多种，采用分块思想

3、对于得数为1-sqrt(n)的都合起来求，对于得数大于sqrt(n)则暴力一个个求

4、那么如何找得数为小于等于sqrt(n)的个数呢，拿1来举例子，当n除以大于自己一半的数得数就为1，即范围在n/2+1 到 n/1 ,同理得数为2的数的范围在n/3+1 到 n/2，以此类推

完整代码：

class Solution {
public:
    /**
     * 
     * @param n long长整型 
     * @return int整型
     */
    int work(long long n) {
        // write code here
        int ans = 0;
        long long maxn, minn;
        for (int i=1; 1LL*i*i<=n; i++)
        {
            maxn = n/i;
            minn = n/(i+1)+1;
            ans = (ans + (maxn-minn+1)*i%998244353)%998244353;
        }
        for (int i=1; i<minn; i++)
            ans = (ans+n/i%998244353)%998244353;
        return ans;
    }
};