hdu 5750

AC:

#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 100005
using namespace std;
bool check[maxn];
int pr[maxn / 5], mn[maxn];
int tt[maxn];
int tot;
int T;
void euler() {
    memset(check,0,sizeof(check));
    tot=0;
    //x=1时的函数值 
    for(int i=2;i<maxn;i++) {
        if(!check[i]) {
            pr[++tot]=i;
            mn[i] = i;
            tt[i] = 1;
            //质数函数值 
        }
        for(int j=1;j<=tot;j++) {
            if(i*pr[j]>maxn) break;
            check[i*pr[j]]=1;
            
            mn[i*pr[j]] = pr[j];
                
            if(i%pr[j]==0) {
                break;
            } 
        }
    } 
}
int main()
{
    euler();
    tt[0] = 0;
    for (int i = 1; i < maxn; ++ i) tt[i] += tt[i - 1];            //预处理出小于i的素数的个数 
    scanf("%d", &T);
    while (T --)
    {
        int n, d, ans;
        scanf("%d%d", &n, &d); n --;
        if (d >= maxn)
        {
            ans = n / d;
            for (int i = 1; pr[i] < ans; ++ i)
                if (d % pr[i] == 0) {ans = pr[i]; break;}
            printf("%d\n", tt[ans]);
        }
        else printf("%d\n", tt[min(n / d, mn[d])]);
    }
}

几乎一模一样但WA:

#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 100005
using namespace std;
bool isprime[maxn];
int pr[maxn / 5], mn[maxn];
int tt[maxn];
int tot;
int T;
void euler() {
    memset(isprime,1,sizeof(isprime));
    tot=0;
    //x=1时的函数值 
    for(int i=2;i<maxn;i++) {
        if(isprime[i]) {
            pr[++tot]=i;
            mn[i] = i;
            tt[i] = 1;
            //质数函数值 
        }
        for(int j=1;j<=tot;j++) {
            if(i*pr[j]>maxn) break;
            isprime[i*pr[j]]=0;
            
            mn[i*pr[j]] = pr[j];
                
            if(i%pr[j]==0) {
                break;
            } 
        }
    } 
}
int main()
{
    euler();
    tt[0] = 0;
    for (int i = 1; i < maxn; ++ i) tt[i] += tt[i - 1];            //预处理出小于i的素数的个数 
    scanf("%d", &T);
    while (T --)
    {
        int n, d, ans;
        scanf("%d%d", &n, &d); n --;
        if (d >= maxn)
        {
            ans = n / d;
            for (int i = 1; pr[i] < ans; ++ i)
                if (d % pr[i] == 0) {ans = pr[i]; break;}
            printf("%d\n", tt[ans]);
        }
        else printf("%d\n", tt[min(n / d, mn[d])]);
    }
}

 

转载于:https://www.cnblogs.com/zhangjialu2015/p/5700331.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值