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

• 問題: No.1838 Modulo Straight
• ユーザー: 👑 potato167
• 提出日時: 2022-02-11 22:04:25
• 言語: C++17 (gcc 13.2.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 1,539 ms / 2,000 ms
• コード長: 3,400 bytes
• コンパイル時間: 4,114 ms
• コンパイル使用メモリ: 215,700 KB
• 実行使用メモリ: 308,916 KB
• 最終ジャッジ日時: 2023-09-10 02:42:36
• 合計ジャッジ時間: 25,651 ms
• ジャッジサーバーID （参考情報）: judge12 / judge13

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 4,380 KB
testcase_01	AC	2 ms; 4,380 KB
testcase_02	AC	2 ms; 4,376 KB
testcase_03	AC	2 ms; 4,376 KB
testcase_04	AC	2 ms; 4,376 KB
testcase_05	AC	3 ms; 4,376 KB
testcase_06	AC	2 ms; 4,376 KB
testcase_07	AC	3 ms; 4,376 KB
testcase_08	AC	3 ms; 4,376 KB
testcase_09	AC	4 ms; 4,376 KB
testcase_10	AC	6 ms; 5,684 KB
testcase_11	AC	1,425 ms; 308,888 KB
testcase_12	AC	1,435 ms; 308,768 KB
testcase_13	AC	1,343 ms; 308,868 KB
testcase_14	AC	1,355 ms; 308,752 KB
testcase_15	AC	1,539 ms; 308,916 KB
testcase_16	AC	968 ms; 161,888 KB
testcase_17	AC	1,087 ms; 161,824 KB
testcase_18	AC	1,032 ms; 161,872 KB
testcase_19	AC	865 ms; 113,052 KB
testcase_20	AC	850 ms; 113,052 KB
testcase_21	AC	842 ms; 113,000 KB
testcase_22	AC	131 ms; 18,280 KB
testcase_23	AC	148 ms; 17,752 KB
testcase_24	AC	160 ms; 17,032 KB
testcase_25	AC	169 ms; 16,848 KB
testcase_26	AC	182 ms; 16,576 KB
testcase_27	AC	274 ms; 17,520 KB
testcase_28	AC	273 ms; 17,772 KB
testcase_29	AC	299 ms; 18,696 KB
testcase_30	AC	346 ms; 22,500 KB
testcase_31	AC	391 ms; 29,644 KB
testcase_32	AC	564 ms; 61,024 KB
testcase_33	AC	967 ms; 161,784 KB
testcase_34	AC	956 ms; 161,784 KB
testcase_35	AC	693 ms; 88,388 KB
testcase_36	AC	487 ms; 44,312 KB
testcase_37	AC	358 ms; 25,104 KB
testcase_38	AC	288 ms; 17,696 KB
testcase_39	AC	207 ms; 97,440 KB
testcase_40	AC	2 ms; 4,380 KB

ソースコード

#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";
}