Submission #7024642

---

Source Code Expand

#include<bits/stdc++.h>
using namespace std;

long long mod = 1'000'000'007;

long long modpow(long long x, long long y){
    if (y == 0) return 1;
    if (y % 2) return x * modpow(x, y - 1) % mod;
    return modpow(x, y / 2) * modpow(x, y / 2) % mod;
}

long long cmb(long long x, long long y){
    if (x < 0 || y < 0) return 0;
    long long acc = 1;
    long long div = 1;
    for (int i = 0; i < y; i++){
        acc *= x - i;
        acc %= mod;
        div *= i + 1;
        div %= mod;
    }
    return acc * modpow(div, mod - 2) % mod;
}

int main(){
    int r, c, x, y, d, l;
    cin >> r >> c >> x >> y >> d >> l;
    //包除原理をつかう
    long long ans;
    ans = 2 * (cmb((x - 1) * y, d) * cmb((x - 1) * y - d, l) % mod) % mod; ans %= mod;
    ans += (2 * (cmb(x * (y - 1), d) * cmb(x * (y - 1) - d, l) % mod)) % mod; ans %= mod;
    ans += mod - cmb((x - 2) * y, d) * cmb((x - 2) * y - d, l) % mod; ans %= mod;
    ans += mod - cmb(x * (y - 2), d) * cmb(x * (y - 2) - d, l) % mod; ans %= mod;
    ans += mod - cmb((x - 2) * (y - 2), d) * cmb((x - 2) * (y - 2) - d, l) % mod; ans %= mod;
    ans += mod - 4 * cmb((x - 1) * (y - 1), d) * cmb((x - 1) * (y - 1) - d, l) % mod; ans %= mod;
    ans += 2 * cmb((x - 1) * (y - 2), d) * cmb((x - 1) * (y - 2) - d, l) % mod; ans %= mod;
    ans += 2 * cmb((x - 2) * (y - 1), d) * cmb((x - 2) * (y - 1) - d, l) % mod; ans %= mod;
    ans = mod + cmb(x * y, d) * cmb(x * y - d, l) - ans; ans %= mod;
    ans *= r - x + 1; ans %= mod;
    ans *= c - y + 1; ans %= mod;
    cout << ans << endl;
}

Submission Info

Judge Result