結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 103 ms
10,940 KB
testcase_03 AC 30 ms
7,296 KB
testcase_04 AC 88 ms
17,920 KB
testcase_05 AC 10 ms
5,376 KB
testcase_06 AC 71 ms
10,496 KB
testcase_07 AC 77 ms
18,520 KB
testcase_08 AC 101 ms
10,936 KB
testcase_09 AC 71 ms
7,424 KB
testcase_10 AC 60 ms
17,624 KB
testcase_11 AC 77 ms
17,712 KB
testcase_12 AC 50 ms
5,504 KB
testcase_13 AC 45 ms
17,868 KB
testcase_14 AC 75 ms
7,040 KB
testcase_15 AC 46 ms
18,648 KB
testcase_16 AC 78 ms
11,392 KB
testcase_17 AC 97 ms
17,356 KB
testcase_18 AC 82 ms
18,300 KB
testcase_19 AC 38 ms
7,040 KB
testcase_20 AC 56 ms
17,888 KB
testcase_21 AC 68 ms
10,640 KB
testcase_22 AC 127 ms
18,592 KB
testcase_23 AC 126 ms
18,524 KB
testcase_24 AC 124 ms
18,704 KB
testcase_25 AC 125 ms
18,688 KB
testcase_26 AC 128 ms
18,688 KB
testcase_27 AC 125 ms
18,688 KB
testcase_28 AC 124 ms
18,652 KB
testcase_29 AC 127 ms
18,660 KB
testcase_30 AC 125 ms
18,612 KB
testcase_31 AC 123 ms
18,672 KB
testcase_32 AC 100 ms
18,728 KB
testcase_33 AC 98 ms
18,676 KB
testcase_34 AC 96 ms
18,688 KB
testcase_35 AC 96 ms
18,688 KB
testcase_36 AC 96 ms
18,688 KB
testcase_37 AC 127 ms
18,568 KB
testcase_38 AC 95 ms
18,688 KB
testcase_39 AC 97 ms
18,660 KB
testcase_40 AC 96 ms
18,692 KB
testcase_41 AC 102 ms
18,488 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