Pagini recente » Cod sursa (job #667799) | Cod sursa (job #1025114) | Cod sursa (job #1246269) | Cod sursa (job #1009645) | Cod sursa (job #1993192)
#include <fstream>
using namespace std;
ifstream fin("kperm.in");
ofstream fout("kperm.out");
int n,k,c,r;
long long f[5010];
long long pw(int a, int b)
{
if (b == 0)
return 1;
long long nr = pw(a, b / 2);
if (b % 2 == 0)
return (nr * nr) % 666013;
return (a * nr * nr) % 666013;
}
int main()
{
fin >> n >> k;
c = n / k;
r = n % k;
if (k % 2 == 0)
{
fout << 0 << '\n';
return 0;
}
f[0] = 1;
for (int i = 1;i <= n;++i)
f[i] = (f[i - 1] * i) % 666013;
long long x = f[r];
x = (x * f[k - r]) % 666013;
x = (x* pw(f[c + 1], r)) % 666013;
x = (x * pw(f[c],k - r)) % 666013;
fout << x << '\n';
return 0;
}