Cod sursa(job #1455883)

Utilizator GeiGeiGeorge Cioroiu GeiGei Data 29 iunie 2015 12:56:30
Problema Kperm Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.92 kb
#include <cstdio>
#include <iostream>
#include <vector>
#include <set>
#include <cmath>
#define mod 666013

using namespace std;

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

int put (int b, int e) {
    int 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);

    int n, k;
    cin >> n >> k;
    int a = (n - 1) % k + 1;
    int b = n / k;
    int ans = 1;
    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;
}