Project Euler 012 Highly divisible triangular number

最新推荐文章于 2020-09-13 17:00:38 发布

Owaski

最新推荐文章于 2020-09-13 17:00:38 发布

阅读量328

点赞数

分类专栏： PE 文章标签：数学

本文链接：https://blog.csdn.net/qq_20118433/article/details/55224474

版权

PE 专栏收录该内容

14 篇文章 0 订阅

订阅专栏

题意：求第一个满足约数个数 $>500$ 个的三角形数。
分析：考虑第 $n$ 个三角形数是 $\frac{n(n+1)}{2}$ ，而 $n$ 和 $n+1$ 互质，因此有

d (n ( n + 1 ) 2) = ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ d (n 2) \cdot d (n + 1), n 为 偶 数 d (n) \cdot d (n + 1 2), n 为 偶 数

$d\left(\frac{n(n+1)}{2}\right)= \left\{\begin{matrix} d(\frac{n}{2})\cdot d(n+1), n为偶数 \\ d(n)\cdot d(\frac{n+1}{2}), n为偶数 \end{matrix}\right.$
因此先预处理出

d $d$ ，然后枚举

n $n$ 即可。

#include <bits/stdc++.h>

#define ll long long

#define pii std::pair<int,int>
#define mp std::make_pair
#define fi first
#define se second

#define SZ(x) (int)(x).size()
#define pb push_back

template<class T>inline void chkmax(T &x, const T &y) {if(x < y) x = y;}
template<class T>inline void chkmin(T &x, const T &y) {if(x > y) x = y;}

template<class T>
inline void read(T &x) {
    char c;int f = 1;x = 0;
    while(((c=getchar()) < '0' || c > '9') && c != '-');
    if(c == '-') f = -1;else x = c-'0';
    while((c=getchar()) >= '0' && c <= '9') x = x*10+c-'0';
    x *= f;
}
static int outn;
static char out[(int)2e7];
template<class T>
inline void write(T x) {
    if(x < 0) out[outn++] = '-', x = -x;
    if(x) {
        static int tmpn;
        static char tmp[20];
        tmpn = 0;
        while(x) tmp[tmpn++] = x%10+'0', x /= 10;
        while(tmpn) out[outn++] = tmp[--tmpn];
    }
    else out[outn++] = '0';
}
const int N = 1e7;

int p[N + 9], pn;
bool is[N + 9];
int d[N + 9], g[N + 9];
// d[i]表示i的约数个数，g[i]表示i的最小质因子的次数

inline void sieve(int n) {
    d[1] = 1;
    for(int i = 2; i <= n; ++i) {
        if(!is[i]) p[++pn] = i, d[i] = 2, g[i] = 1;
        for(int j = 1; j <= pn && p[j] * i <= n; ++j) {
            int nex = p[j] * i;
            is[nex] = true;
            if(i % p[j] == 0) {
                g[nex] = g[i] + 1;
                d[nex] = d[i] / (g[i] + 1) * (g[nex] + 1);
                break;
            }
            d[nex] = d[i] * 2, g[nex] = 1;
        }
    }
}

int main() {
    sieve(N);
    for(int n = 1; ; ++n) {
        int u = n, v = n + 1;
        if(u & 1) v >>= 1;
        else u >>= 1;
        if(d[u] * d[v] > 500) {
            printf("%lld\n", (ll)u * v);
            break;
        }
    }
    return 0;
}