[曾经的Ynoi2012] NOIP2015洋溢着希望

人生第一道$Ynoi$当然要写博客纪念一下啦$qwq$

---

update 2020.4.9:这玩意已经被lxl从ynoi除名了真香

---

这肯定是$Ynoi$里面最水的一个了
题意就是让你维护区间$\sin$和
那么根据三角函数的公式可以得到
$$\begin{gathered} \sin (\alpha +\beta )=\sin \alpha \cos \beta +\cos \alpha \sin \beta \\ \cos (\alpha +\beta )=\cos \alpha \cos \beta -\sin \alpha \sin \beta \end{gathered}$$
然后分别维护区间$\sin$,$\cos$就可以了
然后经过若干卡常就过了。。。
有一个显然的优化是下传标记的时候先把$tag$的$\sin$和$\cos$算出来然后再更新，这样可以减少一半时间
记得开longlong

#include <bits/stdc++.h>
using namespace std;

# define Rep(i,a,b) for(int i=a;i<=b;i++)
# define _Rep(i,a,b) for(int i=a;i>=b;i--)
# define RepG(i,u) for(int i=head[u];~i;i=e[i].next)

typedef long long ll;

const int N=2e5+5;

template<typename T> void read(T &x){
   x=0;int f=1;
   char c=getchar();
   for(;!isdigit(c);c=getchar())if(c=='-')f=-1;
   for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+c-'0';
    x*=f;
}

int n,m;
ll a[N];

struct segment_tree{
    int l,r;
    double sin,cos;
    ll tag;
}seg[N<<2];

# define lc (u<<1)
# define rc (u<<1|1)

void pushup(int u){
    seg[u].sin=seg[lc].sin+seg[rc].sin;
    seg[u].cos=seg[lc].cos+seg[rc].cos;
}

void renew(int u,ll k){
    double lstsin=seg[u].sin,lstcos=seg[u].cos,nsin=sin(k),ncos=cos(k);
    seg[u].sin=lstsin*ncos+lstcos*nsin;
    seg[u].cos=lstcos*ncos-lstsin*nsin;
}

void pushdown(int u){
    renew(lc,seg[u].tag);
    renew(rc,seg[u].tag);
    seg[lc].tag+=seg[u].tag;
    seg[rc].tag+=seg[u].tag;
    seg[u].tag=0;
}

void build(int u,int l,int r){
    seg[u].l=l,seg[u].r=r;
    if(l==r){
        seg[u].sin=sin(a[l]);
        seg[u].cos=cos(a[l]);
        return;
    }
    int mid=l+r>>1;
    build(lc,l,mid);
    build(rc,mid+1,r);
    pushup(u);
}

void update(int u,int l,int r,ll k){
    if(seg[u].l>=l&&seg[u].r<=r){
        seg[u].tag+=k;
        renew(u,k);
        return;
    }
    if(seg[u].tag)pushdown(u);
    int mid=seg[u].l+seg[u].r>>1;
    if(l<=mid)update(lc,l,r,k);
    if(r>mid)update(rc,l,r,k);
    pushup(u);
}

double query(int u,int l,int r){
    if(seg[u].l>=l&&seg[u].r<=r)return seg[u].sin;
    if(seg[u].tag)pushdown(u);
    int mid=seg[u].l+seg[u].r>>1;
    double res=0;
    if(l<=mid)res+=query(lc,l,r);
    if(r>mid)res+=query(rc,l,r);
    return res;
}

int main()
{
    read(n);
    Rep(i,1,n)read(a[i]);
    read(m);
    build(1,1,n);
    Rep(i,1,m){
        int opt,x,y;ll k;
        read(opt),read(x),read(y);
        if(opt==1)read(k),update(1,x,y,k);
        else printf("%.1lf\n",query(1,x,y));
    }
    return 0;
}

据说程序长度刚好是100行有助于卡常呢