Pagini recente » Cod sursa (job #80769) | Cod sursa (job #1544680) | Cod sursa (job #1402829) | Cod sursa (job #2446570) | Cod sursa (job #26238)
Cod sursa(job #26238)
#include <stdio.h>
long long nr(int n, int p)
{
long long rez = 0;
int tmp;
for (tmp = p; n >= tmp; tmp *= p)
rez = rez + (long long)p * (n / tmp) * (n / tmp - 1) / 2 + ((long long)n % tmp + 1) * (n / tmp);
return rez;
}
int main()
{
int n, b, i, step, tmp, a[16], put[16];
long long crt[16], sol;
freopen("zero2.in", "r", stdin);
freopen("zero2.out", "w", stdout);
for (step = 1; step <= 10; ++step)
{
scanf("%d%d", &n, &b);
for (a[0] = 0, tmp = b, i = 2; i * i <= b; ++i)
if (tmp % i == 0)
{
a[++a[0]] = i;
put[a[0]] = 0;
while(tmp % i == 0)
{
tmp /= i;
++put[a[0]];
}
}
if (tmp > 1)
{
a[++a[0]] = tmp;
put[a[0]] = 1;
}
for (i = 1; i <= a[0]; ++i)
crt[i] = nr(n, a[i]);
sol = -1;
for (i = 1; i <= a[0]; ++i)
if (sol == -1 || sol > crt[i] / put[i])
sol = crt[i] / put[i];
printf("%d\n", sol);
}
return 0;
}