Pagini recente » Cod sursa (job #3306311) | Cod sursa (job #3321705) | Cod sursa (job #666335) | Cod sursa (job #3356693) | Cod sursa (job #1014866)
#include<stdio.h>
#include<string.h>
#define NMAX 1007
#define Baza 10000000
int n, D[2][NMAX][150], h[10], GCD[NMAX][NMAX];
inline int cmmdc(int a, int b){
if(!b)
return a;
return cmmdc(b, a % b);
}
void Adun(int a[], int b[]){
if(a[0] < b[0])
a[0] = b[0];
int t = 0;
for(int i = 1; i <= a[0]; ++ i){
a[i] += t + b[i];
t = a[i] / Baza;
a[i] %= Baza;
}
while(t != 0){
a[++ a[0]] = t % 10;
t /= 10;
}
}
int main(){
freopen("indep.in", "r", stdin);
freopen("indep.out", "w", stdout);
scanf("%d", &n);
for(int i = 1; i <= 1000; ++ i)
for(int j = 1; j <= 1000; ++ j)
GCD[i][j] = cmmdc(i, j);
h[0] = h[1] = 1;
int a;
for(int i = 1; i <= n; ++ i){
scanf("%d", &a);
Adun(D[1][a], h);
for(int j = 1; j <= 1000; ++ j){
Adun(D[1][j], D[0][j]);
Adun(D[1][GCD[j][a]], D[0][j]);
///D[i + 1][j] += D[i][j];
///D[i + 1][cmmdc(j, a[i])] += D[i][j];
}
memcpy(D[0], D[1], sizeof(D[0]));
memset(D[1], 0, sizeof(D[1]));
}
printf("%d", D[0][1][D[0][1][0]]);
for(int i = D[0][1][0] - 1; i >= 1; -- i)
printf("%07d", D[0][1][i]);
return 0;
}