题意:
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;
}