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

人生第一道$Ynoi$當然要寫博客紀念一下啦$qwq$

---

update 2020.4.9:這玩意已經被lxl從ynoi除名了真香

---

這肯定是$Ynoi$裏面最水的一個了
題意就是讓你維護區間$\sin$和
那麼根據三角函數的公式可以得到
$\sin(\alpha+\beta)=\sin\alpha\cos\beta+\cos\alpha\sin\beta\\ \cos(\alpha+\beta)=\cos\alpha\cos\beta-\sin\alpha\sin\beta$
然後分別維護區間$\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行有助於卡常呢