码迷,mamicode.com
首页 > 其他好文 > 详细

对你的爱深不见底

时间:2017-08-05 19:51:51      阅读:249      评论:0      收藏:0      [点我收藏+]

标签:gui   roo   uil   ovs   mos   pfx   a10   hnoi   oci   

对你的爱深不见底

技术分享

技术分享

这题的话,其实可以先打个表找一下规律,然后会发现,n是没什么卵用的,你只要关注m就行了.然而,m很大很大很大,不得不开高精,这就为解题带来了麻烦.

我们先把答案的序列打出来:

技术分享

乍一看,很有规律的样子(别说你没看出来).我们再细分一下:

技术分享

我们发现,划分成的子区间都是公差为1的等差数列,且区间长都是fib数.我们在观察一下:

技术分享

我们发现,每个区间的第一行元素+1也是fib数!

那么,如果我们知道一个数m,我们要知道ans,我们可以通过两个fib数列推出!

我们先可以求出m>sum(fib[1~i])(fib[1]=1,fib[2]=2)中i的最大值.也就是说,我们已经推到了数m对应ans所在的区间的上一个区间的末行.

而且,我们又可以在同时求出ans所在区间的第一行的元素,那么我们就可以轻松知道ans了.

我们设数m对应ans所在的区间的上一个区间的末行的元素为x,ans所在区间的第一行的元素为y

则ans=y+m-x-1.(可以自己去推一推)

最后,注意输出的数要对那个万恶的数(为什么万恶?额,没什么)取模再输出.

技术分享
 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #define LL long long
 5 using namespace std;
 6 const int maxl=1005;
 7 const int TT=258280327;
 8 int n,ans;
 9 struct bigint{
10     int len,a[maxl];
11     bigint(){len=0,memset(a,0,sizeof a);}
12     void read(){
13         char c[maxl]; scanf("%s",c);
14         len=strlen(c); for (int i=1; i<=len; i++) a[i]=c[i-1]-0;
15         for (int i=1; i<=len/2; i++) swap(a[i],a[len-i+1]);
16     }
17     bigint operator + (const bigint b){
18         bigint c; c.len=max(len,b.len);
19         for (int i=1; i<=c.len; i++){
20             c.a[i]+=a[i]+b.a[i];
21             c.a[i+1]=c.a[i]/10;
22             c.a[i]%=10;
23         }
24         while (c.a[c.len+1]) c.len++;
25         return c;
26     }
27     bigint operator - (const bigint b){
28         bigint c; c.len=len;
29         for (int i=1; i<=c.len; i++){
30             c.a[i]+=a[i]-b.a[i];
31             if (c.a[i]<0) c.a[i]+=10,c.a[i+1]--;
32             c.a[i+1]+=c.a[i]/10;
33             c.a[i]%=10;
34         }
35         while (c.len>1&&!c.a[c.len]) c.len--;
36         return c;
37     }
38     bool operator > (const bigint b){
39         if (len>b.len) return 1; else if (len<b.len) return 0;
40         for (int i=len; i; i--) if (a[i]>b.a[i]) return 1;
41         else if (a[i]<b.a[i]) return 0;
42         return 0;
43     }
44     void write(){
45         for (int i=len; i>=1; i--) a[i-1]=(a[i-1]+(LL)a[i]*10%TT)%TT;
46         printf("%d",a[1]); putchar(\n);
47     }
48 }m,x,y,z,a,b,c,final,dif,one,tot;
49 int main(){
50     int T; cin>>T; one.len=1,one.a[1]=1;
51     for (; T; T--){
52         cin>>n,ans=0,m.read();
53         if (m.len==1&&m.a[1]<7){
54             if (m.a[1]==1) puts("0"); else
55             if (m.a[1]==2) puts("0"); else
56             if (m.a[1]==3) puts("1"); else
57             if (m.a[1]==4) puts("1"); else
58             if (m.a[1]==5) puts("2"); else
59             if (m.a[1]==6) puts("3");
60             continue;
61         }
62         m=m-one;
63         x.len=1,x.a[1]=2,y.len=1,y.a[1]=3,z=x+y;
64         a.len=1,a.a[1]=1,b.len=1,b.a[1]=2,c=a+b;
65         tot=z;
66         while (m>tot){
67             c=a+b; a=b; b=c;
68             final=tot+z;
69             if (!(m>final)){final=tot; break;} else tot=final;
70             x=y; y=z; z=x+y;
71         }
72         dif=m-final;
73         c=c-one; final=c+dif; final=final-one;
74         final.write();
75     }
76     return 0;
77 }
View Code

 

对你的爱深不见底

标签:gui   roo   uil   ovs   mos   pfx   a10   hnoi   oci   

原文地址:http://www.cnblogs.com/whc200305/p/7291106.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!