結果

問題 No.1604 Swap Sort:ONE
ユーザー 👑 potato167
提出日時 2021-07-16 21:21:51
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 3,178 bytes
コンパイル時間 1,962 ms
コンパイル使用メモリ 199,848 KB
最終ジャッジ日時 2025-01-23 01:17:28
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
#define _GLIBCXX_DEBUG
using namespace std;
using std::cout;
using std::cin;
using std::endl;
using ll=long long;
using ld=long double;
ll I=1167167167167167167;
ll mod=1e9+7;
#define rep(i,a) for (ll i=0;i<a;i++)
template<class T> using _pq = priority_queue<T, vector<T>, greater<T>>;
template<class T> ll LB(vector<T> &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();}
template<class T> ll UB(vector<T> &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();}
template<class T> bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;}
template<class T> bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}else return 0;}
template<class T> void So(vector<T> &v) {sort(v.begin(),v.end());}
template<class T> void Sore(vector<T> &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});}
template<class T> void print_tate(vector<T> &v) {rep(i,v.size()) cout<<v[i]<<"\n";}
void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";}


namespace po167{
template <class T>
struct seg_sum{
	int n;
	int seg_size;
	std::vector<T> seg;
	T e;
	//要素数、単位元 typeは+演算が定義されているもの
	seg_sum(int k,T e_in){
		e=e_in;
		n=k;
		seg_size=1;
		while(seg_size<n) seg_size*=2;
		seg.resize(seg_size*2-1);
		for(int i=0;i<seg_size*2-1;i++) seg[i]=e;
	}
 
	T sum_sub(int a,int b,int l,int r,int k){
		if(r<=a||b<=l) return e;
		if(a<=l&&r<=b) return seg[k];
		T X=sum_sub(a,b,l,(l+r)/2,k*2+1);
		T Y=sum_sub(a,b,(l+r)/2,r,k*2+2);
		return X+Y;
	}
	T sum(int a,int b){
		//assert(0<=a&&b<n);
		return sum_sub(a,b,0,seg_size,0);
	}
	T get(int index){
		assert(0<=index&&index<n);
		return seg[seg_size-1+index];
	}
	void insert(int index,T a){
		assert(0<=index&&index<n);
		int now_index=index+seg_size-1;
		seg[now_index]=a;
		while(now_index!=0){
			now_index-=1;
			now_index/=2;
			seg[now_index]=seg[now_index*2+1]+seg[now_index*2+2];
		}
	}
};
 
 
template <class T>
//p,qのなかみは、重複ありで内容が同じならおけ
long long inversion(std::vector<T> &p,std::vector<T> &q){
	int n=p.size();
	assert(p.size()==q.size());
	po167::seg_sum<int> S(n,0);
	std::vector<int> p_index(n,-1);
	std::map<T,std::queue<int>> m;
	for(int i=0;i<n;i++) m[p[i]].push(i);
	for(int i=0;i<n;i++){
		if(m[q[i]].empty()) return -1;
		int a=m[q[i]].front();
		m[q[i]].pop();
		p_index[i]=a;
	}
	long long ans=0;
	for(int i=0;i<n;i++){
		ans+=S.sum(p_index[i],n);
		S.insert(p_index[i],1);
	}
	return ans;
}
 
template <class T>
//qのなかみは、0~nの順列の並び替えである必要がある
long long inversion(std::vector<T> &q){
	int n=q.size();
	po167::seg_sum<int> S(n,0);
	std::vector<int> p_index(n);
	for(int i=0;i<n;i++){
		p_index[i]=i;
	}
	long long ans=0;
	for(int i=0;i<n;i++){
		assert(0<=q[i]&&q[i]<n);
		assert(S.get(q[i])==0);
		ans+=S.sum(q[i],n);
		S.insert(q[i],1);
	}
	return ans;
}
}
using po167::inversion;
using po167::seg_sum;


//おちこんだりもしたけれど、私はげんきです。
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	
	int N;
	cin>>N;
	vector<int> p(N);
	rep(i,N){
		cin>>p[i];
		p[i]--;
	}
	cout<<inversion(p)<<endl;
}
0