K. Lonely Numbers 【数论】1600

K. Lonely Numbers

图片说明

题意

给定1-n,问里面有多少个孤独的数,如果一个数a,b可以使得a,a/gcd(a,b),b/gcd(a,b)构成三角形,a,b就不会孤独。

解法

对于非平方数的合数x,必定有一个因子是小于,令其为y,则
的gcd为y,三角形就变成了 ,其中, 因为
所以就能构成合法的三角形,y + z-1>=Z 而这里的y*(z-1) 和 z都是合数,所以是合数就不是孤独的

对于质数x,只有他的 出现才可以,构成这样的三角形

#include <bits/stdc++.h>
#define ios ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define debug_in  freopen("in.txt","r",stdin)
#define debug_out freopen("out.txt","w",stdout);
#define pb push_back
#define all(x) x.begin(),x.end()
#define fs first
#define sc second
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
const int maxn = 1e6+10;
const int maxM = 1e6+10;
const int inf = 1e8;
const ll inf2 = 1e17;

template<class T>void read(T &x){
    T s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
    x = s*w;
}
template<class H, class... T> void read(H& h, T&... t) {
    read(h);
    read(t...);
}
template <typename ... T>
void DummyWrapper(T... t){}

template <class T>
T unpacker(const T& t){
    cout<<' '<<t;
    return t;
}
template <typename T, typename... Args>
void pt(const T& t, const Args& ... data){
    cout << t;
    DummyWrapper(unpacker(data)...);
    cout << '\n';
}


//--------------------------------------------
int T,N;
bool vis[maxn];
int P[maxn],tail;
int ans[maxn];
void init(int N){
    for(int i = 2;i<=N;i++){
        if(!vis[i]) P[tail++] = i;
        for(int j = 0;j<tail && P[j] <= N/i;j++){
            vis[P[j] * i] = 1;
            if(i%P[j] == 0) break;
        }
    }

    for(int i = 1;i<=N;i++){ //预处理所有的答案
        ans[i] += ans[i-1];
        if(!vis[i]) ans[i]++;
        int ok = 0;
        int x = sqrt(i);
        if(x * x == i && x != 1){
            if(!vis[i]) ans[i]--,vis[i] = 0;
            if(!vis[x]) ans[i]--,vis[x] = 0;
        }
    }
}

int main(){
    // debug_in;

    read(T);
    init(1000000);
    while(T--){
        read(N);
        printf("%d\n",ans[N]);
    }

    return 0;
}
全部评论

相关推荐

2024-12-29 11:08
湖南工业大学 Java
程序员牛肉:简历没什么大问题了。 而且不要再换项目了。三月份就开暑期实习了,现在都一月份了。实在来不及重新开一下项目了。把一个项目写完或许很快,但是把一个项目搞懂吃透并不简单。所以不要换项目了,把你简历上面的两个项目好好挖一挖吧。 具体 体现在:你能不能流利的说出你的项目的每一个功能点代码实现?你能不能说出在这块除了A技术之外,还有其他技术能够实现嘛?如果有其他技术能够实现,那你这块为什么选择了你当前用的这个技术?
投递牛客等公司10个岗位
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客企业服务