結果

問題 No.1181 Product Sum for All Subsets
ユーザー moririn2528_cmoririn2528_c
提出日時 2020-08-21 21:51:54
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 76 ms / 2,000 ms
コード長 1,924 bytes
コンパイル時間 1,294 ms
コンパイル使用メモリ 93,712 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-05 07:59:30
合計ジャッジ時間 2,861 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,384 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 1 ms
4,380 KB
testcase_11 AC 2 ms
4,380 KB
testcase_12 AC 63 ms
4,380 KB
testcase_13 AC 15 ms
4,376 KB
testcase_14 AC 60 ms
4,376 KB
testcase_15 AC 49 ms
4,380 KB
testcase_16 AC 66 ms
4,380 KB
testcase_17 AC 29 ms
4,380 KB
testcase_18 AC 59 ms
4,380 KB
testcase_19 AC 54 ms
4,380 KB
testcase_20 AC 65 ms
4,376 KB
testcase_21 AC 36 ms
4,380 KB
testcase_22 AC 22 ms
4,380 KB
testcase_23 AC 10 ms
4,376 KB
testcase_24 AC 20 ms
4,380 KB
testcase_25 AC 34 ms
4,376 KB
testcase_26 AC 4 ms
4,380 KB
testcase_27 AC 1 ms
4,380 KB
testcase_28 AC 75 ms
4,380 KB
testcase_29 AC 76 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#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;
}
0