結果
問題 | No.1181 Product Sum for All Subsets |
ユーザー |
![]() |
提出日時 | 2022-01-11 12:39:00 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,721 bytes |
コンパイル時間 | 775 ms |
コンパイル使用メモリ | 84,316 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-14 11:37:05 |
合計ジャッジ時間 | 1,931 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include <iostream>#include <cstdio>#include <string>#include <algorithm>#include <utility>#include <cmath>#include <vector>#include <stack>#include <queue>#include <deque>#include <set>#include <map>#include <tuple>#include <numeric>#include <functional>using namespace std;typedef long long ll;typedef vector<ll> vl;typedef vector<vector<ll>> vvl;typedef pair<ll, ll> P;#define rep(i, n) for(ll i = 0; i < n; i++)#define exrep(i, a, b) for(ll i = a; i <= b; i++)#define out(x) cout << x << endl#define exout(x) printf("%.10f\n", x)#define chmax(x, y) x = max(x, y)#define chmin(x, y) x = min(x, y)#define all(a) a.begin(), a.end()#define rall(a) a.rbegin(), a.rend()#define pb push_back#define re0 return 0const ll mod = 1000000007;const ll INF = 1e16;const ll MAX = 10;ll fact[MAX]; // fact[i] : iの階乗のmodll inv[MAX]; // inv[i] : iの逆数のmodll invfact[MAX]; // invfact[i] : iの階乗の逆数のmodvoid init() {fact[0] = 1;inv[0] = inv[1] = 1;invfact[0] = 1;for(ll i = 1; i < MAX; i++) {fact[i] = i * fact[i-1] % mod;if(i >= 2) {inv[i] = mod - inv[mod % i] * (mod / i) % mod;}invfact[i] = invfact[i-1] * inv[i] % mod;}}// a^n (mod.MOD)を求める。計算量はO(logn)ll modpow(ll a, ll n, ll MOD = mod) {if(n == 0) {return 1;}if(n%2 == 1) {return (a * modpow(a, n-1, MOD)) % MOD;}return (modpow(a, n/2, MOD) * modpow(a, n/2, MOD)) % MOD;}int main() {ll n, k;cin >> n >> k;k %= mod;init();out((modpow(k * (k + 3) % mod * inv[2] % mod, n) + mod - modpow(k * (k + 1) % mod * inv[2] % mod, n)) % mod);re0;}