QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#88722 | #5668. Cell Nuclei Detection | buns_out | AC ✓ | 5233ms | 309544kb | C++14 | 3.4kb | 2023-03-16 21:58:46 | 2023-03-16 21:58:47 |
Judging History
answer
/*
嵌套变量重复
特殊样例 0 1 2 n
数组越界
开long long
*/
#include<bits/stdc++.h>
using namespace std;
#define endl "\n"
#define pp(x) array<int,x>
using ull=unsigned long long;
using ll=long long;
using pii=pair<int,int>;
const int mod=998244353;
const int inf=0x3f3f3f3f;
const int INF=1e9+7;
const int maxn=3e6+100;
const int N=3222;
struct Edge{
int to,next,w;
}edge[maxn];
int head[maxn],cnt=1;
int dep[maxn];
vector<int>v[N][N];
pii a1[maxn],a2[maxn];
pii b1[maxn],b2[maxn];
int s,t,tot;
int n,m;
void init(){
for(int i=0;i<=2000;i++)
for(int j=0;j<=2000;j++)
v[i][j].clear();
memset(head+1,0,sizeof(int)*tot);
cnt=1;
tot=0;
}
void add(int from,int to,int w)
{
edge[++cnt].w=w;
edge[cnt].to=to;
edge[cnt].next=head[from];
head[from]=cnt;
}
bool bfs(int s,int t)
{
memset(dep+1,0,sizeof(int)*tot);
queue<int>q;
dep[s]=1;
q.push(s);
while(!q.empty())
{
int x=q.front();q.pop();
if(x==t)return true;
for(int i=head[x];i;i=edge[i].next)
{
int y=edge[i].to;
int w=edge[i].w;
if(dep[y]==0&&w>0)
q.push(y),dep[y]=dep[x]+1;
}
}
return false;
}
int dfs(int x,int flow,int t)
{
if(x==t)return flow;
int out=flow;
for(int i=head[x];i;i=edge[i].next)
{
int y=edge[i].to;
int w=edge[i].w;
if(w!=0&&dep[y]==dep[x]+1)
{
int tmp=dfs(y,min(out,w),t);
edge[i].w-=tmp;
edge[i^1].w+=tmp;
out-=tmp;
if(!out)break;
}
}
if(out==flow)dep[x]=0;
return flow-out;
}
int DINIC(int s,int t)
{
int ans=0;
while(bfs(s,t))
ans+=dfs(s,INF,t);
return ans;
}
int in1[maxn];
int in2[maxn];
void solve()
{
init();
cin>>n>>m;
s=++tot;
t=++tot;
for(int i=1;i<=n;i++)
{
in1[i]=++tot;
// cout<<"I "<<i<<" id "<<in1[i]<<endl;
add(in1[i],t,1);
add(t,in1[i],0);
cin>>a1[i].first>>a1[i].second>>a2[i].first>>a2[i].second;
for(int x=a1[i].first;x<=a2[i].first;x++)
for(int y=a1[i].second;y<=a2[i].second;y++)
v[x][y].push_back(i);
}
int x,y,xx,yy;
auto check=[&](int i)->bool{
// cout<<x<<" "<<y<<" "<<xx<<" "<<yy<<endl;
// cout<<a1[i].first<<" "<<a1[i].second<<" "<<a2[i].first<<" "<<a2[i].second<<endl;
int x1=max(x,a1[i].first);
int x2=min(xx,a2[i].first);
int y1=max(y,a1[i].second);
int y2=min(yy,a2[i].second);
// cout<<x1<<" "<<y1<<" "<<x2<<" "<<y2<<endl;
int S1=(a2[i].first-a1[i].first)*(a2[i].second-a1[i].second);
int S2=(x2-x1)*(y2-y1);
// cout<<S1<<" "<<S2<<endl;
// cout<<(2*S2>=S1)<<endl;
// cout<<endl;
if(x1>=x2||y1>=y2)return 0;
return 2*S2>=S1;
};
for(int i=1;i<=m;i++)
{
in2[i]=++tot;
// cout<<"J "<<i<<" id "<<in2[i]<<endl;
add(s,in2[i],1);
add(in2[i],s,0);
cin>>x>>y>>xx>>yy;
map<int,bool>mp;
for(int dx=x;dx<=xx;dx++)
for(int dy=y;dy<=yy;dy++)
for(auto j:v[dx][dy])
if(check(j)==1&&!mp.count(j)){
add(in2[i],in1[j],1);
add(in1[j],in2[i],0);
// cout<<"edge "<<in2[i]<<" "<<in1[j]<<endl;
mp[j]=1;
}
}
cout<<DINIC(s,t)<<endl;
}
signed main(){
// freopen("C:\\Users\\tob\\Desktop\\P1000_0.in.txt","r",stdin);
// freopen("C:\\Users\\tob\\Desktop\\P1000_0.out.txt","w",stdout);
ios::sync_with_stdio(false);
cin.tie(nullptr);cout.tie(nullptr);
int __;cin>>__;
while(__--)
solve();
// fclose(stdin);
// fclose(stdout);
return 0;
}
/*
2 2
1 1 3 3
3 3 5 5
2 2 4 4
4 4 6 6
*/
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 59ms
memory: 248064kb
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: 68ms
memory: 247604kb
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: 0
Accepted
time: 1444ms
memory: 298324kb
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:
50000 50000 0 50000 3150
result:
ok 5 lines
Test #4:
score: 0
Accepted
time: 516ms
memory: 309544kb
input:
5 50000 50000 0 0 1 1 1 0 2 1 2 0 3 1 3 0 4 1 4 0 5 1 5 0 6 1 6 0 7 1 7 0 8 1 8 0 9 1 9 0 10 1 10 0 11 1 11 0 12 1 12 0 13 1 13 0 14 1 14 0 15 1 15 0 16 1 16 0 17 1 17 0 18 1 18 0 19 1 19 0 20 1 20 0 21 1 21 0 22 1 22 0 23 1 23 0 24 1 24 0 25 1 25 0 26 1 26 0 27 1 27 0 28 1 28 0 29 1 29 0 30 1 30 0 ...
output:
50000 25050 12500 16000 8000
result:
ok 5 lines
Test #5:
score: 0
Accepted
time: 474ms
memory: 277356kb
input:
5 50000 50000 0 0 2 4 4 0 7 1 8 0 10 1 12 0 15 3 16 0 19 1 20 0 22 2 24 0 26 4 28 0 30 4 32 0 36 3 36 0 40 1 40 0 44 1 44 0 47 2 48 0 49 3 52 0 54 1 56 0 59 4 60 0 64 3 64 0 68 3 68 0 70 1 72 0 76 4 76 0 80 3 80 0 84 4 84 0 87 2 88 0 90 1 92 0 94 4 96 0 98 1 100 0 104 1 104 0 107 2 108 0 110 4 112 0...
output:
10594 10779 10618 10381 10779
result:
ok 5 lines
Test #6:
score: 0
Accepted
time: 5233ms
memory: 274992kb
input:
5 50000 50000 0 0 4 4 1 0 5 4 2 0 6 4 3 0 7 4 4 0 8 4 5 0 9 4 6 0 10 4 7 0 11 4 8 0 12 4 9 0 13 4 10 0 14 4 11 0 15 4 12 0 16 4 13 0 17 4 14 0 18 4 15 0 19 4 16 0 20 4 17 0 21 4 18 0 22 4 19 0 23 4 20 0 24 4 21 0 25 4 22 0 26 4 23 0 27 4 24 0 28 4 25 0 29 4 26 0 30 4 27 0 31 4 28 0 32 4 29 0 33 4 30...
output:
50000 50000 50000 50000 49600
result:
ok 5 lines
Test #7:
score: 0
Accepted
time: 68ms
memory: 247912kb
input:
1 4 4 1 1 3 3 2 1 4 3 1 2 3 4 2 2 4 4 2 1 4 3 3 2 5 4 1 2 3 4 2 3 4 5
output:
3
result:
ok single line: '3'