結果

問題 No.3495 ２変数半二項展開
コンテスト
ユーザー askr58
提出日時 2026-04-03 23:14:52
言語 C++23
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
RE  
実行時間 -
コード長 1,075 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,718 ms
コンパイル使用メモリ 345,944 KB
実行使用メモリ 7,972 KB
最終ジャッジ日時 2026-04-03 23:15:08
合計ジャッジ時間 5,589 ms
ジャッジサーバーID
（参考情報） 		judge4_0 / judge5_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 7 RE * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
using ll=long long;
#include <atcoder/modint>
using mint=atcoder::modint;

int main(){
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	int ttt;
	cin>>ttt;
	while(ttt--){
		ll n,m,l,b;
		cin>>n>>m>>l>>b;
		mint::set_mod(b);
		mint x=m,y=l;
		vector<vector<mint>> a(4,vector<mint>(4));
		a[0][1]=2*y;
		a[0][2]=x-y;
		a[1][0]=2;
		a[1][3]=x-y;
		a[2][0]=1;
		a[3][1]=1;
		vector<vector<mint>> res(4,vector<mint>(4));
		for(int i=0;i<4;i++)res[i][i]=1;
		ll N=n;
		while(N>0){
			if(N&1){
				vector<vector<mint>> nres(4,vector<mint>(4));
				for(int i=0;i<4;i++)for(int j=0;j<4;j++)for(int k=0;k<4;k++){
					nres[i][j]+=a[i][k]*res[k][j];
				}
				res=move(nres);
			}
			vector<vector<mint>> na(4,vector<mint>(4));
			
			for(int i=0;i<4;i++)for(int j=0;j<4;j++)for(int k=0;k<4;k++){
				na[i][j]+=a[i][k]*a[k][j];
			}
			a=move(na);
			N>>=1;
		}
		mint p=2*res[2][1]+2*res[2][2],q=2*res[3][1]+2*res[3][2];
		//cout<<p.val()<<" "<<q.val()<<endl;
		if(n%2==0)cout<<(p/2).val()<<endl;
		else cout<<(q/2).val()<<endl;
	}
}
0