P7385 "EZEC - 6" 문제 풀이 점프
Code:
#include
#define ll unsigned long long
using namespace std;
const ll p=1000000007;
inline ll Read(){
ll dx=0,fh=1;
char c=getchar();
while(c>'9'||c<'0'){
if(c=='-') fh=-1;
c=getchar();
}
while(c>='0'&&c<='9'){
dx=dx*10+c-'0';
c=getchar();
}
return dx*fh;
}
ll quick_pow(ll x,ll y){
ll r=1,base=x,yy=y;
while(y){
if(y&1) r=(r*base)%p;
base=(base*base)%p;
y>>=1;
}
return r%p;
}
ll n,a,b,m,spe;
struct Matricx{
ll val[5][5];
};
Matricx Mul(Matricx x,Matricx y){
Matricx Ans;
for(int i=1;i<=4;++i)
for(int j=1;j<=4;++j){
Ans.val[i][j]=0;
for(int k=1;k<=4;++k) Ans.val[i][j]=(Ans.val[i][j]+((x.val[i][k]%p)*(y.val[k][j]%p))%p)%p;
Ans.val[i][j]%=p;
}
return Ans;
}
Matricx quickpow(Matricx x,ll y){
Matricx r,base=x;
for(int i=1;i<=4;++i){
for(int j=1;j<=4;++j) r.val[i][j]=0;
r.val[i][i]=1;
}
while(y){
if(y&1) r=Mul(r,base);
base=Mul(base,base);
y>>=1;
}
return r;
}
Matricx A,B;
int main(){
n=Read(),a=Read(),b=Read(),m=Read();
a=a*570000004%p,b=b*570000004%p;
for(int i=1;i<=m;++i){
ll spe_num=Read(),spe_score=Read();
spe=(spe+quick_pow(a+b,spe_num-1)*b%p*spe_score%p)%p;
}
A.val[1][1]=1,A.val[1][2]=1,A.val[1][3]=1;
A.val[2][1]=0,A.val[2][2]=a+b,A.val[2][3]=a;
A.val[3][1]=0,A.val[3][2]=0,A.val[3][3]=2*b;
B=quickpow(A,n-1);
ll s1=(a+b*2)%p,f2=( ((a+b)*a)%p + (a*b*2)%p )%p;
printf("%u
",(((B.val[1][1]*s1%p + B.val[1][2]*f2%p)%p + B.val[1][3]*((4*b*b)%p)%p)%p +spe%p)%p);
return 0;
}
코드 를 다 친 후에 각 식 에 모델 이 빠 졌 는 지 확인 해 야 한 다 는 것 을 기억 하 세 요. 본인 은 이 물건 때문에 아침 9 시 부터 오후 3 시 반 까지 알 아 냈 습 니 다!
이 내용에 흥미가 있습니까?
현재 기사가 여러분의 문제를 해결하지 못하는 경우 AI 엔진은 머신러닝 분석(스마트 모델이 방금 만들어져 부정확한 경우가 있을 수 있음)을 통해 가장 유사한 기사를 추천합니다:
[HDU 5608] functionProblem Description There is a function f(x),which is defined on the natural numbers set N,satisfies the following eqaut...
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
CC BY-SA 2.5, CC BY-SA 3.0 및 CC BY-SA 4.0에 따라 라이센스가 부여됩니다.