Pagini recente » Cod sursa (job #1311481) | Cod sursa (job #416106) | Cod sursa (job #385781) | Cod sursa (job #1328573) | Cod sursa (job #612548)
Cod sursa(job #612548)
#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];
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;
fact[0]=1;
for(i=1;i<=n;i++)
fact[i]=(i*fact[i-1])%N;
if(k%2==0)
g<<"0";
else
{c=n/k,r=n%k;
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;}