2015 ACM/ICPC 沈阳赛区网络赛 1010.Jesus Is Here (HDOJ5459)

把s[i]分为s[i-1]的内部组合,s[i-2]的内部组合,以及s[i-1]和s[i-2]的内部组合这三部分,接下来推出递推公式就可以了

#include<stdio.h>
#define N 250000
#define mod 530600414
long long a[N],b[N],c[N],d[N],s[N],k[N],v[N];
int main(){
	s[3]=0;
	s[4]=0;
	s[5]=5;	a[5]=5;	b[5]=0;	c[5]=2;	d[5]=5; k[5]=0,v[5]=5;
	s[6]=16;a[6]=11; b[6]=5;	c[6]=3;	d[6]=13; k[6]=2,v[6]=8;
	for(int i=7;i<=201314;i++){
		b[i]=a[i-1]*c[i-2]%mod+d[i-2]*c[i-1]%mod+(k[i-1]+3)*c[i-2]*c[i-1]%mod;
		s[i]=(s[i-1]+s[i-2]+b[i])%mod;

		a[i]=(a[i-2]+a[i-1]+v[i-2]*c[i-1]%mod+(k[i-1]+3)*c[i-1]%mod)%mod;
		d[i]=(d[i-1]+d[i-2]+v[i-1]*c[i-2]%mod+(k[i-1]+3)*c[i-2])%mod;
		v[i]=(v[i-1]+v[i-2]+k[i-1]+3)%mod;
		k[i]=k[i-2];
		c[i]=(c[i-1]+c[i-2])%mod;

		//printf("%I64d\n",i,s[i]);
	}
	int T;
	scanf("%d",&T);
	int n;
	int cas=1;
	while(T--){
        scanf("%d",&n);
		printf("Case #%d: %I64d\n",cas++,s[n]);
	}
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值