#700643 (C++17) No.1684 Find Brackets

• 問題: No.1684 Find Brackets
• ユーザー: 沙耶花
• 提出日時: 2021-09-17 23:00:04
• 言語: C++17 (gcc 13.2.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 124 ms / 2,000 ms
• コード長: 1,346 bytes
• コンパイル時間: 4,625 ms
• コンパイル使用メモリ: 265,620 KB
• 実行使用メモリ: 48,040 KB
• 最終ジャッジ日時: 2023-09-12 09:03:36
• 合計ジャッジ時間: 8,028 ms
• ジャッジサーバーID （参考情報）: judge14 / judge13

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	99 ms; 43,804 KB
testcase_01	AC	99 ms; 43,860 KB
testcase_02	AC	124 ms; 47,696 KB
testcase_03	AC	99 ms; 43,888 KB
testcase_04	AC	99 ms; 43,956 KB
testcase_05	AC	99 ms; 43,852 KB
testcase_06	AC	99 ms; 44,016 KB
testcase_07	AC	114 ms; 47,848 KB
testcase_08	AC	114 ms; 47,324 KB
testcase_09	AC	105 ms; 46,772 KB
testcase_10	AC	101 ms; 45,052 KB
testcase_11	AC	102 ms; 44,376 KB
testcase_12	AC	110 ms; 46,128 KB
testcase_13	AC	115 ms; 47,228 KB
testcase_14	AC	122 ms; 47,704 KB
testcase_15	AC	107 ms; 47,732 KB
testcase_16	AC	122 ms; 48,040 KB
testcase_17	AC	123 ms; 47,936 KB
testcase_18	AC	124 ms; 47,896 KB
testcase_19	AC	99 ms; 44,016 KB
testcase_20	AC	99 ms; 44,020 KB
testcase_21	AC	106 ms; 47,752 KB
testcase_22	AC	107 ms; 47,908 KB
testcase_23	AC	107 ms; 47,700 KB

ソースコード

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint1000000007;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000001

struct combi{
	deque<mint> kaijou;
	deque<mint> kaijou_;
	
	combi(int n){
		kaijou.push_back(1);
		for(int i=1;i<=n;i++){
			kaijou.push_back(kaijou[i-1]*i);
		}
		
		mint b=kaijou[n].inv();
		
		kaijou_.push_front(b);
		for(int i=1;i<=n;i++){
			int k=n+1-i;
			kaijou_.push_front(kaijou_[0]*k);
		}
	}
	
	mint combination(int n,int r){
		if(r>n)return 0;
		mint a = kaijou[n]*kaijou_[r];
		a *= kaijou_[n-r];
		return a;
	}
	
	mint junretsu(int a,int b){
		mint x = kaijou_[a]*kaijou_[b];
		x *= kaijou[a+b];
		return x;
	}
	
	mint catalan(int n){
		return combination(2*n,n)/(n+1);
	}
	
};
 
int main(){
	
	int N,M;
	cin>>N>>M;
	
	combi C(5000000);
	
	vector<mint> calc(N+5,0);
	rep(i,N+5){
		calc[i] = C.combination(N,i);
		if(i!=0)calc[i] += calc[i-1];
	}
	
	mint ans = 0;
	mint t = 0;
	for(int i=M;i<=N;i++){
		t = C.combination(N,i);
		t *= i;
		ans += t;
		//mint temp = calc[N]-calc[i];
		//cout<<ans.val()<<',';
		if(N-i > i){
			ans += C.combination(N,i) * abs(N-i - i);
		}
		//cout<<ans.val()<<',';
		ans += calc[N] - calc[max(i,N-i)];
		//cout<<ans.val()<<endl;
	}
	ans *= 2;
	
	cout<<ans.val()<<endl;

	
	return 0;
}