QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#131672#5463. Range Closest Pair of Points Querysumi007TL 6411ms247460kbC++141.9kb2023-07-27 20:18:152023-07-27 20:18:17

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-07-27 20:18:17]
  • 评测
  • 测评结果:TL
  • 用时:6411ms
  • 内存:247460kb
  • [2023-07-27 20:18:15]
  • 提交

answer

#include<bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
using namespace std;
using namespace __gnu_pbds;
#define ll long long
#define i64 __int128
#define fi first
#define se second
#define pii pair<int,int>
#define pb push_back
#define d(i,j) 1ll*j*998244353+i
const ll N = 3e5+5,inf = 1e18;
ll n,Q,B,nodex,mn[N],mn_B[N],b[N],R[N],ans[N],L[N];
int dx[10]={-1,-1,-1,0,0,0,1,1,1},dy[10]={-1,0,1,-1,0,1,-1,0,1};
struct node{
	ll x,y;
}p[N];
vector<pii> q[N];
vector<int> pos[N*30];
gp_hash_table<ll,int> num[30];
ll calc(ll sx,ll sy,ll ex,ll ey){
	return (sx-ex)*(sx-ex)+(sy-ey)*(sy-ey);
}
ll ask(ll l,ll r){
	ll res = inf;
	for(int i=l;i<=min(R[b[l]],r);i++) res = min(res,mn[i]);
	for(int i=b[l]+1;i<=b[r];i++){
		if(mn_B[i]==inf) break;
		res = min(res,mn_B[i]);
	}
	return res;
}
int main(){
	ios::sync_with_stdio(0);
	scanf("%d %d",&n,&Q);
	B = sqrt(n);
	for(int i=1;i<=n;i++){
		scanf("%d %d",&p[i].x,&p[i].y);
		b[i] = (i-1)/B+1;
		R[b[i]] = i;
		if(b[i]!=b[i-1]) L[b[i]] = i;
	}
	for(int i=1,l,r;i<=Q;i++){
		scanf("%d %d",&l,&r);
		q[r].pb({l,i});
	}
	for(int i=1;i<=n;i++) mn[i] = mn_B[i] = inf;
	for(int i=1;i<=n;i++){
		for(int j=0;(1ll<<j)<=1e8;j++){
			ll tx = (p[i].x)>>j,ty = (p[i].y)>>j,lim = (int)(1e8)>>j;
			if(num[j].find(d(tx,ty))==num[j].end()) num[j][d(tx,ty)] = ++nodex;
			for(int k=0;k<9;k++){
				ll nx = tx+dx[k],ny = ty+dy[k]; 
				if(nx<0 || ny<0 || nx>lim || ny>lim || num[j].find(d(nx,ny))==num[j].end()) continue;
				int id = num[j][d(nx,ny)];
				vector<int> tmp;
				for(int l:pos[id]){
					ll dis = calc(p[l].x,p[l].y,p[i].x,p[i].y);
					mn[l] = min(mn[l],dis),mn_B[b[l]] = min(mn_B[b[l]],dis);
					if(dis>=(1ll<<(2*j-2))) tmp.pb(l);
				}
				pos[id] = tmp;
				if(k==4) pos[id].pb(i);
			}
		}
		for(pii t:q[i]) ans[t.se] = ask(t.fi,i);
	}
	for(int i=1;i<=Q;i++) printf("%lld\n",ans[i]);
	return 0;
}

详细

Test #1:

score: 100
Accepted
time: 27ms
memory: 234436kb

input:

5 5
2 4
1 1
3 3
5 1
4 2
1 5
2 3
2 4
3 5
1 3

output:

2
8
8
2
2

result:

ok 5 number(s): "2 8 8 2 2"

Test #2:

score: 0
Accepted
time: 24ms
memory: 234416kb

input:

2 1
1 1
1 1
1 2

output:

0

result:

ok 1 number(s): "0"

Test #3:

score: 0
Accepted
time: 27ms
memory: 235772kb

input:

2 1
1 100000000
100000000 1
1 2

output:

19999999600000002

result:

ok 1 number(s): "19999999600000002"

Test #4:

score: 0
Accepted
time: 140ms
memory: 240812kb

input:

20000 250000
3 10
5 4
4 7
1 5
2 1
10 6
2 3
8 4
2 1
8 5
9 8
7 7
4 5
2 7
9 4
9 10
3 2
9 5
10 2
9 2
3 1
9 9
6 5
9 5
9 10
9 1
1 2
8 8
3 4
7 6
6 2
6 8
6 6
8 4
10 2
1 1
10 2
8 3
4 4
5 5
5 1
4 9
7 6
6 8
6 4
1 6
10 3
3 2
4 10
6 8
9 7
2 10
7 8
10 7
3 2
5 1
6 4
7 9
1 3
4 9
4 8
9 4
5 2
2 2
9 2
9 2
9 6
6 9
8 7
...

output:

0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
...

result:

ok 250000 numbers

Test #5:

score: 0
Accepted
time: 1397ms
memory: 239420kb

input:

20000 250000
72 45
72 34
34 10
20 96
79 39
43 5
72 49
56 85
1 72
44 70
73 89
69 76
49 89
57 38
39 9
33 47
22 3
96 41
90 82
25 6
72 92
73 38
53 21
16 88
59 9
54 2
14 6
7 94
99 68
27 82
70 50
81 81
60 81
2 98
33 19
98 9
35 36
49 66
86 7
3 95
32 89
62 42
68 88
65 16
94 6
85 10
51 69
90 36
70 87
13 79
4...

output:

0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0...

result:

ok 250000 numbers

Test #6:

score: 0
Accepted
time: 6411ms
memory: 247460kb

input:

20000 250000
116 165
150 677
951 676
552 324
267 222
739 755
705 663
235 142
152 714
735 641
414 201
313 663
683 300
403 739
37 521
42 833
553 733
886 449
86 913
55 637
731 932
143 161
500 891
719 79
916 237
431 992
804 210
643 332
165 362
178 332
821 510
762 34
188 83
283 357
743 407
750 487
19 658...

output:

0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
45
1
0
1
0
0
0
2
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
4
0
0
0
0
0
0
0
0
0
2
0
0
0
0
0
0
0
0
0
0
52
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
2
0
0
0
0
5
0
0
0
0
0
1
0
0
0
1
1
0
0
0
2
0
0
0
0
0
0
...

result:

ok 250000 numbers

Test #7:

score: -100
Time Limit Exceeded

input:

20000 250000
193144 901630
894860 3728
802840 155641
625273 792794
433205 942335
223506 874548
425235 550629
341169 950916
49296 305542
709463 512381
9742 994078
106664 845553
38691 373010
184728 331946
344567 438182
854583 291040
94405 555036
56330 494420
928479 123077
796593 798314
300374 490072
2...

output:


result: