• 求最大公约数和小于n的所有质数


    //algorithm.h
    enum SWAP_TYPE{MEMORY, COMPLEX};
    struct SIntArray
    {
        int *pData;
        int  num;
        SIntArray():pData(NULL),num(0){}
        void Clear(){delete pData; pData = NULL; num = 0;}
    };
    void      Wswap(int &m, int &n, SWAP_TYPE name = MEMORY);
    int       Wgcd(int m, int n);
    SIntArray Wprime(int m);
    
    //algorithm.cpp
    void Wswap(int &m, int &n, SWAP_TYPE name)//交换
    {
        switch(name)
        {
        case MEMORY:
            m = m + n;
            n = m - n;
            m = m - n;
            break;
        case COMPLEX:
            int temp = m;
            m = n;
            n = temp;
            break;
        }
    }
    
    int GCD(int m, int n)
    {
    
        if( n == 0)
        {
            return m;
        }
        else
        {
            int r = m % n;
            GCD(n, r);
        }
    }
    
    int Wgcd(int m, int n)//寻找最大公约数
    {
        if(m > n)
        {
            Wswap(m, n);//make sure m > n    
        }
        int result = GCD(m, n);
        return result;
    }
    
    SIntArray Wprime(int m)//找所有小于m的质数
    {
        SIntArray result;
        if(m <= 0)
        {
            return result;
        }
        int *number = new int [m+1];
        for(int i = 0; i <= m; i++)
        {
            number[i] = i;
        }
    
        int limit = sqrt((float)m);
        for(int j = 2; j <= limit; j++)
        {
            if(number[j] != 0)
            {
                int temp = j*j;// j 2*j   (j-1)*j    1:(j-1) have been processed, so only process j ... n/j
    while(temp <= m) // eliminate { number[temp] = 0; temp = temp + j; } } } int count = 0; for(int j = 2; j <= m; j++) { if(number[j]) { number[count] = number[j]; count++; } } int *out = new int [count]; memcpy(out, number, count*sizeof(int)); delete number; number = NULL; result.pData = out; result.num = count; return result; }
    //main.cpp
    int _tmain(int argc, _TCHAR* argv[]) { /* Wgcd 最大公约数 int m,n; cout<<"input two number:"<<endl; cin>>m>>n; cout<<"greatest common divisor:"<<endl; cout<<Wgcd(m, n)<<endl; */ int m; //小于n的所有质数 cout<<"input one number:"<<endl; cin>>m; cout<<"prime number:"<<endl; SIntArray result; result = Wprime(m); for(int i = 0; i <result.num; i++) { cout<<result.pData[i]<<" "; } result.Clear(); return 0; }
  • 相关阅读:
    有点难度的二分
    请教神牛_字符串hash
    引水进城
    dp的斜率优化
    关于学习oi的一些事项
    永续债
    消费税
    增值税
    BSC交流
    钉钉吐槽功能点
  • 原文地址:https://www.cnblogs.com/welen/p/3856144.html
Copyright © 2020-2023  润新知