結果

問題 No.1193 Penguin Sequence
ユーザー HIR180HIR180
提出日時 2020-08-22 14:37:29
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 178 ms / 2,000 ms
コード長 2,676 bytes
コンパイル時間 3,783 ms
コンパイル使用メモリ 247,544 KB
実行使用メモリ 9,052 KB
最終ジャッジ日時 2024-10-15 08:52:39
合計ジャッジ時間 10,407 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 139 ms
8,796 KB
testcase_01 AC 178 ms
8,924 KB
testcase_02 AC 170 ms
9,052 KB
testcase_03 AC 174 ms
8,928 KB
testcase_04 AC 169 ms
8,924 KB
testcase_05 AC 173 ms
8,928 KB
testcase_06 AC 170 ms
9,028 KB
testcase_07 AC 167 ms
9,052 KB
testcase_08 AC 170 ms
8,920 KB
testcase_09 AC 175 ms
8,924 KB
testcase_10 AC 174 ms
8,928 KB
testcase_11 AC 110 ms
7,952 KB
testcase_12 AC 111 ms
7,960 KB
testcase_13 AC 148 ms
8,692 KB
testcase_14 AC 139 ms
8,384 KB
testcase_15 AC 159 ms
8,752 KB
testcase_16 AC 130 ms
8,160 KB
testcase_17 AC 30 ms
6,784 KB
testcase_18 AC 43 ms
6,912 KB
testcase_19 AC 168 ms
9,048 KB
testcase_20 AC 135 ms
8,496 KB
testcase_21 AC 118 ms
8,124 KB
testcase_22 AC 44 ms
7,040 KB
testcase_23 AC 101 ms
7,784 KB
testcase_24 AC 94 ms
7,544 KB
testcase_25 AC 62 ms
7,272 KB
testcase_26 AC 39 ms
6,784 KB
testcase_27 AC 147 ms
8,552 KB
testcase_28 AC 111 ms
7,968 KB
testcase_29 AC 150 ms
8,688 KB
testcase_30 AC 73 ms
7,400 KB
testcase_31 AC 66 ms
7,424 KB
testcase_32 AC 126 ms
8,160 KB
testcase_33 AC 98 ms
7,748 KB
testcase_34 AC 85 ms
7,552 KB
testcase_35 AC 97 ms
7,884 KB
testcase_36 AC 85 ms
7,544 KB
testcase_37 AC 145 ms
8,536 KB
testcase_38 AC 33 ms
6,784 KB
testcase_39 AC 32 ms
6,784 KB
testcase_40 AC 31 ms
6,656 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