题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=6222
题目大意:
寻找大于n的最小数t
使得t-1,t,t+1这三个数构成的三角形为整数
n<10^30
思路:
看到n这么大,显然是瞎jiba找规律然后打表?
1.首先打表,打出所有的满足条件的t;
2.然后我们美滋滋的发现了规律;
3.按得出的递推式构造出10^30以内所有满足条件的数t;
4.对于每个n,在t的表里面顺序查找即可;
打表代码:
利用海伦公式
#include<cstdio>
#include<cmath>
#include<iostream>
#include<cstring>
using namespace std;
int i,t,s,pd,delta;
int main()
{
for (t=1; t<=100000; t++)
{
int delta=3*(t*t-4);
int pd=floor(sqrt(delta));
if (pd*pd!=delta) continue;
delta=pd;
if ((t*delta)%4!=0) continue;
s=(t/4)*delta;
if (s==0) continue;
cout<<t<<' '<<s<<endl;
}
}
打表后看前面几个数(t和面积)
4 6
14 72
52 1170
194 16128
724 226974
2702 3159000
10084 44031786
然后发现 f[i]=4*f[i-1]+f[i-2]
主代码:
#include<cstdio>
#include<string>
#include<iostream>
#include<cstring>
using namespace std;
int a[1000],b[1000],c[1000];
string st[1000];
char S[1000];
int T,i,la,lb,j,L;
string sta,stb;
int main()
{
scanf("%d\n",&T);
st[1]="4";
st[2]="14";
for (i=3; i<=100; i++)
{
sta=st[i-2]; stb=st[i-1];
la=sta.size(); lb=stb.size();
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof(c));
for (j=0; j<la; j++) a[la-j]=sta[j]-'0';
for (j=0; j<lb; j++) b[lb-j]=stb[j]-'0';
for (j=1; j<=lb; j++) b[j]*=4;
for (j=1; j<=lb; j++) b[j+1]+=b[j]/10,b[j]=b[j]%10;
while (b[lb+1]!=0) lb++;
L=lb;
for (j=1; j<=L; j++) c[j]=b[j]-a[j];
for (j=1; j<=L; j++) if (c[j]<0) { c[j]=c[j]+10; c[j+1]--; }
while (c[L]==0) L--;
string temp="";
for (j=L; j>=1; j--) temp+='0'+c[j];
st[i]=temp;
}
// 利用公式做高精度运算,然后制表
for (i=1; i<=T; i++)
{
gets(S);
int p=1;
string temp="";
for (j=0; j<strlen(S); j++) temp+=S[j];
bool have=0;
for (j=1; j<=100; j++) if (st[j].size()>temp.size() || (st[j].size()==temp.size() && st[j]>=temp)) { have=1; break; }
if (have==0) cout<<-1<<endl; else cout<<st[j]<<endl;
}
// 顺序查找即可
}