#738481 (C++17) No.1838 Modulo Straight

提出ソース

結果

問題	No.1838 Modulo Straight
ユーザー	👑 potato167
提出日時	2022-02-11 22:04:25
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	1,451 ms / 2,000 ms
コード長	3,400 bytes
コンパイル時間	2,419 ms
コンパイル使用メモリ	219,060 KB
実行使用メモリ	308,820 KB
最終ジャッジ日時	2024-06-27 18:51:23
合計ジャッジ時間	22,386 ms
ジャッジサーバーID （参考情報）	judge2 / judge1

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	1 ms 5,376 KB
testcase_02	AC	2 ms 5,376 KB
testcase_03	AC	3 ms 5,376 KB
testcase_04	AC	3 ms 5,376 KB
testcase_05	AC	2 ms 5,376 KB
testcase_06	AC	3 ms 5,376 KB
testcase_07	AC	4 ms 5,376 KB
testcase_08	AC	4 ms 5,376 KB
testcase_09	AC	4 ms 5,376 KB
testcase_10	AC	6 ms 5,632 KB
testcase_11	AC	1,349 ms 308,772 KB
testcase_12	AC	1,371 ms 308,776 KB
testcase_13	AC	1,197 ms 308,736 KB
testcase_14	AC	1,262 ms 308,736 KB
testcase_15	AC	1,451 ms 308,820 KB
testcase_16	AC	881 ms 161,856 KB
testcase_17	AC	954 ms 161,920 KB
testcase_18	AC	946 ms 161,920 KB
testcase_19	AC	743 ms 112,916 KB
testcase_20	AC	746 ms 112,896 KB
testcase_21	AC	747 ms 112,896 KB
testcase_22	AC	128 ms 18,200 KB
testcase_23	AC	143 ms 17,792 KB
testcase_24	AC	152 ms 16,972 KB
testcase_25	AC	162 ms 16,768 KB
testcase_26	AC	175 ms 16,768 KB
testcase_27	AC	266 ms 17,536 KB
testcase_28	AC	267 ms 17,472 KB
testcase_29	AC	289 ms 18,672 KB
testcase_30	AC	321 ms 22,400 KB
testcase_31	AC	371 ms 29,748 KB
testcase_32	AC	501 ms 60,728 KB
testcase_33	AC	857 ms 161,968 KB
testcase_34	AC	864 ms 162,048 KB
testcase_35	AC	614 ms 88,448 KB
testcase_36	AC	441 ms 44,448 KB
testcase_37	AC	333 ms 24,700 KB
testcase_38	AC	280 ms 17,664 KB
testcase_39	AC	199 ms 97,280 KB
testcase_40	AC	2 ms 5,376 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#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 ILL=1167167167167167167;
const int INF=2100000000;
const ll mod=998244353;
#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;});}
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;




void solve();

//  rainy ~ 雨に打たれて ~
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	
	int t=1;
	//cin>>t;
	rep(i,t) solve();
}

void solve(){
	ll M,K;
	cin>>M>>K;
	vector<ll> q(K),p(M*K),base(M*K);
	vector<ll> diff(M),ind(M);
	rep(i,M*K){
		int a;
		cin>>a;
		diff[a]+=M-1-2*q[ind[a]];
		q[ind[a]]++;
		ind[a]++;
		p[i]=a;
		base[i]=i%M;
	}
	ll tmp=inversion(p,base);
	ll ans=tmp;
	rep(i,M){
		tmp+=diff[i];
		chmin(ans,tmp);
		//cout<<tmp<<"\n";
	}
	cout<<ans<<"\n";
}

yukicoder

結果

テストケース

ソースコード