[ACM实例代码] 01背包 完全背包 快速幂 筛法求素数

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

const int max_n=2300;
const int max_m=5000;
const int primen=1000000;

bool prime[primen];


typedef long long LL;

int max(int a,int b)
{
    return a>b?a:b;
}

void package_all(void)                  //完全背包
{
    int w[max_n]={7,4,10,2,9,8};
    int v[max_n]={9,9,20,5,7,6};
    int dp[max_m];
    int N,V;
    N=6;
    V=200;
    memset(dp,0,sizeof(dp));
    for(int i=0;i<N;i++)
        for(int vi=v[i];vi<=V;vi++)
            dp[vi]=max(dp[vi],dp[vi-v[i]]+w[i]);
    for(int i=0;i<=V;i++) printf("%d
",dp[i]);
}

void package01(void)                    //01背包
{
    int w[max_n]={7,4,10,2,9,8};
    int v[max_n]={9,9,20,5,7,6};
    int dp[max_m];
    int N,V;
    N=6;
    V=200;
    memset(dp,0,sizeof(dp));
    for(int i=0;i<N;i++)
        for(int vi=V;vi>=v[i];vi--)
            dp[vi]=max(dp[vi],dp[vi-v[i]]+w[i]);
    for(int i=0;i<=V;i++) printf("%d
",dp[i]);
}

LL quick_pow(LL a,LL b,LL c) //a^b mod c    快速幂
{
    LL res=1;
    printf("%lld^%lld mod %lld is ",a,b,c);
    while(b)
    {
        if(b&1)
        {
            res=res*a%c;
        }
        a=a*a%c;
        b>>=1;
    }
    printf("%lld
",res);
    return res;
}

void make_prime(void)                   //经典筛法求素数
{
    prime[0]=prime[1]=0;
    prime[2]=1;
    for(int i=3;i<=primen;i++)
        prime[i]=i%2;
    for(int i=3;i<=(int)sqrt(primen*1.0);i++)
    {
        if(prime[i])
        {
            for(int j=i+i;j<=primen;j+=i)
                prime[j]=0;
        }
    }
    for(int i=0;i<=1000;i++)
        if(prime[i]) printf("%d ",i);
}

int main(void)
{
    make_prime();
    printf("
");
    quick_pow(2,2147482,2147481);
    package_all();
    package01();
}