Submission #1837480

---

Source Code Expand

#include<bits/stdc++.h>
#define range(i,a,b) for(int i = (a); i < (b); i++)
#define rep(i,b) for(int i = 0; i < (b); i++)
#define all(a) (a).begin(), (a).end()
#define show(x)  cerr << #x << " = " << (x) << endl;
//const int INF = 1e8;
using namespace std;

const int M = 1000000007;

//x^n mod M
typedef long long ull;
ull power(ull x, ull n, ull M){
	ull res = 1;
	if(n > 0){
		res = power(x, n / 2, M);
		if(n % 2 == 0) res = (res * res) % M;
		else res = (((res * res) % M) * x ) %M;
	}
	return res;
}

int main(){
	int n;
	cin >> n;

	vector<pair<int, bool>> p(n * 2); //1 : PC, 0 : power supply
	rep(i,n){
		int a;
		cin >> a;
		p[i] = make_pair(a, true);
	}
	rep(i,n){
		int b;
		cin >> b;
		p[i + n] = make_pair(b, false);
	}


	sort(all(p));

	int cnt = 0, cur = -1;
	vector<int> v;
	rep(i,n * 2){
		cnt += p[i].second ? 1 : -1;
		//show(cnt)
		if(cnt == 0){
			v.emplace_back((i - cur) / 2);
			cur = i;
		}
	}

	//for(auto i : p){ cout << i.second << ' ' ; } cout << endl;

	long long ans = 1;
	long long fact[100005];
	fact[0] = 1;
	range(i,1,n + 1){
		fact[i] = fact[i - 1] * i;
	}

	for(auto i : v){
		(ans *= fact[i]) %= M;
	}
	cout << ans << endl;
}

Submission Info

Judge Result