SGU 491 约数

题意

传送门 SGU 491

题解

枚举 B B B y y y,枚举总次数为 O ( n log ⁡ n ) O(n\log n) O(nlogn)。若 A x + B y = N Ax+By = N Ax+By=N,则 A A A N − B y N-By NBy 的约数。对于固定的 B y By By,枚举满足 A < B A<B A<B 的约数。预处理 [ 1 , n ] [1,n] [1,n] 的约数,总数是 O ( n log ⁡ n ) O(n\log n) O(nlogn),故枚举约数均摊是 O ( log ⁡ n ) O(\log n) O(logn),总枚举次数 O ( n log ⁡ 2 n ) O(n\log^2 n) O(nlog2n)。对于固定 B B B 需要对 A A A 去重,此时枚举的 A A A 数量均摊意义下 O ( n / B log ⁡ n ) O(n/B\log n) O(n/Blogn),则排序去重约 O ( n log ⁡ 2 n ) − O ( n log ⁡ 3 n ) O(n\log^2 n)-O(n\log^3 n) O(nlog2n)O(nlog3n)

#include <bits/stdc++.h>
using namespace std;
#define pb push_back
const int MAXN = 1E5 + 5;
int N;
vector<int> fac[MAXN];

int main()
{
    ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    cin >> N;
    for (int i = 1; i <= N; ++i)
        for (int j = i; j <= N; j += i)
            fac[j].pb(i);
    int res = 0;
    for (int B = 2; B < N; ++B)
    {
        vector<int> rec;
        for (int y = 1; y <= N / B; ++y)
        {
            int t = N - B * y;
            for (int A : fac[t])
            {
                if (A >= B)
                    break;
                rec.pb(A);
            }   
        }
        sort(rec.begin(), rec.end());
        res += unique(rec.begin(), rec.end()) - rec.begin();
    }
    cout << res << '\n';
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值