2018-icpc沈阳-G-思维

http://codeforces.com/gym/101955/problem/G

　　给出一个6000*6000的坐标系，有四种操作，一是加入放置一个点到某个空格子上，二是从某个有点的格子移走一个点，三是将距离(x,y)距离为根号k的点权值加上w，四是询问距离(x,y)距离为根号k的点权值总和。

　　由于k都是整数，而在圆上的整数点很少，所以想到，A^2+B^2=K^2，处理出所有(A,B)对于每个A^2+B^2。1，2操作就很简单了，3，4操作的话直接暴力从K对应的(A,B)暴力查找合法的点。

　　各种TLE,WA，这个C数组大小6000*6000，每次都memset的话就会T，只能用一个vector记录下本次涉及到的点最后清零。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<map>
 #include<set>
 #include<stack>
 #include<deque>
 #include<bitset>
 #include<unordered_map>
 #include<unordered_set>
 #include<queue>
 #include<cstdlib>
 #include<ctype.h>
 #include<ctime>
 #include<functional>
 #include<algorithm>
 #include<bits/stdc++.h>
 using namespace std;
 #define LL long long
 #define pii pair<int,int>
 #define mp make_pair
 #define pb push_back
 #define fi first
 #define se second
 #define inf 0x3f3f3f3f
 #define debug puts("debug")
 #define mid ((L+R)>>1)
 #define lc (id<<1)
 #define rc (id<<1|1)
 const int maxn=;
 const int maxm=;
 const double PI=acos(-1.0);
 const double eps=1e-;
 const LL mod=1e9+;
 LL gcd(LL a,LL b){return b==?a:gcd(b,a%b);}
 LL lcm(LL a,LL b){return a/gcd(a,b)*b;}
 LL qpow(LL a,LL b,LL c){LL r=; for(;b;b>>=,a=a*a%c)if(b&)r=r*a%c;return r;}
 struct Edge{int v,w,next;};
 
 template<class T>
 ostream & operator<<(ostream &out,vector<T>&v){
     for(auto x:v)cout<<x<<' ';
     return out;
 }
 void read(LL &n){
     n=; char c=getchar();
     while(c<''||c>'')c=getchar();
     while(c>=''&&c<='') n=(n<<)+(n<<)+(c-''),c=getchar();
 }
 LL C[][];
 LL dir[][]={,,,-,-,,-,-};
 bool ojbk(LL x,LL y){
     if(x<||x>||y<||y>||C[x][y]==)return ;
     return ;
 }
 #define pll pair<LL,LL>
 const LL MAX=+;
 vector<pll> cl,a[MAX];
 void AC(){
     LL T,n,m,cas=;
     scanf("%lld",&T);
     while(T--){
         scanf("%lld%lld",&n,&m);
         LL last=;
         LL x,y,w,k,op;
         for(int i=;i<=n;++i){
             scanf("%lld%lld%lld",&x,&y,&w);
             cl.pb(mp(x,y));
             C[x][y]+=w;
         }
         printf("Case #%lld:\n",++cas);
         while(m--){
             scanf("%lld%lld%lld",&op,&x,&y);
             x=(x+last)%+;
             y=(y+last)%+;cl.pb(mp(x,y));
             if(op==){
                 scanf("%lld",&w);
                 C[x][y]=w;
 
             }
             else if(op==){
                 C[x][y]=;
             }
             else if(op==){
                 scanf("%lld%lld",&k,&w);
                 set<pll>S;
                 for(auto v:a[k]){
                     for(int i=;i<;++i){
                         LL sx=x+v.fi*dir[i][],sy=y+v.se*dir[i][];
                         if(ojbk(sx,sy))S.insert(mp(sx,sy));
                     }
 
                 }for(auto o:S){
                         C[o.fi][o.se]+=w;
                     }
             }
             else{
                 scanf("%lld",&k);
                 LL ans=;
                 set<pll>S;
                 for(auto v:a[k]){
 
                     for(int i=;i<;++i){
                     LL sx=x+v.fi*dir[i][],sy=y+v.se*dir[i][];
                         if(ojbk(sx,sy))S.insert(mp(sx,sy));
                     }
 
                 }for(auto o:S){
                         ans+=C[o.fi][o.se];
                     }
                 printf("%lld\n",ans);
                 last=ans;
             }
         }
         for(auto v:cl)C[v.fi][v.se]=;
         cl.clear();
     }
 }
 int main(){
     for(LL i=;i<=;++i){
         for(LL j=;j<=;++j){
             if(i*i+j*j>MAX-)break;
             a[i*i+j*j].pb(mp(i,j));
         }
     }
     AC();
     return ;
 }
 /*
 
 1
 3 6
 2999 3000 1
 3001 3000 1
 3000 2999 1
 1 2999 3000 1
 4 2999 2999 1
 2 2995 2996
 3 2995 2995 1 1
 4 2995 2995 1
 4 3000 3000 1
 
 */