QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#88419 | #5668. Cell Nuclei Detection | CCSU_Wine | TL | 25ms | 98200kb | C++20 | 2.6kb | 2023-03-16 10:47:41 | 2023-03-16 10:47:45 |
Judging History
answer
#include <queue>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 2e6+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<=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[2010][2010];
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 = tot+2;
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+2;j++){
for(int z = max(b-2,0ll);z<=d+2;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: 15ms
memory: 98092kb
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: 25ms
memory: 98200kb
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 ...