Pagini recente » Cod sursa (job #2607206) | Clasamentul arhivei ACM | Cod sursa (job #217143) | Cod sursa (job #932886) | Cod sursa (job #612550)
Cod sursa(job #612550)
#include<fstream.h>
#define N 666013
#define lul unsigned long long
lul i,n,k,c,r,a,b,d,e,t,s,fact[5000],m;
lul put(lul x,lul n)
{lul i,x1=x%N,x2=(x*x)%N,k=0,v[70];
while(n)
v[++k]=n%2,n/=2;
for(i=k-1;i;i--)
if(!v[i])
x2=(x1*x2)%N,x1=(x1*x1)%N;
else
x1=(x1*x2)%N,x2=(x2*x2)%N;
return x1;}
int main()
{ifstream f("kperm.in");
ofstream g("kperm.out");
f>>n>>k;
if(k%2==0)
g<<"0";
else
{c=n/k,r=n%k;
fact[0]=1;
m=r;
if(k-r>m)
m=k-r;
if(c+1>m)
m=c+1;
for(i=1;i<=m;i++)
fact[i]=(i*fact[i-1])%N;
a=fact[r];
b=fact[k-r];
d=fact[c];
e=fact[c+1];
t=(a*b)%N;
d=put(d,k-r);
e=put(e,r);
s=(d*e)%N;
g<<((t*s)%N);}
return 0;}