結果

問題 No.1193 Penguin Sequence
ユーザー HIR180HIR180
提出日時 2020-08-22 14:37:29
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 174 ms / 2,000 ms
コード長 2,676 bytes
コンパイル時間 3,441 ms
コンパイル使用メモリ 245,708 KB
実行使用メモリ 10,888 KB
最終ジャッジ日時 2023-08-05 11:47:55
合計ジャッジ時間 9,539 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 143 ms
10,888 KB
testcase_01 AC 174 ms
10,592 KB
testcase_02 AC 170 ms
10,740 KB
testcase_03 AC 171 ms
10,744 KB
testcase_04 AC 167 ms
10,704 KB
testcase_05 AC 168 ms
10,628 KB
testcase_06 AC 168 ms
10,652 KB
testcase_07 AC 166 ms
10,700 KB
testcase_08 AC 167 ms
10,652 KB
testcase_09 AC 169 ms
10,868 KB
testcase_10 AC 170 ms
10,604 KB
testcase_11 AC 106 ms
10,068 KB
testcase_12 AC 108 ms
10,048 KB
testcase_13 AC 150 ms
10,440 KB
testcase_14 AC 139 ms
10,564 KB
testcase_15 AC 160 ms
10,780 KB
testcase_16 AC 132 ms
10,120 KB
testcase_17 AC 32 ms
9,516 KB
testcase_18 AC 43 ms
9,676 KB
testcase_19 AC 169 ms
10,772 KB
testcase_20 AC 140 ms
10,260 KB
testcase_21 AC 118 ms
10,112 KB
testcase_22 AC 45 ms
9,696 KB
testcase_23 AC 102 ms
10,036 KB
testcase_24 AC 91 ms
9,864 KB
testcase_25 AC 58 ms
9,996 KB
testcase_26 AC 37 ms
9,760 KB
testcase_27 AC 145 ms
10,536 KB
testcase_28 AC 111 ms
10,024 KB
testcase_29 AC 144 ms
10,464 KB
testcase_30 AC 71 ms
9,796 KB
testcase_31 AC 62 ms
9,828 KB
testcase_32 AC 118 ms
10,304 KB
testcase_33 AC 93 ms
10,036 KB
testcase_34 AC 86 ms
9,920 KB
testcase_35 AC 94 ms
10,040 KB
testcase_36 AC 79 ms
9,960 KB
testcase_37 AC 143 ms
10,392 KB
testcase_38 AC 32 ms
9,656 KB
testcase_39 AC 32 ms
9,440 KB
testcase_40 AC 30 ms
9,492 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