Cod sursa(job #72550)

Utilizator info_arrandrei gigea info_arr Data 14 iulie 2007 12:21:04
Problema Kperm Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.09 kb
using namespace std;

#define MOD 666013
#define ll long long

#include <stdio.h>

FILE *fin=fopen("kperm.in","r"),
     *fout=fopen("kperm.out","w");

int n,k,r,c;
ll sol;

int main()
{
  fscanf(fin,"%lld %lld",&n,&k);
  if (k%2==0)
   {
      fprintf(fout,"0\n");
      return 0;
   }
  else 
   {
      int i;
      ll fact=1; sol=1;
      c=n/k; r=n%k;
      for (i=2; i<=r; i++)
       fact=(fact*i)%MOD;
      sol=(long long)(sol*fact)%MOD;
      fact=1;
      for (i=2; i<=k-r; i++)
       fact=(fact*i)%MOD;
      sol=(long long)(sol*fact)%MOD;
      fact=1;
      for (i=2; i<=c+1; i++)
       fact=(fact*i)%MOD;
      ll p=fact; 
      for (i=2; i<=r; i++)  
       fact=(fact*p)%MOD;
      if (r==0) fact=1;  
      sol=(long long)(sol*fact)%MOD;
      fact=1; 
      for (i=1; i<=c; i++)
       fact=(fact*i)%MOD;
      p=fact; 
      for (i=1; i<=k-r; i++)
       fact=(fact*p)%MOD;
  //    if (k-r==0) fact=1; 
      sol=(long long)(sol*fact)%MOD;   
      sol%=MOD;
      fprintf(fout,"%lld\n",sol);
    }
return 0;
}