Cod sursa(job #1455898)

Utilizator GeiGeiGeorge Cioroiu GeiGei Data 29 iunie 2015 13:23:32
Problema Kperm Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.05 kb
#include <cstdio>
#include <iostream>
#include <vector>
#include <set>
#include <cmath>
#define mod 666013

using namespace std;

long long fact (long long a) {
    long long ans = 1;
    for (long long i = 2; i <= a; i++) {
        ans = (ans * i) % mod;
    }
    return ans;
}

long long put (long long b, long long e) {
    long long ans = 1;
    while (e > 0) {
        if (e % 2 == 1) {
            ans = (ans * b) % mod;
        }
        b = (b * b) % mod;
        e /= 2;
    }
    return ans;
}

int main() {
    freopen("kperm.in", "r", stdin);
    freopen("kperm.out", "w", stdout);

    ios_base::sync_with_stdio(false);
    cin.tie(0);

    long long n, k;
    cin >> n >> k;
    if (k % 2 == 0) {
        cout << 0;
        return 0;
    }
    long long ans = 1;
    long long a = n % k;
    long long b = n / k;
    ans = (ans * fact(a)) % mod;
    ans = (ans * fact(k - a)) % mod;
    ans = (ans * put(fact(b + 1), a)) % mod;
    ans = (ans * put(fact(b), k - a)) % mod;
    cout << ans;

    return 0;
}