Problem Description
作为CNCS的半壁江山,狗哥常常在宇宙中心邵阳眺望黄浦江,夜晚的星空总是迷人,有时候还能见到彗星滑落。
狗哥是幸运的,他在两秒钟内看到了十七颗彗星划过天际,作为打ACM的学者,自然不会有「稳定-1」情况。他开始研究彗星运动的轨迹,发现他们都遵照斐波那契螺旋线在运动着。
尤里卡!狗哥觉得这就是找寻「生命,宇宙和一切的终极答案」的精要所在,但是怎么表示呢?狗哥觉得求取斐波那契螺旋线经过的一个个方格的面积之和就是公式的表现。
例如下图,螺旋线每划过一个方格,都转过了四分之一圈。如果我们以四分之一圈为单位,那么我们用类似带分数的形式表示螺旋线转动的起点和终点。例如,0+0 到 0 + 1 意即从第一个方格转到第二个方格,划过了前两个方格,他们的面积之和为2(1+1)。同理,0+0 到 1+0 划过了前五个方格,他们的面积之和为40(1+1+4+9+25)。
但是聪明的狗哥需要一个程序去获得指定范围内的螺旋线面积之和,狗哥给了你一首「希望之花」的时间,而他需要利用这个时间去打出四暗刻单骑。如果你能完成这个程序,狗哥会封你为格拉摩根伯爵
Input
不定组数据。
首先输入一个整数Q,代表狗哥询问次数。
接下来Q行,每行四个整数a,b,c,d,代表狗哥想求 a+b 到 c+d 之间的螺旋线面积之和。
1<= Q <= 10000
0<= a,c <= 10000
0 <= b,d <= 3
结果对192600817取模。
Output
一个数字,表示螺旋线面积之和。
Sample Input
4
0 0 0 1
0 0 1 0
1 2 2 1
1 1 0 3
4
0 0 0 1
0 0 1 0
1 2 2 1
1 1 0 3
Sample Output
2
40
4791
98
2
40
4791
98
#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int N=4e4+100;
ll m[N];
ll f[N];
int main()
{
m[1]=1;
f[1]=1;
for(ll i=2;i<=4e4+50;i++){
f[i]=(f[i-1]+f[i-2])%192600817;
m[i]=(f[i]*f[i])%192600817;
//cout<<m[i]<<endl;
m[i]=m[i]+m[i-1];
//cout<<m[i]<<endl;
}
int n,a,b,c,d;
while(~scanf("%d",&n)){
while(n--){
scanf("%d%d%d%d",&a,&b,&c,&d);
a=a*4+b+1;c=c*4+d+1;
if(a>c) swap(a,c);
printf("%I64d\n",(m[c]-m[a-1])%192600817);
}
}
//cout << "Hello world!" << endl;
return 0;
}
原文:https://www.cnblogs.com/-yjun/p/10544063.html