結果

問題 No.1253 雀見椪
ユーザー firiexpfiriexp
提出日時 2020-10-09 21:47:04
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 33 ms / 2,000 ms
コード長 2,557 bytes
コンパイル時間 863 ms
コンパイル使用メモリ 102,672 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-27 16:13:03
合計ジャッジ時間 2,217 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 4 ms
4,380 KB
testcase_04 AC 32 ms
4,376 KB
testcase_05 AC 31 ms
4,380 KB
testcase_06 AC 31 ms
4,380 KB
testcase_07 AC 32 ms
4,380 KB
testcase_08 AC 33 ms
4,376 KB
testcase_09 AC 32 ms
4,376 KB
testcase_10 AC 32 ms
4,376 KB
testcase_11 AC 32 ms
4,380 KB
testcase_12 AC 31 ms
4,376 KB
testcase_13 AC 32 ms
4,376 KB
testcase_14 AC 1 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <numeric>
#include <bitset>
#include <cmath>

static const int MOD = 1000000007;
using ll = long long;
using u32 = unsigned;
using u64 = unsigned long long;
using namespace std;

template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;

template <u32 M>
struct modint {
    u32 val;
public:
    static modint raw(int v) { modint x; x.val = v; return x; }
    modint() : val(0) {}
    template <class T>
    modint(T v) { ll x = (ll)(v%(ll)(M)); if (x < 0) x += M; val = u32(x); }
    modint(bool v) { val = ((unsigned int)(v) % M); }
    modint& operator++() { val++; if (val == M) val = 0; return *this; }
    modint& operator--() { if (val == 0) val = M; val--; return *this; }
    modint operator++(int) { modint result = *this; ++*this; return result; }
    modint operator--(int) { modint result = *this; --*this; return result; }
    modint& operator+=(const modint& b) { val += b.val; if (val >= M) val -= M; return *this; }
    modint& operator-=(const modint& b) { val -= b.val; if (val >= M) val += M; return *this; }
    modint& operator*=(const modint& b) { u64 z = val; z *= b.val; val = (u32)(z % M); return *this; }
    modint& operator/=(const modint& b) { return *this = *this * b.inv(); }
    modint operator+() const { return *this; }
    modint operator-() const { return modint() - *this; }
    modint pow(long long n) const { modint x = *this, r = 1; while (n) { if (n & 1) r *= x; x *= x; n >>= 1; } return r; }
    modint inv() const { return pow(M-2); }
    friend modint operator+(const modint& a, const modint& b) { return modint(a) += b; }
    friend modint operator-(const modint& a, const modint& b) { return modint(a) -= b; }
    friend modint operator*(const modint& a, const modint& b) { return modint(a) *= b; }
    friend modint operator/(const modint& a, const modint& b) { return modint(a) /= b; }
    friend bool operator==(const modint& a, const modint& b) { return a.val == b.val; }
    friend bool operator!=(const modint& a, const modint& b) { return a.val != b.val; }
};
using mint = modint<MOD>;

void solve(){
    ll n; int a, b, c, d, e, f;
    scanf("%lld %d %d %d %d %d %d", &n, &a, &b, &c, &d, &e, &f);
    mint p = mint(a)/mint(b), q = mint(c)/mint(d), r = mint(e)/mint(f);
    printf("%d\n", (1-(p+q).pow(n)-(q+r).pow(n)-(r+p).pow(n) + p.pow(n)*2 + q.pow(n)*2 + r.pow(n)*2).val);
}

int main() {
    int t;
    cin >> t;
    while(t--){
        solve();
    }
    return 0;
}
0