QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#106415#6308. MagicPure_FuriesRE 1ms5620kbC++144.7kb2023-05-17 18:09:372023-05-17 18:09:42

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-05-17 18:09:42]
  • 评测
  • 测评结果:RE
  • 用时:1ms
  • 内存:5620kb
  • [2023-05-17 18:09:37]
  • 提交

answer

#include<bits/stdc++.h>

namespace atcoder {

namespace internal {

template <class T> struct simple_queue {
	std::vector<T> payload;
	int pos = 0;
	void reserve(int n) { payload.reserve(n); }
	int size() const { return int(payload.size()) - pos; }
	bool empty() const { return pos == int(payload.size()); }
	void push(const T& t) { payload.push_back(t); }
	T& front() { return payload[pos]; }
	void clear() {
		payload.clear();
		pos = 0;
	}
	void pop() { pos++; }
};

}  // namespace internal
template <class Cap> struct mf_graph {
  public:
	mf_graph() : _n(0) {}
	mf_graph(int n) : _n(n), g(n) {}

	int add_edge(int from, int to, Cap cap) {
		assert(0 <= from && from < _n);
		assert(0 <= to && to < _n);
		assert(0 <= cap);
		int m = int(pos.size());
		pos.push_back({from, int(g[from].size())});
		int from_id = int(g[from].size());
		int to_id = int(g[to].size());
		if (from == to) to_id++;
		g[from].push_back(_edge{to, to_id, cap});
		g[to].push_back(_edge{from, from_id, 0});
		return m;
	}

	struct edge {
		int from, to;
		Cap cap, flow;
	};

	edge get_edge(int i) {
		int m = int(pos.size());
		assert(0 <= i && i < m);
		auto _e = g[pos[i].first][pos[i].second];
		auto _re = g[_e.to][_e.rev];
		return edge{pos[i].first, _e.to, _e.cap + _re.cap, _re.cap};
	}
	std::vector<edge> edges() {
		int m = int(pos.size());
		std::vector<edge> result;
		for (int i = 0; i < m; i++) {
			result.push_back(get_edge(i));
		}
		return result;
	}
	void change_edge(int i, Cap new_cap, Cap new_flow) {
		int m = int(pos.size());
		assert(0 <= i && i < m);
		assert(0 <= new_flow && new_flow <= new_cap);
		auto& _e = g[pos[i].first][pos[i].second];
		auto& _re = g[_e.to][_e.rev];
		_e.cap = new_cap - new_flow;
		_re.cap = new_flow;
	}

	Cap flow(int s, int t) {
		return flow(s, t, std::numeric_limits<Cap>::max());
	}
	Cap flow(int s, int t, Cap flow_limit) {
		assert(0 <= s && s < _n);
		assert(0 <= t && t < _n);
		assert(s != t);

		std::vector<int> level(_n), iter(_n);
		internal::simple_queue<int> que;

		auto bfs = [&]() {
			std::fill(level.begin(), level.end(), -1);
			level[s] = 0;
			que.clear();
			que.push(s);
			while (!que.empty()) {
				int v = que.front();
				que.pop();
				for (auto e : g[v]) {
					if (e.cap == 0 || level[e.to] >= 0) continue;
					level[e.to] = level[v] + 1;
					if (e.to == t) return;
					que.push(e.to);
				}
			}
		};
		auto dfs = [&](auto self, int v, Cap up) {
			if (v == s) return up;
			Cap res = 0;
			int level_v = level[v];
			for (int& i = iter[v]; i < int(g[v].size()); i++) {
				_edge& e = g[v][i];
				if (level_v <= level[e.to] || g[e.to][e.rev].cap == 0) continue;
				Cap d =
					self(self, e.to, std::min(up - res, g[e.to][e.rev].cap));
				if (d <= 0) continue;
				g[v][i].cap += d;
				g[e.to][e.rev].cap -= d;
				res += d;
				if (res == up) break;
			}
			return res;
		};

		Cap flow = 0;
		while (flow < flow_limit) {
			bfs();
			if (level[t] == -1) break;
			std::fill(iter.begin(), iter.end(), 0);
			while (flow < flow_limit) {
				Cap f = dfs(dfs, t, flow_limit - flow);
				if (!f) break;
				flow += f;
			}
		}
		return flow;
	}

	std::vector<bool> min_cut(int s) {
		std::vector<bool> visited(_n);
		internal::simple_queue<int> que;
		que.push(s);
		while (!que.empty()) {
			int p = que.front();
			que.pop();
			visited[p] = true;
			for (auto e : g[p]) {
				if (e.cap && !visited[e.to]) {
					visited[e.to] = true;
					que.push(e.to);
				}
			}
		}
		return visited;
	}

  private:
	int _n;
	struct _edge {
		int to, rev;
		Cap cap;
	};
	std::vector<std::pair<int, int>> pos;
	std::vector<std::vector<_edge>> g;
};

}  // namespace atcoder


atcoder::mf_graph<int>g(80003);
using namespace std;
int ls[80003],rs[80003],N,rt;
int n,l[5003],r[5003],I[5003];
void newnode(int &k){
	ls[++N]=ls[k];
	rs[N]=rs[k];
	g.add_edge(k+(n+1<<1),N+(n+1<<1),1e9);
	k=N;
}
void add(int &k,int _l,int _r,int p,int i){
	newnode(k);
	g.add_edge(i,k+(n+1<<1),1e9);
	if(_l==_r)return;
	if(p<=(_l+_r>>1))
		add(ls[k],_l,_l+_r>>1,p,i);
	else
		add(rs[k],_l+_r+2>>1,_r,p,i);
}
void query(int k,int _l,int _r,int l,int r,int i){
	if(l>_r||r<_l||!k)return;
	if(l<=_l&&_r<=r){g.add_edge(k+(n+1<<1),i,1e9);return;}
	query(ls[k],_l,_l+_r>>1,l,r,i);
	query(rs[k],_l+_r+2>>1,_r,l,r,i);
}
bool cmp(int i,int j){return l[i]<l[j];}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(NULL);
	cout.tie(NULL);
	cin>>n;
	for(int i=0;i<n;i++){
		cin>>l[i]>>r[i];
		l[i]--;r[i]--;I[i]=i;
		g.add_edge(n<<1,i,1);
		g.add_edge(i+n,n<<1|1,1);
	}sort(I,I+n,cmp);
	for(int i=0;i<n;i++){
		query(rt,0,n+n-1,l[I[i]],r[I[i]],i+n);
		add(rt,0,n+n-1,r[I[i]],i);
	}
	cout<<n+n-g.flow(n<<1,n<<1|1);
}

详细

Test #1:

score: 100
Accepted
time: 1ms
memory: 5620kb

input:

5
2 3
6 7
1 9
5 10
4 8

output:

9

result:

ok 1 number(s): "9"

Test #2:

score: -100
Dangerous Syscalls

input:

5000
7985 7987
42 46
1591 1593
407 410
6305 6306
1456 1457
5874 5875
7135 7137
7041 7046
6813 6815
8868 8871
665 666
4055 4056
9789 9796
7067 7068
4745 4746
5167 5171
1735 1737
2125 2128
1444 1447
1348 1352
6087 6090
1381 1384
1600 1601
5187 5190
2801 2802
8449 8450
9376 9377
4021 4024
2674 2676
490...

output:


result: