#line 2 "library/other/template.hpp" #define _CRT_SECURE_NO_WARNINGS #pragma target("avx2") #pragma optimize("O3") #pragma optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i,n) for(int i=0;i<(n);i++) #define REP(i,n) for(int i=1;i<=(n);i++) #define all(V) V.begin(),V.end() typedef unsigned int uint; typedef long long lint; typedef unsigned long long ulint; typedef std::pair P; typedef std::pair LP; constexpr int INF = INT_MAX/2; constexpr lint LINF = LLONG_MAX/2; constexpr double eps = DBL_EPSILON; constexpr double PI=3.141592653589793238462643383279; template class prique :public std::priority_queue, std::greater> {}; template inline bool chmax(T& lhs, const U& rhs) { if (lhs < rhs) { lhs = rhs; return 1; } return 0; } template inline bool chmin(T& lhs, const U& rhs) { if (lhs > rhs) { lhs = rhs; return 1; } return 0; } inline lint gcd(lint a, lint b) { while (b) { lint c = a; a = b; b = c % b; } return a; } inline lint lcm(lint a, lint b) { return a / gcd(a, b) * b; } bool isprime(lint n) { if (n == 1)return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0)return false; } return true; } template T mypow(T a, lint b) { T res(1); while(b){ if(b&1)res*=a; a*=a; b>>=1; } return res; } lint modpow(lint a, lint b, lint m) { lint res(1); while(b){ if(b&1){ res*=a;res%=m; } a*=a;a%=m; b>>=1; } return res; } template void printArray(std::vector& vec) { rep(i, vec.size()){ std::cout << vec[i]; std::cout<<(i==(int)vec.size()-1?"\n":" "); } } template void printArray(T l, T r) { T rprev = std::prev(r); for (T i = l; i != rprev; i++) { std::cout << *i << " "; } std::cout << *rprev << std::endl; } LP extGcd(lint a,lint b) { if(b==0)return {1,0}; LP s=extGcd(b,a%b); std::swap(s.first,s.second); s.second-=a/b*s.first; return s; } LP ChineseRem(const lint& b1,const lint& m1,const lint& b2,const lint& m2) { lint p=extGcd(m1,m2).first; lint tmp=(b2-b1)*p%m2; lint r=(b1+m1*tmp+m1*m2)%(m1*m2); return std::make_pair(r,m1*m2); } template inline constexpr decltype(auto) lambda_fix(F&& f){ return [f=std::forward(f)](auto&&... args){ return f(f,std::forward(args)...); }; } #line 3 "library/data-structure/BIT.hpp" template class BIT { int n; std::vector bit; private: T query(int a) { T cnt = 0; while (a > 0) { cnt += bit[a]; a -= a & -a; } return cnt; } public: BIT(int n) :n(n) {bit.resize(n + 1);} void add(int a, T x) { a++; while (a <= n) { bit[a] += x; a += a & -a; } } T query(int l, int r) {return query(r) - query(l);} void clear() {bit.assign(n + 1, 0);} int lower_bound(T x){ int p=0,k=1; while(k*2<=n)k*=2; while(k>0){ if(p+k<=n&&bit[p+k]0){ if(p+k<=n&&bit[p+k]<=x){ x-=bit[p+k]; p+=k; } k/=2; } return p; } }; #line 3 "main.cpp" int n,q,b,d[100010]; double angle[100010],bangle[2010],r[2010]; void update(int bid){ double nx=0,ny=0,nowa=0; for(int j=bid*b;j>n>>q; b=std::sqrt(n); std::fill(d,d+n,1); std::fill(r,r+n/b+10,b); BIT anglebit(n); rep(i,q){ int t; std::cin>>t; if(t==0){ int i; double x; std::cin>>i>>x; i--; angle[i]=x*PI/180; anglebit.add(i,angle[i]-anglebit.query(i,i+1)); if(i+b>i>>x; i--; d[i]=x; if(i+b>i; int now=0; double nx=0,ny=0,na=0; while(now+b