結果

問題 No.1548 [Cherry 2nd Tune B] 貴方と私とサイクルとモーメント
ユーザー 沙耶花沙耶花
提出日時 2021-06-11 23:11:12
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 130 ms / 4,500 ms
コード長 2,069 bytes
コンパイル時間 4,830 ms
コンパイル使用メモリ 271,336 KB
実行使用メモリ 18,688 KB
最終ジャッジ日時 2024-12-15 02:17:36
合計ジャッジ時間 10,839 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,820 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 101 ms
10,788 KB
testcase_03 AC 32 ms
7,168 KB
testcase_04 AC 93 ms
17,700 KB
testcase_05 AC 10 ms
6,816 KB
testcase_06 AC 72 ms
10,504 KB
testcase_07 AC 78 ms
18,584 KB
testcase_08 AC 108 ms
11,208 KB
testcase_09 AC 70 ms
7,552 KB
testcase_10 AC 62 ms
17,668 KB
testcase_11 AC 75 ms
17,772 KB
testcase_12 AC 49 ms
6,816 KB
testcase_13 AC 44 ms
17,756 KB
testcase_14 AC 74 ms
6,912 KB
testcase_15 AC 47 ms
18,572 KB
testcase_16 AC 83 ms
11,356 KB
testcase_17 AC 100 ms
17,540 KB
testcase_18 AC 79 ms
18,356 KB
testcase_19 AC 38 ms
7,168 KB
testcase_20 AC 58 ms
17,968 KB
testcase_21 AC 65 ms
10,544 KB
testcase_22 AC 127 ms
18,620 KB
testcase_23 AC 127 ms
18,636 KB
testcase_24 AC 128 ms
18,600 KB
testcase_25 AC 130 ms
18,508 KB
testcase_26 AC 128 ms
18,672 KB
testcase_27 AC 124 ms
18,592 KB
testcase_28 AC 125 ms
18,624 KB
testcase_29 AC 128 ms
18,656 KB
testcase_30 AC 125 ms
18,688 KB
testcase_31 AC 129 ms
18,620 KB
testcase_32 AC 93 ms
18,676 KB
testcase_33 AC 96 ms
18,688 KB
testcase_34 AC 95 ms
18,648 KB
testcase_35 AC 98 ms
18,592 KB
testcase_36 AC 97 ms
18,488 KB
testcase_37 AC 127 ms
18,504 KB
testcase_38 AC 96 ms
18,592 KB
testcase_39 AC 99 ms
18,596 KB
testcase_40 AC 100 ms
18,624 KB
testcase_41 AC 97 ms
18,624 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

using A = array<mint,5>;

A op(A a,A b){
	rep(i,5){
		a[i] += b[i];
	}
	return a;
}

A e(){
	A ret;
	rep(i,5)ret[i] = 0;
	return ret;
}

A mapping(int a,A b){
	if(a==-1)return b;
	mint t = 1;
	rep(i,5){
		b[i] = b[0] * t;
		t *= a;
	}
	return b;
}

int composition(int a,int b){
	if(a==-1)return b;
	if(b==-1)return a;
	return a;
}

int id(){
	return -1;
}

int main(){
	
	int N;
	cin>>N;
	
	vector<A> AA(N);
	rep(i,N){
		int t;
		scanf("%d",&t);
		mint cur = 1;
		rep(j,5){
			AA[i][j] = cur;
			cur *= t;
		}
	}
	
	lazy_segtree<A,op,e,int,mapping,composition,id> seg(AA);
	
	int Q;
	cin>>Q;
	
	rep(_,Q){
		
		int t;
		scanf("%d",&t);
		
		if(t==0){
			int u,v,w,b;
			scanf("%d %d %d %d",&u,&v,&w,&b);
			u--;v--;w--;
			if(u>v)swap(u,v);
			vector<int> l,r;
			if(u<w&&w<v)l.push_back(u),r.push_back(v+1);
			else{
				l.push_back(0),r.push_back(u+1);
				l.push_back(v),r.push_back(N);
			}
			rep(j,l.size()){
				seg.apply(l[j],r[j],b);
			}
		}
		else{
			mint ans = 0;
			int u,v,w;
			scanf("%d %d %d",&u,&v,&w);
			u--;v--;w--;
			if(u>v)swap(u,v);
			vector<int> l,r;
			if(u<w&&w<v)l.push_back(u),r.push_back(v+1);
			else{
				l.push_back(0),r.push_back(u+1);
				l.push_back(v),r.push_back(N);
			}
			A ret = e();
			rep(j,l.size()){
				ret = op(ret,seg.prod(l[j],r[j]));
			}
			mint M = ret[1];
			M /= ret[0];
			if(t==1){
				ans += ret[1];
				ans -= M * ret[0];
			}
			if(t==2){
				ans += ret[2];
				ans -= M * ret[1] * 2;
				ans += M * M * ret[0];
			}
			if(t==3){
				ans += ret[3];
				ans -= M * ret[2] * 3;
				ans += M * M * ret[1] * 3;
				ans -= M * M * M * ret[0];
			}
			if(t==4){
				ans += ret[4];
				ans -= ret[3] * M * 4;
				ans += M * M * ret[2] * 6;
				ans -= M * M * M * ret[1] * 4;
				ans += M * M *M * M * ret[0];
			}
			ans /= ret[0];
			printf("%d\n",ans.val());
			
		}
		
	}
	
	
	return 0;
}
0