結果

問題 No.1193 Penguin Sequence
ユーザー HIR180HIR180
提出日時 2020-08-22 14:37:29
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 181 ms / 2,000 ms
コード長 2,676 bytes
コンパイル時間 3,295 ms
コンパイル使用メモリ 244,876 KB
実行使用メモリ 12,292 KB
最終ジャッジ日時 2024-04-23 08:57:47
合計ジャッジ時間 9,756 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 128 ms
11,256 KB
testcase_01 AC 164 ms
10,960 KB
testcase_02 AC 169 ms
11,240 KB
testcase_03 AC 181 ms
12,272 KB
testcase_04 AC 162 ms
12,268 KB
testcase_05 AC 163 ms
11,380 KB
testcase_06 AC 161 ms
11,896 KB
testcase_07 AC 160 ms
10,832 KB
testcase_08 AC 163 ms
11,772 KB
testcase_09 AC 174 ms
11,636 KB
testcase_10 AC 169 ms
11,640 KB
testcase_11 AC 110 ms
10,316 KB
testcase_12 AC 110 ms
11,232 KB
testcase_13 AC 145 ms
12,292 KB
testcase_14 AC 136 ms
11,204 KB
testcase_15 AC 157 ms
11,536 KB
testcase_16 AC 122 ms
12,152 KB
testcase_17 AC 29 ms
9,964 KB
testcase_18 AC 42 ms
10,516 KB
testcase_19 AC 171 ms
10,956 KB
testcase_20 AC 134 ms
10,444 KB
testcase_21 AC 116 ms
10,816 KB
testcase_22 AC 44 ms
10,380 KB
testcase_23 AC 109 ms
11,236 KB
testcase_24 AC 98 ms
10,768 KB
testcase_25 AC 60 ms
10,920 KB
testcase_26 AC 38 ms
9,692 KB
testcase_27 AC 143 ms
11,656 KB
testcase_28 AC 108 ms
10,568 KB
testcase_29 AC 144 ms
11,156 KB
testcase_30 AC 66 ms
10,952 KB
testcase_31 AC 60 ms
10,164 KB
testcase_32 AC 121 ms
11,392 KB
testcase_33 AC 94 ms
10,320 KB
testcase_34 AC 81 ms
11,372 KB
testcase_35 AC 91 ms
12,040 KB
testcase_36 AC 77 ms
11,436 KB
testcase_37 AC 137 ms
10,572 KB
testcase_38 AC 31 ms
11,368 KB
testcase_39 AC 31 ms
10,604 KB
testcase_40 AC 31 ms
10,220 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

//Let's join Kaede Takagaki Fan Club !!
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#define all(x) x.begin(),x.end()
template<class T>
void dmp(T a){
	rep(i,a.size()) cout << a[i] << " ";
	cout << endl;
}
template<class T>
bool chmax(T&a, T b){
	if(a < b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
bool chmin(T&a, T b){
	if(a > b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
void g(T &a){
	cin >> a;
}
template<class T>
void o(const T &a,bool space=false){
	cout << a << (space?' ':'\n');
}
//ios::sync_with_stdio(false);
const ll mod = 998244353;
template<class T>
void add(T&a,T b){
	a+=b;
	if(a >= mod) a-=mod;
}
ll modpow(ll x,ll n){
	ll res=1;
	while(n>0){
		if(n&1) res=res*x%mod;
		x=x*x%mod;
		n>>=1;
	}
	return res;
}
ll F[200005],R[200005];
void make(){
	F[0] = 1;
	for(int i=1;i<200005;i++) F[i] = F[i-1]*i%mod;
	for(int i=0;i<200005;i++) R[i] = modpow(F[i],mod-2);
}
ll C(int a,int b){
	return F[a]*R[b]%mod*R[a-b]%mod;
}

int n, a[200005];
vector<int>vi;
struct BIT
{
	int bit[(1<<20)+5];
	int f(int x)
	{
		return x&-x;
	}
	void add(int i,int x)
	{
		for(int s=i;s<=(1<<20);s+=f(s)) bit[s]+=x;
	}
	int sum(int i)
	{
		int res=0;
		for(int s=i;s>0;s-=f(s)) res+=bit[s];
		return res;
	}
}bit;
int main(){
	scanf("%d", &n);
	repn(i, n) scanf("%d", &a[i]);
	repn(i, n) vi.pb(a[i]);
	SORT(vi); ERASE(vi);
	ll R = 0, S = 0;
	repn(i, n){
		a[i] = POSL(vi, a[i]) + 1;
		R += i-1-bit.sum(a[i]);
		S += i-1-bit.sum(a[i])+bit.sum(a[i]-1);
		bit.add(a[i], 1);
	}
	make();
	ll ans = 1;
	repn(i, n) ans = ans * C(n, i) % mod;
	ll c1 = 0;
	c1 = 1LL*n*(n+1)/2%mod;
	c1 = c1*c1%mod;
	repn(i, n) c1 -= 1LL*i*i%mod;
	c1 = (c1%mod+mod)%mod;
	c1 = c1 * ((mod+1)/2) % mod;
	c1 = c1 * modpow(n, mod-2) % mod * modpow(n, mod-2) % mod;
	ll c2 = 0;
	for(int i=1;i<=n;i++){
		c2 += 1LL*i*(i-1)%mod * modpow(n, mod-2) % mod * modpow(n-1, mod-2) % mod;
	}
	c2 %= mod;
	c1 = c1 * ans % mod;
	c2 = c2 * ans % mod;
	ll ret = S%mod*c1%mod;
	ret += R%mod*c2%mod;
	cout << (ret%mod) << endl;
}
0