QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#106363 | #6308. Magic | zhouhuanyi | ML | 95ms | 16248kb | C++11 | 3.3kb | 2023-05-17 15:23:48 | 2023-05-17 15:23:53 |
Judging History
answer
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<bitset>
#define N 10002
#define M 5001
#define inf 1e9
using namespace std;
int read()
{
char c=0;
int sum=0;
while (c<'0'||c>'9') c=getchar();
while ('0'<=c&&c<='9') sum=sum*10+c-'0',c=getchar();
return sum;
}
struct reads
{
int l,r;
bool operator < (const reads &t)const
{
return l<t.l;
}
};
reads st[N+1];
int n,s,t,ans,depth[N+1],l[N+1],r[N+1],p[N+1],tong1[N+1],num1[N+1],length1,tong2[N+1],num2[N+1],length2,len=1;
bool cl[N+1];
bitset<M+1>B[N+1];
bitset<M+1>vis[N+1];
bitset<M+1>used1;
bitset<M+1>used2;
bool cmp(int x,int y)
{
return depth[x]<depth[y];
}
void add(int x,int y)
{
if (!cl[x]) B[x][num2[y]]=1;
else B[x][num1[y]]=1;
return;
}
bool bfs()
{
int top;
queue<int>q;
bitset<M+1>S;
for (int i=s;i<=t;++i) depth[i]=0;
used1.set(),used2.set(),q.push(s),used1[num1[s]]=0,depth[s]=1;
while (!q.empty())
{
top=q.front(),q.pop();
if (!cl[top])
{
S=used2&B[top];
for (int i=S._Find_first();i!=S.size();i=S._Find_next(i)) depth[tong2[i]]=depth[top]+1,used2[i]=0,q.push(tong2[i]);
}
else
{
S=used1&B[top];
for (int i=S._Find_first();i!=S.size();i=S._Find_next(i)) depth[tong1[i]]=depth[top]+1,used1[i]=0,q.push(tong1[i]);
}
}
return !used2[num2[t]];
}
int dinic(int x,int flow)
{
if (x==t) return flow;
int k,y;
while (1)
{
if (!cl[x])
{
bitset<M+1>S=B[x]&vis[x];
y=S._Find_first();
if (y==S.size()) return 0;
k=dinic(tong2[y],1);
if (k)
{
B[x][y]=0,B[tong2[y]][num1[x]]=1;
return 1;
}
else vis[x][y]=0;
}
else
{
bitset<M+1>S=B[x]&vis[x];
y=S._Find_first();
if (y==S.size()) return 0;
k=dinic(tong1[y],1);
if (k)
{
B[x][y]=0,B[tong1[y]][num2[x]]=1;
return 1;
}
else vis[x][y]=0;
}
}
return 0;
}
int main()
{
int flow;
n=read(),ans=n<<1,s=0,t=(n<<1)+1;
for (int i=1;i<=n;++i) st[i].l=read(),st[i].r=read();
sort(st+1,st+n+1);
tong1[++length1]=s,num1[s]=length1,cl[s]=0;
for (int i=1;i<=n;++i) tong1[++length1]=i+n,num1[i+n]=length1,cl[i]=0;
tong2[++length2]=t,num2[t]=length2,cl[t]=1;
for (int i=1;i<=n;++i) tong2[++length2]=i,num2[i]=length2,cl[i]=1;
for (int i=1;i<=n;++i) add(s,i),add(i+n,t);
for (int i=1;i<=n;++i)
for (int j=i+1;j<=n;++j)
if (st[i].l<st[j].l&&st[j].l<st[i].r&&st[i].r<st[j].r)
add(i,j+n);
while (bfs())
{
for (int i=s;i<=t;++i) p[i]=i;
sort(p,p+t,cmp),used1.reset(),used2.reset();
for (int i=s;i<=t;++i) r[depth[p[i]]]=i;
for (int i=t;i>=s;--i) l[depth[p[i]]]=i;
for (int i=0;i<=t;++i)
{
if (i+1<=t)
{
for (int j=l[i+1];j<=r[i+1];++j)
{
if (!cl[p[j]]) used1[num1[p[j]]]=1;
else used2[num2[p[j]]]=1;
}
}
for (int j=l[i];j<=r[i];++j)
{
if (!cl[p[j]]) vis[p[j]]=used2;
else vis[p[j]]=used1;
}
if (i+1<=t)
{
for (int j=l[i+1];j<=r[i+1];++j)
{
if (!cl[p[j]]) used1[num1[p[j]]]=0;
else used2[num2[p[j]]]=0;
}
}
}
while (flow=dinic(s,inf)) ans-=flow;
}
printf("%d\n",ans);
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 1ms
memory: 5640kb
input:
5 2 3 6 7 1 9 5 10 4 8
output:
9
result:
ok 1 number(s): "9"
Test #2:
score: 0
Accepted
time: 16ms
memory: 16124kb
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:
8134
result:
ok 1 number(s): "8134"
Test #3:
score: 0
Accepted
time: 29ms
memory: 16196kb
input:
5000 3171 3172 4062 4064 4647 4651 3670 3673 7112 7114 9714 9717 3781 3789 8422 8426 457 460 5450 5454 7113 7122 6313 6320 9969 9973 828 832 6878 6892 4476 4483 892 903 251 259 6304 6315 130 134 9206 9215 2679 2686 9090 9091 8222 8228 9374 9375 2985 2989 3397 3401 4916 4918 6819 6821 883 889 2516 25...
output:
7047
result:
ok 1 number(s): "7047"
Test #4:
score: 0
Accepted
time: 41ms
memory: 16128kb
input:
5000 7269 7286 1979 1990 4225 4241 7866 7872 2052 2067 1508 1514 2366 2370 3488 3493 8979 8987 302 306 6730 6732 7704 7705 5528 5544 7420 7425 4705 4712 593 601 6662 6668 5228 5257 2008 2013 548 562 7949 7950 1017 1020 1025 1028 6 11 4722 4736 9945 9950 8368 8379 6781 6787 4558 4566 400 404 858 864 ...
output:
6191
result:
ok 1 number(s): "6191"
Test #5:
score: 0
Accepted
time: 50ms
memory: 16152kb
input:
5000 3005 3008 7811 7821 2832 2840 9812 9818 3947 3952 3629 3665 7443 7455 1473 1478 1467 1494 5499 5508 229 232 9477 9498 9500 9514 4769 4775 9488 9503 1514 1520 5101 5112 2455 2456 3558 3610 9072 9188 659 666 2286 2301 9735 9782 5959 5984 5823 5844 1827 1835 3658 3681 3494 3503 1016 1018 3418 3420...
output:
5636
result:
ok 1 number(s): "5636"
Test #6:
score: 0
Accepted
time: 95ms
memory: 16248kb
input:
5000 1017 1019 5731 5744 8592 8697 3414 3561 277 356 4421 4458 7969 7989 3733 3759 1975 1986 6895 6898 8580 8657 2320 2325 6494 6510 3574 3616 7721 7780 7756 7835 1744 1748 5085 5102 5428 5588 6823 6847 9348 9405 1969 2064 1152 1242 3676 3679 5569 5713 213 248 8277 8285 3739 3782 5582 5604 9076 9105...
output:
5225
result:
ok 1 number(s): "5225"
Test #7:
score: -100
Memory Limit Exceeded
input:
5000 4269 4768 7697 7804 4289 4316 6641 6742 3782 3871 9785 9858 1845 1979 6300 6412 8194 8616 3 18 8536 8790 7 129 4406 4448 3556 3609 2800 2918 3081 3164 1854 1948 8220 8546 7070 7318 1117 1195 8770 8944 5754 5760 8512 8590 5714 5829 6132 6156 8252 8381 889 1149 1899 1949 4675 4683 4282 4426 2977 ...
output:
5124