三角数因子数

#include <stdio.h>
#include <string.h>
#include <math.h>
#define __STDC_FORMAT_MACROS
#include <inttypes.h>
#define max_num 10000000
int64_t prime[max_num] = { 0 };
int64_t num[max_num] = { 0 };
int64_t sum_min[max_num] = { 0 };
int primes()
{
    for(int64_t i = 2; i <= max_num; i++)
    {
        if(!prime[i])
        {
            prime[++prime[0]] = i;
            num[i] = 2;
            sum_min[i] = 1;
        }
        for(int64_t j = 1; j <= prime[0]; j++)
        {
            if(prime[j] * i > max_num) break;
            prime[i*prime[j]] = 1;
            if(i%prime[j]==0)
            {
                num[i * prime[j]] = num[i] / (sum_min[i] + 1) * (sum_min[i] + 2);
                sum_min[i * prime[j]] = sum_min[i] + 1;
                break;
            }
            else
            {
                num[i * prime[j]] = num[i] * num[prime[j]];
                sum_min[i * prime[j]] = sum_min[prime[j]];
            }
        }
    }
    /*for(int i = 2; i <= 1000; i++)
    {
        printf("%5d:%8"PRId64,i,num[i]);
        if(i%3==0)
        {
            printf("\n");
        }
    }*/
}


int main()
{
    primes();
    int64_t i=1;
    while(i++)
    {
        int64_t sum;


        if(i%2==0)
        {
            sum = num[i/2]*num[i-1];
        }
        else
        {
            sum = num[(i-1)/2]* num[i];
        }
        if(sum > 500) break;
    }
    printf("%"PRId64"\n",i*(i-1)/2);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值