10亿以内和987654互质正整数的和[转自深圳-冒泡]

容斥原理，把与987654不互质的数的和算出来，即可。

#include <iostream>
#include <cstdio>
using namespace std;
long long pri[4]={2,3,97,1697};         //987654的质因子
long long big=1000000000;
int main()
{
    long long sum=0,mat;
    int i,j;
    for(i=1;i<16;i++){
        int cnt=0;
        int temp=i;
        mat=1;
        for(j=0;temp;temp>>=1,j++){
            if(temp%2){
                mat*=pri[j];
                cnt++;
            }
        }
        if(cnt%2)
            sum+=(mat+big/mat*mat)*(big/mat)/2;
        else
            sum-=(mat+big/mat*mat)*(big/mat)/2;
    }
    //printf("%lld
",sum);
    printf("%lld
",(1+big)*big/2-sum);
    return 0;
}


/*

//以下代码计算987654的质因子
#include <iostream>
#include <cstdio>
using namespace std;

int main()
{
    int x=987654,i;
    for(i=2;i<x;i++){
        if(x%i==0){
            printf("%d
",i);
            while(x%i==0)
                x/=i;
        }
    }
    if(x>1)printf("%d
",x);
    return 0;
}
*/