Gym - 101652P B - Fear Factoring (数学题)

原题地址:http://codeforces.com/gym/101652/attachments

可以参考这篇除法分块博客https://www.cnblogs.com/Tony-Double-Sky/p/9490077.html

下面是我的代码以及一些注释

#include <bits/stdc++.h>
#include <cmath>
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#include <vector>
#include <stack>
#include <set>
#include <map>
#include <cctype>
#define eps 1e-8
#define INF 0x3f3f3f3f
#define PI acos(-1)
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define CLR(x,y) memset((x),y,sizeof(x))
#define fuck(x) cerr << #x << "=" << x << endl

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int seed = 131;
const int maxn = 1e5 + 5;
const int mod = 998244353;

ull cal(ull n) {
    ull i, j;
    ull ans = 0;
    for (i = 1; i <= n; i = j + 1) {//i是因子,枚举的是因子
        j = n / (n / i);
        /*
        j是在区间[1,n]内和因子i出现次数相同的最大的因子
        */
        ans += (i + j) * (j - i + 1) * (n / i) / 2;
        /*
        i+j 是首项+尾项
        j-i+1 是项数
        n/i 是出现的次数
        */
    }
    return ans;
}
int main() {
    ull a, b;
    while (cin >> a >> b) {
        cout << cal(b) - cal(a - 1) << endl;
    }
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值