hdu 5212 Code (容斥)

题意：

int calc
{
  
  int res=0;
  
  for(int i=1;i<=n;i++)
    
    for(int j=1;j<=n;j++)
    
    {
      
      res+=gcd(a[i],a[j])*(gcd(a[i],a[j])-1);
      
      res%=10007;
    
    }
  
  return res;

}

这样的程序的结果。

题解：

直接暴力肯定挂，这题用容斥搞，先筛出所有数的约数，并且对每个约数累加能被多少个数分解得到。然后从大到小枚举最大公约数，大的数如果是小的数的倍数，那么大的数对应的gcd累加和肯定能满足小的数，于是小的数要将去这些重复的。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<stack>
#include<string>
#define B(x) (1<<(x))
using namespace std;
typedef long long ll;
void cmax(int &a,int b){ if(b>a)a=b; }
void cmin(int &a,int b){ if(b<a)a=b; }
void cmax(ll &a,ll b){ if(b>a)a=b; }
void cmin(ll &a,ll b){ if(b<a)a=b; }
void add(int &a,int b,int mod){ a=(a+b)%mod; }
void add(ll &a,ll b,ll mod){ a=(a+b)%mod; }
void add(int &a,int b){ a+=b; }
void add(ll &a,ll b){ a+=b; }
const int oo=0x3f3f3f3f;
const ll MOD=1000000007;
const int maxn = 10005;
int a[maxn];
int id[maxn];

int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
        ll ans=0;
        memset(id,0x3f,sizeof id);
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            id[a[i]]=i;
        }
        for(int i=1;i<=n;i++){
            int min_id=n*2;
             for(int j=a[i]*2;j<=10000;j+=a[i]){
                if(id[j]>i&&id[j]<min_id){
                    min_id=id[j];
                }
             }
            if(min_id!=n*2){
                ans+=min_id;
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}