QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#88414#5668. Cell Nuclei DetectionCCSU_WineTL 27ms106748kbC++202.6kb2023-03-16 10:30:222023-03-16 10:30:23

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-03-16 10:30:23]
  • 评测
  • 测评结果:TL
  • 用时:27ms
  • 内存:106748kb
  • [2023-03-16 10:30:22]
  • 提交

answer

#include <queue>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 5e6+7;
typedef long long ll;
#define int long long 
const int inf = 0x3f3f3f3f;
int n,m,s,t;
struct edge{
	int to,next,now,w;
}e[N];
int head[N],ecnt,d[N],vis[N],cur[N];
int dist[N];
void add(int u,int v,int w){
	e[ecnt] = {v,head[u],0,w};
	head[u] = ecnt++;
}
bool bfs(int s,int t){
	for(int i = 0;i<=2*t;i++) d[i] = vis[i] = 0;
	queue<int> q;
	q.push(s);
	vis[s] = 1;
	while(!q.empty()){
		int u = q.front();
		q.pop();
		for(int i = head[u];~i;i = e[i].next){
			int to = e[i].to,now = e[i].now,w = e[i].w;
			if(!vis[to]&&now < w){
				vis[to] = 1;
				d[to] = d[u] + 1;
				q.push(to);
				if(to == t) return 1;
			}
		}
	}
	return vis[t];
}
ll dfs(int u,int t,int imp){
	if(u == t || !imp) return imp;
	ll res = 0;
	for(int &i = cur[u];~i;i=e[i].next){
		int to = e[i].to;
		int now = e[i].now;
		if(d[u]+1!=d[to]) continue;
		int f = dfs(to,t,min(imp,e[i].w-now)) ;
		if(f==0) continue;
		e[i].now += f;
		e[i^1].now -= f;
		res += f;
		imp -= f;
		if(!imp) break; 
	}
	return res;
}
int k;
ll Dinic(int s,int t){
	ll flow = 0;
	while(bfs(s,t)){
		for(int i = 0;i<=k;i++){
			cur[i] = head[i]; 
		}
		flow += dfs(s,t,inf);
	}
	return flow;
}
int x[N],y[N],xx[N],yy[N];
vector<int> G[2100][2100];
int cal(int x,int y,int xx,int yy){
	return (yy-y)*(xx-x);
}
bool check(int x,int y,int xx,int yy,int a,int b,int c,int d){
	if(x==xx&&y==yy){
		return 0;
	}
	int mx = max(x,a);
	int my = max(y,b);
	int mxx = min(xx,c);
	int myy = min(yy,d);
	if(mxx<=mx||myy<=my) return 0;	
	return cal(mx,my,mxx,myy)*2 >= cal(x,y,xx,yy);
}
signed main(){
ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr);
	int T;
	cin>>T;
	while(T--){
		int num1,num2;
		cin>>num1>>num2;
		int tot = num1+num2;
		m = num2*40; 
		k = m;
		ecnt = 0;
		for(int i = 0;i<=m;i++){
			head[i] = -1;
		}
		int s = 0,t = tot + 1;
		for(int i = 1;i<=num1;i++){
			cin>>x[i]>>y[i]>>xx[i]>>yy[i];
			G[x[i]][y[i]].push_back(i);
		} 
		for(int i = num1+1;i<=num2+num1;i++){
			int a,b,c,d;
			cin>>a>>b>>c>>d;
			for(int j = max(a-2,0ll);j<=c;j++){
				for(int z = max(b-2,0ll);z<=d;z++){
					if(G[j][z].size()==0) continue;
					for(int ix : G[j][z]){
						if(check(x[ix],y[ix],xx[ix],yy[ix],a,b,c,d)){
							add(i,ix,1);
							add(ix,i,0);
						}
					}		
				}
			}
		}
		for(int i = 1;i<=num1;i++){
			add(i,t,1);
			add(t,i,0);
			G[x[i]][y[i]].clear();
		} 
		for(int i = num1+1;i<t;i++){
			add(s,i,1);
			add(i,s,0);
		}
		cout << Dinic(s,t) << '\n';
	}
	return 0;
}
//

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 27ms
memory: 106748kb

input:

3
2 2
1 1 3 3
3 3 5 5
2 2 4 4
4 4 6 6
2 3
1 1 3 3
3 3 5 5
1 3 3 5
2 1 4 5
3 1 5 3
3 3
1 1 2 2
2 2 3 3
3 3 4 4
1 1 3 3
2 2 4 4
3 3 5 5

output:

0
1
3

result:

ok 3 lines

Test #2:

score: 0
Accepted
time: 27ms
memory: 106748kb

input:

3
2 2
1 1 3 3
3 3 5 5
2 2 4 4
4 4 6 6
2 3
1 1 3 3
3 3 5 5
1 3 3 5
2 1 4 5
3 1 5 3
3 3
1 1 2 2
2 2 3 3
3 3 4 4
1 1 3 3
2 2 4 4
3 3 5 5

output:

0
1
3

result:

ok 3 lines

Test #3:

score: -100
Time Limit Exceeded

input:

5
50000 50000
0 0 4 4
4 0 8 4
8 0 12 4
12 0 16 4
16 0 20 4
20 0 24 4
24 0 28 4
28 0 32 4
32 0 36 4
36 0 40 4
40 0 44 4
44 0 48 4
48 0 52 4
52 0 56 4
56 0 60 4
60 0 64 4
64 0 68 4
68 0 72 4
72 0 76 4
76 0 80 4
80 0 84 4
84 0 88 4
88 0 92 4
92 0 96 4
96 0 100 4
100 0 104 4
104 0 108 4
108 0 112 4
112 ...

output:


result: