51nod 1491-黄金系统（数论->递推）

1491 黄金系统

题目来源：  CodeForces

基准时间限制：1 秒 空间限制：131072 KB 分值: 40  难度：4级算法题

 收藏

 关注

q=5√+12， 在黄金系统下面 a0a1...an 等于  ∑ni=0ai∗qn−i  ，其中  ai  是0或者1。

现在给出两个黄金系统下面的数字，请比较他们的大小。

Input

单组测试数据。
第一行有一个字符串A。
第二行有一个字符串B。
按照a0到an的顺序输入。
他们都是非空串，可能有前导0，并且只有0和1组成，长度不超过100000。

Output

如果A>B，输出>；
如果A=B，输出=；
如果A<B，输出<；

Input示例

00100
11

Output示例

=

System Message  (题目提供者)

C++的运行时限为：1000 ms ，空间限制为：131072 KB  示例及语言说明请按这里

 允许其他 AC 的用户查看此代码，分享代码才能查看别人的代码并有机会获得勋章

（根号5+1）/2.....和斐波那契数列的通项公式好像啊，这里就不回顾了，直接递推，先求出q^1,发现和q^0没啥关系啊，再往后求，关系就出来了，将q^2拆开，很容易发现等于1+q,咦，这不是第0项和第一项吗？接着往后推发现就是这关系，于是就得出结论了，q^n=q^(n-1)+q^(n-2);(n>=2)。。。。。。。

呢就好办了，每一位要么是0，要么是1，当a串中的第i位为1，则令f[i]+1；当b串中的第i位为1，令f[i]-1;然后从后往前递推到f[1]和f[0]就好啦。。最后一定是f[1]*q+f[0]的形式，若大于0，则a>b,小于0，则a<b;否则a=b....这样问题还不算完全解决，因为len有100000,这么大，递推过去，有可能会爆long long，哇，这里怎么搞，想了很久，就卡这了，于是翻了翻题解，看了半天看懂了2333333，其实好像有个不等式直接能说明问题。dalao：http://blog.csdn.net/drin_e/article/details/52234636

#include<map>          
#include<stack>          
#include<queue>          
#include<vector>          
#include<math.h>          
#include<stdio.h>          
#include<iostream>          
#include<string.h>          
#include<stdlib.h>  
#include<algorithm> 
#include<functional>  
using namespace std;          
typedef long long  ll;         
#define inf  1000000000     
#define MOD 1000000007           
#define maxn  100005
#define lowbit(x) (x&-x)          
#define eps 1e-6
char s1[maxn],s2[maxn];
double f[maxn];
int main(void)
{
	int i,j,len1,len2,flag=0;
	long double x=0,y=0;
	scanf("%s%s",s1+1,s2+1);
	len1=strlen(s1+1);
	len2=strlen(s2+1);
	for(i=1;i<=len1;i++)
		f[len1-i]=s1[i]-'0';
	for(i=1;i<=len2;i++)
		f[len2-i]+='0'-s2[i];
	for(i=max(len1,len2);i>=3;i--)
	{
		if(f[i]>1) flag=1;
		if(f[i]<-1) flag=2;
		f[i-1]+=f[i];f[i-2]+=f[i];
	}
	f[1]+=f[2];f[0]+=f[2];
	double ans=(sqrt((double)5)+1)/2.0;
	if(f[1]*ans+f[0]>0) flag=1;
	else if(f[1]*ans+f[0]<0) flag=2;
	if(flag==1) printf(">\n");
	else if(flag==2) printf("<\n");
	else printf("=\n");
	return 0;
}