結果
問題 |
No.510 二次漸化式
|
ユーザー |
![]() |
提出日時 | 2017-04-08 13:36:52 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 509 ms / 3,000 ms |
コード長 | 2,804 bytes |
コンパイル時間 | 885 ms |
コンパイル使用メモリ | 98,212 KB |
実行使用メモリ | 87,348 KB |
最終ジャッジ日時 | 2024-07-17 23:44:45 |
合計ジャッジ時間 | 11,943 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 34 |
ソースコード
#include<iostream> #include<vector> #include<string> #include<algorithm> #include<map> #include<set> #include<utility> #include<cmath> #include<cstring> #include<queue> #include<stack> #include<cstdio> #include<sstream> #include<iomanip> #include<assert.h> #define loop(i,a,b) for(int i=a;i<b;i++) #define rep(i,a) loop(i,0,a) #define pb push_back #define all(in) in.begin(),in.end() #define shosu(x) fixed<<setprecision(x) using namespace std; //kaewasuretyuui typedef long long ll; typedef pair<int,int> pii; typedef vector<ll> vi; typedef vector<vi> vvi; typedef vector<pii> vp; typedef vector<vp> vvp; typedef vector<string> vs; typedef vector<double> vd; typedef vector<vd> vvd; typedef pair<int,pii> pip; typedef vector<pip>vip; const double PI=acos(-1); const double EPS=1e-7; const int inf=1e8; const ll INF=1e16; int dx[4]={0,1,0,-1}; int dy[4]={1,0,-1,0}; #define MOD 1000000007 typedef vvi SegT; vvi defvalue; vvi defvalue2; SegT f(SegT A,SegT B){ vvi C(A.size(),vi(B[0].size())); rep(i,A.size())rep(k,B.size())rep(j,B[0].size()) (C[i][j]+=A[i][k]*B[k][j])%=MOD; return C; } class SegTree{ private: vector<SegT>val; int n; SegT combine(SegT a,SegT b){return f(b,a);} public: SegTree(int size){ n=1; while(n<size)n<<=1; val=vector<SegT>(2*n,defvalue); } SegTree(const vector<SegT> &in){ n=1; while(n<in.size())n<<=1; val=vector<SegT>(2*n,defvalue); for(int i=n-1+in.size()-1;i>=0;i--){ if(n-1<=i)val[i]=in[i-(n-1)]; else val[i]=combine(val[i*2+1],val[i*2+2]); } } void update(int q,int i,ll a){ i+=n-1; if(!q)val[i][0][2]=a; else{ val[i][1][1]=a; val[i][2][1]=2*a; val[i][2][2]=a*a%MOD; } while(i>0){ i=(i-1)/2; val[i]=combine(val[i*2+1],val[i*2+2]); } } SegT query(int a,int b,int k=0,int l=0,int r=-1){//[a,b) if(r==-1)r=n; if(r<=a||b<=l)return defvalue2; if(a<=l&&r<=b)return val[k]; else return combine(query(a,b,k*2+1,l,(l+r)/2),query(a,b,k*2+2,(l+r)/2,r)); } void tmp(){ rep(i,val.size()){ cout<<i<<endl; rep(j,4){rep(k,4)cout<<" "<<val[i][j][k];cout<<endl;} cout<<endl; } } }; int main(){ vvi r=vvi(4,vi(4)); r[0][0]=r[1][3]=r[2][3]=r[3][3]=1; defvalue=r; r=vvi(4,vi(4)); r[0][0]=r[1][1]=r[2][2]=r[3][3]=1; defvalue2=r; int n,q; cin>>n>>q; SegTree st(n+1); assert(1<=n&&n<=100000); assert(1<=q&&q<=20000); int x,y; while(q--){ // st.tmp(); int a,b; char c; cin>>c; if(c=='x'){ cin>>a>>b; assert(0<=a&&a<n); assert(0<=b&&b<=100000); st.update(0,a,b); }else if(c=='y'){ cin>>a>>b; assert(0<=a&&a<n); assert(0<=b&&b<=100000); st.update(1,a,b); }else if(c=='a'){ cin>>a; assert(0<=a&&a<=n); vvi t=st.query(0,a); cout<<(t[0][0]+t[0][1]+t[0][2]+t[0][3])%MOD<<endl; }else assert(-1); } }