結果
問題 |
No.1099 Range Square Sum
|
ユーザー |
|
提出日時 | 2020-06-27 01:04:30 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,191 ms / 2,000 ms |
コード長 | 2,368 bytes |
コンパイル時間 | 1,087 ms |
コンパイル使用メモリ | 102,600 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-05 03:31:22 |
合計ジャッジ時間 | 11,187 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 30 |
ソースコード
#include<iostream> #include<vector> #include<string> #include<map> #include<set> #include<queue> #include<algorithm> #include<cmath> #include<iomanip> #include<random> #include<stdio.h> using namespace std; typedef long long ll; long long Mod(long long val, long long m) { long long res = val % m; if (res < 0) res += m; return res; } vector<ll> A(200000); ll Adata[200000]; ll RSQdata[200000]; ll RSSQdata[200000]; ll n; ll RSQ(ll l,ll r){ ll id=l; ll L,q,R; ll res=0; while (r>=id){ L=id; q=id/n; R=min((q+1)*n-1,r); if (R-L+1==n){ res+=RSQdata[q]; } else{ for (ll i=L;i<R+1;i++){ res+=A[i]+Adata[q]; } } id=R+1; } return res; } void update(ll l,ll r,ll x){ ll id=l; ll L,q,R; while (r>=id){ L=id; q=id/n; R=min((q+1)*n-1,r); if (R-L+1==n){ Adata[q]+=x; RSQdata[q]+=x*n; RSSQdata[q]+=2*x*RSQ(L,R)-n*x*x; } else{ for (ll i=L;i<R+1;i++){ A[i]+=x; } RSQdata[q]+=(R-L+1)*x; RSSQdata[q]+=2*x*RSQ(L,R)-(R-L+1)*x*x; } id=R+1; } } ll RSSQ(ll l,ll r){ ll id=l; ll L,q,R; ll res=0; while (r>=id){ L=id; q=id/n; R=min((q+1)*n-1,r); if (R-L+1==n){ res+=RSSQdata[q]; } else{ for (ll i=L;i<R+1;i++){ res+=(A[i]+Adata[q])*(A[i]+Adata[q]); } } id=R+1; } return res; } int main(){ ll N; scanf("%lld",&N); for (ll i=0;i<N;i++){ scanf("%lld",&A[i]); } while (N>=n*n){ n+=1; } n-=1; for (ll i=0;i<N;i++){ ll q=i/n; RSQdata[q]+=A[i]; RSSQdata[q]+=A[i]*A[i]; } ll Q; scanf("%lld",&Q); while (Q--){ int q; scanf("%lld",&q); if (q==1){ ll l,r,x; scanf("%lld",&l); scanf("%lld",&r); scanf("%lld",&x); l-=1;r-=1; update(l,r,x); } else{ ll l,r; scanf("%lld",&l); scanf("%lld",&r); l-=1;r-=1; ll ans=RSSQ(l,r); cout<<ans<<endl; } } }