Pagini recente » Cod sursa (job #246600) | Rating Fejer Magdolna (magdifej) | Cod sursa (job #1652571) | Cod sursa (job #1645403) | Cod sursa (job #1472666)
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxV 1000
#define maxL 10000
#define base 1000000000
#define ct 9
using namespace std;
int n, x, i, j, k, val, w[maxL];
int dp[2][maxV + 2][maxL];
void add(int A[], int B[])
{
int i, t = 0;
for (i = 1; i <= A[0] || i <= B[0] || t; i++, t /= base)
A[i] = (t += A[i] + B[i]) % base;
A[0] = i - 1;
}
int gcd(int d, int i)
{
int r;
if (d < i)
swap(d, i);
r = d % i;
while (r)
d = i,
i = r,
r = d % i;
return i;
}
void read()
{
freopen("indep.in", "r", stdin);
freopen("indep.out", "w", stdout);
scanf("%d", &n);
w[++ w[0]] = 1;
for (i = 1; i <= n; ++ i)
{
scanf("%d", &val);
add(dp[i & 1][val], w);
for (j = 1; j <= maxV; ++ j)
{
k = gcd(j, val);
add(dp[i & 1][k], dp[(i - 1) & 1][j]);
add(dp[i & 1][j], dp[(i - 1) & 1][j]);
}
memset(dp[(i - 1) & 1], 0, sizeof(dp[(i - 1) & 1]));
}
}
void write()
{
int x, Pow;
freopen("indep.out", "w", stdout);
n &= 1;
printf("%d", dp[n][1][dp[n][1][0]]);
for (i = dp[n][1][0] - 1; i >= 1; -- i)
{
x = 1; Pow = 0;
while (dp[n][1][i] > x)
x *= 10,
++ Pow;
for (j = Pow + 1; j <= ct; ++ j)
printf("0");
printf("%d", dp[n][1][i]);
}
}
int main()
{
read();
write();
return 0;
}