Project-Euler-021线性求约数和

021题：

题意：

​ 记d(n)为n的所有真因数（小于n且整除n的正整数）之和
​ 如果d(a) = b且d(b) = a，且a ≠ b，那么a和b构成一个亲和数对，a和b被称为亲和数

​ 例如，220的真因数包括1、2、4、5、10、11、20、22、44、55和100，因此d(220) = 284；而284的真因数包括1、2、4、71和142，因此d(284) = 220

​ 求所有小于10000的亲和数的和

思路：

​ 在数论中我们称d(n)这样的函数为约数和函数。具体有以下几个性质：

​ 1.约数和函数为积性函数，因此满足任意两个互质的数a，b，都有d(a✖️b) = d(a) ✖️d(b)

​ 2.若b = a ✖️ a中的最小素因子p1 ，也就是

a = p1^a1 * p2^a2 * ... * pn^an
b = p1 * p1^a1 * p2^a2 * ... * pn^an

则：

d(b) = d(a) * (1 - p1^(a1 + 2)) / (1 - p1^(a1 + 1))	

​

​ 3.对于素数来说，它的约数只有1和它本身，所以：

d(n) = n + 1;    //当n为素数

​

​ 4.对于求解一般的约数和函数，我们有一个约数和定理：

d(n)=(p1^0+p1^1+p1^2+…p1^a1)(p2^0+p2^1+p2^2+…p2^a2)…(pk^0+pk^1+pk^2+…pk^ak)
//n = p1^a1 * p2^a2 *...* pk^ak且p都为素数

​ 通过等比数列求和公式化简一下，为：

d(n)=Π (1 - pi^(ai + 1)) / (1 - pi)		
//pi为它的素因子，ai为该种素因子的个数

​ 我们就可以根据n的素因子种类pi和每种的个数ai来求约数和了。关于pi和ai我们可以利用之前的线性筛求因子个数框架来求取 抛个链接：走你

​

代码：

#include <stdio.h>
#include <string.h>
#include <inttypes.h>
#define ll int64_t
#define d int32_t
#define f double
#define r return
#define N (d)1e4
#define mem(a) memset(a, 0, sizeof(a))
#define scanfd(a) scanf("%d", &a)
#define scanfl(a) scanf("%lld", &a)
#define printfl(a) printf("%lld\n", a)
#define printfd(a) printf("%d\n", a)
#define For(i, star, endd) for (d i = star; i <= endd; i++)
#define Forr(i, endd, star) for (d i = endd; i >= star; i--)

ll prime[N + 5];     //存全部的素数
ll dd[N + 5];        //用来存约数和
ll book[N + 5];     //book[i]用来存储i的最小素因子幂次方值，如book[24] = 8

void init() {
    mem(prime);
    mem(dd);
    mem(book);
    For(i, 2, N) {
        if(!book[i]) {
            prime[++prime[0]] = i;
            dd[i] = i + 1;              //素数的约数只有1和它本身
            book[i] = i;                //素数的最小素因子幂次方值等于它本身
        }
        For(j, 1, prime[0]) {
            ll t = prime[j] * i;
            if (t > N) break;
            if (i % prime[j] == 0) {        //当为性质2的时候
                book[t] = book[i] * prime[j];
                dd[t] = dd[i] * (book[i] * prime[j] * prime[j] - 1) / (book[i] * prime[j] - 1);
                break;
            } else {                        //当为性质1的时候
                dd[t] = dd[prime[j]] * dd[i];
                book[t] = prime[j];
            }
        }
    }
    r;
}

d main () {
    init();
    ll ans = 0;
    For(i, 2, N) {
        dd[i] -= i;
    }
    For(i, 0, N) {
        if (dd[i] != i && dd[i] <= N && i == dd[dd[i]]) {
            ans += i;
        }
    }
    printfl(ans);
    r 0;
}

最后结果是：31626

转载请注明出处！！！

如果有写的不对或者不全面的地方 可通过主页的联系方式进行指正，谢谢