結果

問題 No.1548 [Cherry 2nd Tune B] 貴方と私とサイクルとモーメント
ユーザー 沙耶花沙耶花
提出日時 2021-06-11 23:11:12
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 126 ms / 4,500 ms
コード長 2,069 bytes
コンパイル時間 5,030 ms
コンパイル使用メモリ 268,168 KB
実行使用メモリ 18,608 KB
最終ジャッジ日時 2023-08-21 14:04:38
合計ジャッジ時間 10,686 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 96 ms
10,672 KB
testcase_03 AC 28 ms
6,976 KB
testcase_04 AC 86 ms
17,684 KB
testcase_05 AC 10 ms
4,376 KB
testcase_06 AC 69 ms
10,196 KB
testcase_07 AC 75 ms
18,144 KB
testcase_08 AC 98 ms
10,824 KB
testcase_09 AC 69 ms
7,396 KB
testcase_10 AC 61 ms
17,332 KB
testcase_11 AC 71 ms
17,268 KB
testcase_12 AC 48 ms
5,120 KB
testcase_13 AC 42 ms
17,364 KB
testcase_14 AC 72 ms
6,792 KB
testcase_15 AC 45 ms
18,300 KB
testcase_16 AC 77 ms
11,252 KB
testcase_17 AC 96 ms
17,104 KB
testcase_18 AC 79 ms
18,156 KB
testcase_19 AC 37 ms
6,868 KB
testcase_20 AC 54 ms
17,676 KB
testcase_21 AC 61 ms
10,148 KB
testcase_22 AC 123 ms
18,368 KB
testcase_23 AC 123 ms
18,468 KB
testcase_24 AC 123 ms
18,364 KB
testcase_25 AC 126 ms
18,392 KB
testcase_26 AC 124 ms
18,392 KB
testcase_27 AC 125 ms
18,416 KB
testcase_28 AC 125 ms
18,608 KB
testcase_29 AC 125 ms
18,420 KB
testcase_30 AC 124 ms
18,420 KB
testcase_31 AC 123 ms
18,344 KB
testcase_32 AC 92 ms
18,472 KB
testcase_33 AC 95 ms
18,516 KB
testcase_34 AC 95 ms
18,468 KB
testcase_35 AC 95 ms
18,476 KB
testcase_36 AC 95 ms
18,476 KB
testcase_37 AC 125 ms
18,416 KB
testcase_38 AC 93 ms
18,420 KB
testcase_39 AC 95 ms
18,492 KB
testcase_40 AC 94 ms
18,424 KB
testcase_41 AC 94 ms
18,392 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