結果
問題 | No.1181 Product Sum for All Subsets |
ユーザー |
![]() |
提出日時 | 2020-07-24 05:37:27 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,894 bytes |
コンパイル時間 | 2,051 ms |
コンパイル使用メモリ | 195,844 KB |
最終ジャッジ日時 | 2025-01-12 04:53:48 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include <bits/stdc++.h>/*#include <boost/multiprecision/cpp_dec_float.hpp>#include <boost/multiprecision/cpp_int.hpp>*/#pragma GCC target("avx")#pragma GCC optimize("O3")#pragma GCC optimize("unroll-loops")using namespace std;/*namespace mp = boost::multiprecision;using Bint = mp::cpp_int;using Real = mp::number<mp::cpp_dec_float<1024>>;*/#define lli long long int#define uli unsigned long long int#define INF 999999999999#define rep(i,m,n) for(lli i = m;i < n;i++)#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)#define pb(n) push_back(n)#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());#define Sort(n) sort(n.begin(), n.end())#define Rev(n) reverse(n.begin(),n.end())#define Out(S) cout << S << endl#define NeOut(S) cout << S#define HpOut(S) cout << setprecision(50) << S << endl#define Vec(K,L,N,S) vector<L> K(N,S)#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))#define pint pair<lli,lli>#define paf(L,R) pair<L,R>#define mod 1000000007#define MAX 10000000#define ALL(a) a.begin(),a.end()#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))long long N,K,sum=0,num=0;lli modpow(lli n,lli k,lli m){if(k==0)return 1;else if(k%2==1)return modpow(n,k-1,m)*n%m;else{lli temp=modpow(n,k/2,m);return temp*temp%m;}}void dfs(int n,vector<int>&P){if(n)rep(i,0,K){P[n-1]=i+1;dfs(n-1,P);}else{rep(p,0,(1<<N)-1){num=1;rep(q,0,N)if(p&(1<<(N-q-1)))num=num*P[q]%mod;sum=(sum+num)%mod;}}}int main(){ios::sync_with_stdio(false);cin.tie(0);cin >> N >> K;/*Vec(P,int,N,0);dfs(N,P);Out(sum);*/K%=mod;Out(((modpow(K*(K+3)%mod,N,mod)-modpow(K*(K+1)%mod,N,mod))*modpow(modpow(2,mod-2,mod),N,mod)%mod+10LL*mod)%mod);}