XTU-OJ 1396-函数

题目描述

根据质因子唯一分解定理可知

n=∏i=1mpkii

,其中 pi 都是素数。

定义 g(n)=∑mi=1ki , 求 ∑bi=ag(i) 。

输入

第一行是一个整数T(1≤T≤10000),表示样例的个数。 以后每行两个整数a,b(1≤a≤b≤106)

输出

每行输出一个样例的结果。

样例输入

2
1 2
1 1000000

样例输出

1
3626619

解题思路:这题就和 1377 同根同源,几乎没什么变化,代码只需要改一行。

AC代码:

#include <stdio.h>

const int N = 1e6;
int vis[N+2] = {0};
int prime[80000] = {0};
void is_prime()
{
    for (int i = 2; i < N; i ++)                // 欧筛
    {
        if (!vis[i])
        {
            prime[++prime[0]] = i;
            vis[i] = 1;
        }
        for (int j = 1; j <= prime[0] && i <= N/prime[j]; j ++)
        {
            vis[i*prime[j]] = vis[i]+1;    
            if (i % prime[j] == 0)        break;
        }
    }
    for (int i = 2; i <= N; i ++)               // 前缀和
        vis[i] += vis[i-1];
}

int main()
{
    is_prime();
    int T;
    scanf("%d",&T);
    while ( T --)
    {
        int a,b;
        scanf("%d %d",&a,&b);
        printf("%d\n",vis[b]-vis[a-1]);
    }
    return 0;
}

<think>我们参考多个引用中的实现方式,编写一个判断素数函数,并且在主函数中调用它。 根据引用[3]和引用[4],我们可以选择两种方式: 方式1:在判断素数函数中直接输出结果(如引用[1]、引用[4]、引用[5]) 方式2:在判断素数函数中返回一个值(如0或1,表示是否为素数),然后在主函数根据返回值输出(如引用[3]) 这里我们采用第二种方式,因为这样可以使函数的功能更纯粹(只判断,不输出),更符合模块化设计。 素数判断的优化:只需要判断到sqrt(n)即可(如引用[3]),因为如果n不是素数,则它一定有一个因子小于等于sqrt(n)。 具体步骤: 1. 如果n小于2,则不是素数。 2. 2是素数。 3. 对于大于2的数,判断从2到sqrt(n)(取整)之间是否有因子,有则不是素数,否则是素数。 注意:需要包含头文件<math.h>以便使用sqrt函数。 代码实现: 函数原型:int is_prime(int n); 返回1表示素数,0表示非素数。 主函数中: 输入整数n,调用is_prime(n),根据返回值输出相应信息。 参考引用[3]的实现,但稍作改进(例如,将sqrt(n)提取出来,避免在循环中重复计算;同时注意循环条件等)。 注意:在循环中,如果n能被2整除,那么除了2以外,其他偶数都不是素数,但我们循环从2开始,所以可以处理。另外,我们可以在函数开始时单独判断2。 具体实现如下: </think>### C语言实现素数判断函数 下面是完整的C语言代码,实现判断素数函数并在主函数中输出结果: ```c #include <stdio.h> #include <math.h> // 包含sqrt函数 // 声明判断素数函数 int is_prime(int n); int main() { int num; printf("请输入一个整数: "); scanf("%d", &num); if (is_prime(num)) { printf("%d 是素数\n", num); } else { printf("%d 不是素数\n", num); } return 0; } // 定义判断素数函数 int is_prime(int n) { // 处理特殊情况: // 1. 小于等于1的数都不是素数 // 2. 2是素数 if (n <= 1) return 0; if (n == 2) return 1; // 偶数(除了2)都不是素数 if (n % 2 == 0) return 0; // 检查奇数因子,只需检查到sqrt(n)即可 int limit = (int)sqrt(n) + 1; // 计算检查上限 for (int i = 3; i <= limit; i += 2) { if (n % i == 0) { return 0; // 找到因子,不是素数 } } return 1; // 没有找到因子,是素数 } ``` ### 代码说明 1. **算法优化**: - 使用平方根优化:只需检查到$\sqrt{n}$即可,大大减少计算量 - 跳过偶数检查(除了2),减少一半的迭代次数 - 时间复杂度为$O(\sqrt{n})$ 2. **函数设计**: - `is_prime()`函数返回1表示素数,0表示非素数 - 函数独立完成素数判断,不负责输出 -函数负责输入输出逻辑,符合模块化设计原则 3. **关键处理**: - 特殊处理数字1(非素数- 特殊处理数字2(素数- 单独处理偶数以提高效率 ### 测试示例 ``` 请输入一个整数: 17 17 是素数 请输入一个整数: 100 100 不是素数 请输入一个整数: 2 2 是素数 请输入一个整数: 1 1 不是素数 ``` ### 算法优化原理 数学定理:如果$n$不是素数(合数),则它必定有一个质因子$\leq \sqrt{n}$[^3]。因此只需检查$2$到$\sqrt{n}$范围内的整数即可确定$n$是否为素数,无需检查到$n-1$。 时间复杂度:$O(\sqrt{n})$,空间复杂度:$O(1)$ ***
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值