結果

問題 No.2717 Sum of Subarray of Subsequence
ユーザー t9unkubjt9unkubj
提出日時 2025-01-15 18:05:55
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 121 ms / 2,000 ms
コード長 5,644 bytes
コンパイル時間 5,451 ms
コンパイル使用メモリ 312,796 KB
実行使用メモリ 54,016 KB
最終ジャッジ日時 2025-01-15 18:06:05
合計ジャッジ時間 9,324 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#ifdef t9unkubj
#include"debug.cpp"
//#include"template_no_debug.h"
#else
#define dbg(...) 199958
#endif
#undef _GLIBCXX_DEBUG
#pragma GCC optimize("O3")
using namespace std;
#include<bits/stdc++.h>
using ll=long long;
using ull=unsigned long long;
template<class T>using vc=vector<T>;
template<class T>using vvc=vc<vc<T>>;
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
#define REP(i,j,n) for(ll i=(j);i<(ll)(n);i++)
#define DREP(i,n,m) for(ll i=(n);i>=(m);i--)
#define drep(i,n) for(ll i=((n)-1);i>=0;i--)
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
template<class T,class F>
bool chmin(T &x, F y){
if(x>y){
x=y;
return true;
}
return false;
}
template<class T, class F>
bool chmax(T &x, F y){
if(x<y){
x=y;
return true;
}
return false;
}
double pass_time=0;
template<class T,class G,class F>
F modpow(T x,G p,F m){
T ret=1%m;
x%=m;
while(p){
if(p&1)ret=(1ll*ret*x)%m;
x=(1ll*x*x)%m;
p>>=1;
}
return ret;
}
template<class T>
T extgcd(T a,T b,T&x,T&y){//ax+by=gcd(a,b)
if(b==0){
x=1;
y=0;
return a;
}else{
T res=extgcd(b,a%b,y,x);
y-=(a/b)*x;
return res;
}
}
template<class T>
pair<T,T> inv(T x,T m){
T a1,a2;
T res=extgcd(x,m,a1,a2);
return {a1,m/res};
}
constexpr int mod=998244353;
struct mint{
long long val;
inline long long fast(long long x){
if(x<mod&&x>=0)return x;
x%=mod;
if(x<0)x+=mod;
return x;
}
mint():val(0){}
mint(long long val):val(fast(val)){}
mint power(long long m)const {
mint res(1);
mint ret(*this);
while(m){
if(m&1)res*=ret;
ret*=ret;
m>>=1;
}
return res;
}
mint& operator++() {
val++;
if (val == mod) val = 0;
return *this;
}
mint& operator--() {
if (val == 0)val=mod;
val--;
return *this;
}
mint operator++(int) {
mint result = *this;
++*this;
return result;
}
mint operator--(int) {
mint result = *this;
--*this;
return result;
}
mint operator-() const {
return mint(-val);
}
friend mint operator +(const mint&a,const mint&b) noexcept{
return mint(a)+=b;
}
friend mint operator -(const mint&a,const mint&b) noexcept{
return mint(a)-=b;
}
friend mint operator *(const mint&a,const mint&b) noexcept{
return mint(a)*=b;
}
friend mint operator /(const mint&a,const mint&b) noexcept{
return mint(a)/=b;
}
mint& operator+=(const mint&a)noexcept{
val+=a.val;
if(val>=mod)val-=mod;
return *this;
}
mint& operator-=(const mint&a)noexcept{
val-=a.val;
if(val<0)val+=mod;
return *this;
}
mint& operator*=(const mint&a){
val*=a.val;
val=fast(val);
return *this;
}
mint& operator/=(const mint&a){
val*=inv<long long>(a.val,mod).first;
val=fast(val);
return *this;
}
bool operator == (const mint&x)const noexcept{
return this->val==x.val;
}
bool operator != (const mint&x)const noexcept{
return this->val!=x.val;
}
friend ostream& operator << (ostream &os, const mint &x) noexcept {
return os << x.val;
}
friend istream& operator >> (istream &is, mint &x) noexcept {
long long v;
is >> v;
x=mint(v);
return is;
}
};
vector<mint>fact(1,1),invfact(1,1);
void build(int n){
if(n<(int)fact.size())return;
fact=invfact=vector<mint>(n+1);
fact[0]=1;
for(int i=1;i<=n;i++)fact[i]=fact[i-1]*i;
invfact[n]=(1/fact[n]);
for(int i=n-1;i>=0;i--)invfact[i]=invfact[i+1]*(i+1);
}
mint C(int a,int b){//aCb
if(a<0||b<0||a-b<0)return mint(0);
while((int)fact.size()<=a){
fact.push_back(fact.back()*(fact.size()));
}
while((int)invfact.size()<=a){
invfact.push_back(invfact.back()/invfact.size());
}
return fact[a]*invfact[b]*invfact[a-b];
}
mint P(int a,int b){
if(a<b||b<0)return 0;
return fact[a]*invfact[a-b];
}
//ab
mint H(int a,int b){
return C(a+b-1,b);
}
void print(mint a) { cout << a; }
template<class T>
pair<T,T> mod_solve(T a,T b,T m){//ax=b mod mx
a%=m,b%=m;if(a<0)a+=m;if(b<0)b+=m;
T g=gcd(gcd(a,b),m);
a/=g,b/=g,m/=g;
if(gcd(a,m)>1)return {-1,-1};
return {(inv(a,m).first*b)%m,inv(a,m).second};
}
//
//
//
void solve(){
int n;
cin>>n;
vc dp(n+1,vc(3,vc<mint>(2)));
dp[0][0][0]=1;
rep(i,n){
rep(j,3)dp[i+1][j][0]+=dp[i][j][0]+dp[i][j][1];
rep(j,2)dp[i+1][0][1]+=dp[i][0][j];
rep(j,2)dp[i+1][1][1]+=dp[i][0][j];
rep(j,2)dp[i+1][1][1]+=dp[i][1][j];
rep(j,2)dp[i+1][2][1]+=dp[i][1][j];
rep(j,2)dp[i+1][2][1]+=dp[i][2][j];
}
vvc<mint>dp2(n+1,vc<mint>(3));
dp2[0][1]=1;
rep(i,n){
dp2[i+1]=dp2[i];
dp2[i+1][1]+=dp2[i][1];
dp2[i+1][2]+=dp2[i][1];
dp2[i+1][2]+=dp2[i][2];
}
vc<int>a(n);
rep(i,n)cin>>a[i];
mint ans=0;
rep(i,n)
ans+=a[i]*(dp[i+1][1][1])*(dp2[n-i-1][1]+dp2[n-i-1][2]);
cout<<ans<<"\n";
}
signed main(){
cin.tie(0)->sync_with_stdio(0);
pass_time=clock();
int t=1;
//cin>>t;
while(t--)solve();
pass_time=clock()-pass_time;
dbg(pass_time/CLOCKS_PER_SEC);
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0