Cod sursa(job #282371)

Utilizator alecmanAchim Ioan Alexandru alecman Data 17 martie 2009 16:24:48
Problema Kperm Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.95 kb
#include<stdio.h>

#define INPUT "kperm.in"
#define OUTPUT "kperm.out"
#define LL long long

const long MOD = 666013;

FILE *fin = fopen(INPUT, "r"), *fout = fopen(OUTPUT, "w");

LL N, K, R, C;
LL Sum, P;

int main()
{
  fscanf(fin, "%lld %lld", &N, &K);

  Sum = 1;

  if( K % 2 )
  {
    R = N % K;
    C = N / K;
    P = 1;

    for(LL i = 1; i <= R; ++i)
      P = (P * (i % MOD)) % MOD;

    Sum = (Sum * P) % MOD;
    P = 1;

    for(LL i = 1; i <= K-R; ++i)
      P = (P * (i % MOD)) % MOD;
    
    Sum = (Sum * P) % MOD;
    P = 1;

    for(LL i = 1; i <= C+1; ++i)
      P = (P * (i % MOD)) % MOD;
    for(LL i = 2; i <= R; ++i)
      P = ((P % MOD) * (P % MOD)) % MOD;

    Sum = (Sum * P) % MOD;
    P = 1;

    for(LL i = 1; i <= C; ++i)
      P = (P * (i % MOD)) % MOD;
    for(LL i = 2; i <= K-R; ++i)
      P = ((P % MOD) * (P % MOD)) % MOD;

    Sum = (Sum * P) % MOD;
  }
  
  fprintf(fout, "%lld\n", Sum);

  fclose(fin);
  fclose(fout);

  return 0;
}