转载请注明出处:http://tokitsukaze.live/

题目链接:


题意:

6000x6000的平面上有n个点,每个点有点权。m次操作。
操作1:插入一个点。
操作2:删除一个点。
操作3:给出一个圆心(x,y),半径为sqrt(k)的圆,当前圆上存在的点,点权+w。
操作4:询问圆心(x,y),半径为sqrt(k)的圆,当前圆上存在的点的点权和。
强制在线。


题解:
由于题目给的都是整数点,在圆上的整数格点很少,所以对于操作3和操作4暴力枚举圆上的整数点即可。

那么先预处理出所有的a^2 +b^2=c^2,对应每个c^2,把符合条件的(a,b)存到vector里。

对于操作3和4,即(x-x0)^2+(y-y0)^2=k,展开|x-x0|和|y-y0|的绝对值,判断存不存在(x0,y0)即可。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
#include <bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:1024000000,1024000000")
#define mem(a,b) memset((a),(b),sizeof(a))
#define MP make_pair
#define pb push_back
#define fi first
#define se second
#define sz(x) (int)x.size()
#define all(x) x.begin(),x.end()
using namespace __gnu_cxx;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
typedef vector<int> VI;
typedef vector<ll> VL;
void go();
int main(){
#ifdef tokitsukaze
freopen("TEST.txt","r",stdin);
#endif
go();return 0;
}
const int INF=0x3f3f3f3f;
const ll LLINF=0x3f3f3f3f3f3f3f3f;
const double PI=acos(-1.0);
const double eps=1e-6;
const int MAX=1e7+10;
const ll mod=1e9+7;
/********************************* head *********************************/
vector<PII > res[MAX],cl;
void init()
{
int i,j;
for(i=0;i<=6000;i++)
{
for(j=0;j<=6000;j++)
{
if(i*i+j*j<=MAX-10) res[i*i+j*j].pb(MP(i,j));
else break;
}
}
}
bool check(int x,int y)
{
if(x<0||y<0||x>6000||y>6000) return 0;
return 1;
}
ll pre,ans;
void change(int &x,int &y)
{
x=(x+pre)%6000+1;
y=(y+pre)%6000+1;
}
bool exist[6010][6010]={0};
int mp[6010][6010]={0};
int dir[4][2]={1,1,1,-1,-1,1,-1,-1};
void go()
{
int i,t,cas=1,x,y,w,n,q,op,k,tx,ty;
init();
read(t);
while(t--)
{
cl.clear();
read(n,q);
for(i=1;i<=n;i++)
{
read(x,y,w);
exist[x][y]=1;
mp[x][y]=w;
cl.pb(MP(x,y));
}
pre=0;
set<PII > s;
printf("Case #%d:\n",cas++);
while(q--)
{
read(op,x,y);
change(x,y);
if(op==1)
{
read(w);
exist[x][y]=1;
mp[x][y]=w;
cl.pb(MP(x,y));
}
else if(op==2)
{
exist[x][y]=0;
mp[x][y]=0;
}
else if(op==3)
{
read(k,w);
s.clear();
for(auto it:res[k])
{
for(i=0;i<4;i++)
{
tx=x-it.fi*dir[i][0];
ty=y-it.se*dir[i][1];
if(!check(tx,ty)) continue;
if(exist[tx][ty]) s.insert(MP(tx,ty));
}
}
for(auto it:s) mp[it.fi][it.se]+=w;
}
else
{
read(k);
s.clear();
for(auto it:res[k])
{
for(i=0;i<4;i++)
{
tx=x-it.fi*dir[i][0];
ty=y-it.se*dir[i][1];
if(!check(tx,ty)) continue;
if(exist[tx][ty]) s.insert(MP(tx,ty));
}
}
ans=0;
for(auto it:s) ans+=mp[it.fi][it.se];
printf("%lld\n",ans);
pre=ans;
}
}
for(auto it:cl)
{
exist[it.fi][it.se]=0;
mp[it.fi][it.se]=0;
}
}
}