数论--素数

#include <stdio.h>
#include <math.h>

// Eratosthenes筛选法
void sieveofe(int p[], int n)
{
    int i, j;

    p[0] = 0;
    p[1] = 0;
    p[2] = 1;

    // 初始化
    for(i=3; i<=n; i++) {
        p[i++] = 1;
        p[i] = 0;
    }
    int max = sqrt(n);
    for(i=3; i<=max; i++){
        if(p[i]) {
            for(j=i+i; j < n; j+=i)    //进行筛选
                p[j]=0;
        }
    }
}

#define MAX 10000

int main(void)
{
    int p[MAX+1];  //放素数的数组
    int sum, count, i;

    sieveofe(p, MAX);

    for(;;) {
        scanf("%d", &sum);
        if(sum == 0)
            break;

        count = 0;
        for(i=2; i<=sum / 2; i++) {
            if(p[i] && p[sum-i])
                if(i != sum-i) //两个数不相同
                    count++;
        }

        printf("%d\n", count);
    }

    return 0;
}










/* UVA583 UVALive5406 Prime Factors */

#include <bits/stdc++.h>

using namespace std;

const int N31 = 2147483647;           // 2 ^ 31
const int N = ceil(sqrt((double) N31));
const int SQRTN = ceil(sqrt((double) N));
bool isPrime[N + 1];
int prime[N + 1], pcount;

// Eratosthenes筛选法
void esieve(void)
{
    memset(isPrime, true, sizeof(isPrime)); //做标记,假设都是素数

    isPrime[0] = isPrime[1] = false;  //0和1不是素数,先标记好
    pcount = 0;                        //计数工具
    for(int i=2; i<=N; i++) {  //判断是不是素数,从2开始
        if(isPrime[i]) {
            prime[pcount++] = i;

            for(int j=i*i; j<=N; j+=i)  //筛选
                isPrime[j] = false;
        }
    }
}
  //true是正的,false是负的
int main()
{
    esieve();

    int n;
    while(~scanf("%d", &n) && n) {
        bool flag = true;

        printf("%d =", n);
        if(n < 0) {
            printf(" -1");
            n = -n;
            flag = false;
        }

        for(int i=0; i<pcount; i++) {  //所有的都试一遍
            while(n % prime[i] == 0 && n != 1) {
                if(flag) {
                    printf(" %d", prime[i]);
                    flag = false;
                } else
                    printf(" x %d", prime[i]);
                n /= prime[i];
            }
        }

        if(n != 1) {
            if(flag)
                printf(" %d", n);
            else
                printf(" x %d", n);
        }

        printf("\n");
    }

    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值