結果

問題 No.510 二次漸化式
ユーザー ixmelixmel
提出日時 2017-04-08 13:33:54
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,798 bytes
コンパイル時間 1,119 ms
コンパイル使用メモリ 98,420 KB
実行使用メモリ 87,288 KB
最終ジャッジ日時 2024-07-17 23:42:24
合計ジャッジ時間 15,508 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 AC 305 ms
13,568 KB
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 AC 326 ms
87,188 KB
testcase_17 AC 323 ms
87,236 KB
testcase_18 AC 324 ms
87,288 KB
testcase_19 AC 324 ms
87,216 KB
testcase_20 AC 323 ms
87,084 KB
testcase_21 AC 321 ms
87,004 KB
testcase_22 AC 327 ms
87,036 KB
testcase_23 AC 533 ms
87,044 KB
testcase_24 AC 532 ms
87,200 KB
testcase_25 AC 531 ms
87,120 KB
testcase_26 AC 536 ms
87,288 KB
testcase_27 AC 536 ms
87,092 KB
testcase_28 AC 546 ms
87,152 KB
testcase_29 AC 538 ms
87,228 KB
testcase_30 AC 544 ms
87,016 KB
testcase_31 AC 587 ms
87,168 KB
testcase_32 AC 593 ms
87,124 KB
testcase_33 AC 625 ms
87,128 KB
testcase_34 AC 412 ms
87,132 KB
testcase_35 AC 337 ms
87,028 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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]<<endl;
		}else assert(-1);
	}
}
0