結果
問題 | No.1181 Product Sum for All Subsets |
ユーザー |
![]() |
提出日時 | 2020-08-21 21:51:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 76 ms / 2,000 ms |
コード長 | 1,924 bytes |
コンパイル時間 | 855 ms |
コンパイル使用メモリ | 95,052 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-15 05:33:21 |
合計ジャッジ時間 | 2,527 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include<iostream> #include<cstdio> #include<cstring> #include<string> #include<vector> #include<cmath> #include<algorithm> #include<map> #include<queue> #include<deque> #include<iomanip> #include<tuple> #include<cassert> #include<set> using namespace std; typedef long long int LL; typedef pair<int,int> P; typedef pair<LL,int> LP; const int INF=1<<30; const LL MAX=1e9+7; void array_show(int *array,int array_n,char middle=' '){ for(int i=0;i<array_n;i++)printf("%d%c",array[i],(i!=array_n-1?middle:'\n')); } void array_show(LL *array,int array_n,char middle=' '){ for(int i=0;i<array_n;i++)printf("%lld%c",array[i],(i!=array_n-1?middle:'\n')); } void array_show(vector<int> &vec_s,int vec_n=-1,char middle=' '){ if(vec_n==-1)vec_n=vec_s.size(); for(int i=0;i<vec_n;i++)printf("%d%c",vec_s[i],(i!=vec_n-1?middle:'\n')); } void array_show(vector<LL> &vec_s,int vec_n=-1,char middle=' '){ if(vec_n==-1)vec_n=vec_s.size(); for(int i=0;i<vec_n;i++)printf("%lld%c",vec_s[i],(i!=vec_n-1?middle:'\n')); } long long int pow_mod(long long int p_a,long long int p_n,long long int p_p=1e9+7){ //p_a^p_n mod p_p long long int p_b=1,p_t=1; for(;p_b<=p_n;p_b<<=1); for(p_b>>=1;p_b>0;p_b>>=1){ p_t*=p_t; if(p_t>=p_p)p_t%=p_p; if(p_n&p_b)p_t*=p_a; if(p_t>=p_p)p_t%=p_p; } return p_t; } long long int gcd(long long int g_a,long long int g_b){ if(g_a<g_b)gcd(g_b,g_a); if(g_b==0)return g_a; return gcd(g_b,g_a%g_b); } long long int divide(long long int d_a,long long int d_b,long long int d_p=1e9+7){ // d_a/d_b mod d_p // prime:d_p is prime if(d_a>=d_p)d_a%=d_p; if(d_a<0)d_a+=d_p; if(d_b>=d_p)d_b%=d_p; if(d_b<0)d_b+=d_p; d_a*=pow_mod(d_b,d_p-2,d_p); return d_a%d_p; } int main(){ LL n,m; int i,j,k; LL a,b,c; LL t=0; LL s=0; cin>>n>>m; m%=MAX; t=divide((1+m)*m,2); a=1,b=pow_mod(m,n); for(i=0;i<n;i++){ s+=a*b; s%=MAX; b=divide(b*t,m); a*=n-i; a=divide(a,i+1); } cout<<s<<endl; }