QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#106380 | #6308. Magic | zhouhuanyi | ML | 100ms | 16380kb | C++11 | 3.2kb | 2023-05-17 15:49:17 | 2023-05-17 15:49:19 |
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,dque[N+1],sl,sr,leng,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;
bitset<M+1>S;
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;
for (int i=s;i<=t;++i) depth[i]=0;
used1.set(),used2.set(),dque[sl=sr=1]=s,used1[num1[s]]=0,depth[s]=1;
while (sl<=sr)
{
top=dque[sl],sl++;
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,dque[++sr]=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,dque[++sr]=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])
{
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
{
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: 2ms
memory: 5584kb
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: 22ms
memory: 16196kb
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: 26ms
memory: 16364kb
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: 42ms
memory: 16084kb
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: 16096kb
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: 100ms
memory: 16244kb
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: 0
Accepted
time: 79ms
memory: 16380kb
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
result:
ok 1 number(s): "5124"
Test #8:
score: -100
Memory Limit Exceeded
input:
5000 1247 1337 3192 3474 7887 8008 2912 3065 4356 4573 9268 9376 987 1165 6972 7155 580 668 331 718 455 640 4815 4854 3783 3988 8426 8579 7459 7563 1991 2482 9585 9973 7562 7627 7472 8428 86 320 8319 8331 8596 8953 5729 6468 7209 7364 1026 1246 6575 6608 5048 5050 1957 2531 9408 9489 3092 3096 845 1...
output:
5079