QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#186231 | #5668. Cell Nuclei Detection | BUET_POTATOES# | AC ✓ | 5660ms | 323140kb | C++20 | 3.6kb | 2023-09-23 14:30:29 | 2023-09-23 14:30:29 |
Judging History
answer
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt")
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int nmax = 2e3+2;
vector<int>bam[nmax][nmax];
vector<int>dan[nmax][nmax];
const int N = 5e4+5;
int X1[2][N], X2[2][N], Y1[2][N], Y2[2][N];
int interSecc(int L1, int R1, int L2, int R2){
int L = max(L1, L2), R = min(R1, R2);
if(L<=R) return R-L;
return 0;
}
bool hasEdge(int i, int j){
int aria = (X2[0][i] - X1[0][i]) * (Y2[0][i] - Y1[0][i]);
int interAria = interSecc(X1[0][i], X2[0][i], X1[1][j], X2[1][j]) * interSecc(Y1[0][i], Y2[0][i], Y1[1][j], Y2[1][j]);
if(interAria*2 >= aria){
return true;
}
return false;
}
namespace HopcroftKarp { // 1-indexed. L,R indep
const int maxN=1e6+7,maxM=1e6+7; int n, m, match;
int vis[maxN], level[maxN], ml[maxN], mr[maxM];
vector<int> edge[maxN];
void init(int N, int M) { //N=left nodes, M = right
n = N, m = M;
for (int i=1;i<=n;i++) edge[i].clear(),ml[i]=-1;
for (int i=1;i<=m;i++) mr[i] = -1;
match = 0;
}
void add(int u, int v) {
edge[u].push_back(v);
if(ml[u]==-1&&mr[v]==-1)ml[u]=v,mr[v]=u,match++;
}
bool dfs(int u) {
vis[u] = true;
for (int x: edge[u]) {
int v = mr[x];
if (v == -1 || (!vis[v]
&& level[u] < level[v] && dfs(v))) {
ml[u] = x; mr[x] = u; return true;
}
}
return false;
}
int matching() {
while (true) {
queue<int> q;
for (int i = 1; i <= n; ++i) {
if (ml[i] == -1) level[i] = 0, q.push(i);
else level[i] = -1;
}
while (!q.empty()) {
int u = q.front(); q.pop();
for (int x: edge[u]) {
int v = mr[x];
if (v != -1 && level[v] < 0) {
level[v] = level[u] + 1; q.push(v);
}
}
}
for (int i = 1; i <= n; ++i) vis[i] = false;
int d = 0;
for(int i=1;i<=n;++i)if(ml[i]==-1&&dfs(i))++d;
if (d == 0) return match; match += d;
}
} }
void solv(int ces){
int m,n;
cin>>m>>n;
for(int i=1; i<=m; i++){
int x1,y1,x2,y2;
cin>>x1>>y1>>x2>>y2;
X1[0][i] = x1;
X2[0][i] = x2;
Y1[0][i] = y1;
Y2[0][i] = y2;
for(int x=x1; x<=x2; x++){
for(int y=y1; y<=y2; y++){
bam[x][y].push_back(i);
}
}
}
for(int i=1; i<=n; i++){
int x1,y1,x2,y2;
cin>>x1>>y1>>x2>>y2;
X1[1][i] = x1;
X2[1][i] = x2;
Y1[1][i] = y1;
Y2[1][i] = y2;
for(int x=x1; x<=x2; x++){
for(int y=y1; y<=y2; y++){
dan[x][y].push_back(i);
}
}
}
set<pii>allP;
for(int i=0; i<nmax; i++){
for(int j=0; j<nmax; j++){
for(int a : bam[i][j]){
for(int b : dan[i][j]){
if(hasEdge(a, b)){
allP.emplace(a, b);
}
}
}
bam[i][j].clear();
dan[i][j].clear();
}
}
// sort(allP.begin(), allP.end());
// allP.erase(unique(allP.begin(), allP.end()), allP.end());
HopcroftKarp::init(m, n);
for(auto [i,j] : allP){
// if(hasEdge(i, j)){
HopcroftKarp::add(i, j);
// }
}
cout<< HopcroftKarp::matching() <<"\n";
}
int main(){
ios::sync_with_stdio(false);
cin.tie(NULL);
int tc = 1;
cin>>tc;
for(int ces=1; ces<=tc; ces++){
solv(ces);
}
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 82ms
memory: 225960kb
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: 101ms
memory: 224100kb
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: 1770ms
memory: 323140kb
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: 942ms
memory: 311904kb
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: 476ms
memory: 282608kb
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: 5660ms
memory: 272316kb
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: 53ms
memory: 226108kb
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'