bzoj1002 [FJOI2007]轮状病毒——找规律+高精度
2024-10-20 09:22:19
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1002
打表找规律,似乎是这样:https://blog.csdn.net/fzhvampire/article/details/46389897
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n;
struct N{
int a[];
N(){memset(a,,sizeof a);}
}f[],ans;
N add(N x,N y)
{
N ret;
ret.a[]=max(x.a[],y.a[]);
for(int i=;i<=ret.a[];i++)
ret.a[i]=x.a[i]+y.a[i];
for(int i=;i<=ret.a[];i++)
{
int t=ret.a[i]/; ret.a[i]%=;
ret.a[i+]+=t;
}
if(ret.a[ret.a[]+])ret.a[]++;
return ret;
}
N pw(N x)
{
N ret;
ret.a[]=*x.a[];
for(int i=;i<=x.a[];i++)
for(int j=;j<=x.a[];j++)
ret.a[i+j-]+=x.a[i]*x.a[j];
for(int i=;i<=ret.a[];i++)
{
int t=ret.a[i]/; ret.a[i]%=;
ret.a[i+]+=t;
}
if(ret.a[ret.a[]+])ret.a[]++;
while(ret.a[ret.a[]]==&&ret.a[])ret.a[]--;
return ret;
}
void sub(int k)
{
ans.a[]-=k; int i=;
while(ans.a[i]<)ans.a[i]+=,ans.a[++i]--;
while(ans.a[ans.a[]]==&&ans.a[])ans.a[]--;
}
void print()
{
for(int i=ans.a[];i;i--)printf("%d",ans.a[i]);
}
int main()
{
scanf("%d",&n);
f[].a[]=; f[].a[]=;
f[].a[]=; f[].a[]=;
for(int i=;i<=n;i++)
f[i]=add(f[i-],f[i-]);
ans=pw(f[n]);
if(n%==)sub();
print();
return ;
}
还有大家都在说的结论:http://vfleaking.blog.163.com/blog/static/17480763420119685112649/
也就是:f[i] = ( f[i-1]*3 - f[i-2] + 2 )
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n;
struct N{
int a[];
// N(){memset(a,0,sizeof a);}
}f[],ans;
N sub(N x,N y)
{
x.a[]+=;
int j=;
while(x.a[j]>=){x.a[j]%=;x.a[j+]++;j++;}
for(int i=;i<=x.a[];i++)
{
x.a[i]-=y.a[i];
if(x.a[i]<)x.a[i]+=,x.a[i+]--;
}
// for(int i=1;i<=x.a[0];i++)
// if(x.a[i]<0)x.a[i]+=10,x.a[i+1]--;
while(x.a[x.a[]]==&&x.a[])x.a[]--;
return x;
}
//N add(N x,int k)
//{
// x.a[1]+=k; int i=1;
// while(x.a[i]>=10)x.a[i+1]+=x.a[i]/10,x.a[i]%=10,i++;
// if(x.a[x.a[0]+1])x.a[0]++;
// return x;
//}
N mul(N x,int k)
{
for(int i=;i<=x.a[];i++)x.a[i]*=k;
for(int i=;i<=x.a[];i++)
{
x.a[i+]+=x.a[i]/; x.a[i]%=;
}
if(x.a[x.a[]+])x.a[]++;
return x;
}
void print()
{
for(int i=f[n].a[];i;i--)printf("%d",f[n].a[i]);
}
int main()
{
scanf("%d",&n);
f[].a[]=; f[].a[]=;
f[].a[]=; f[].a[]=;
for(int i=;i<=n;i++)
{
// f[i]=mul(f[i-1],3);
// f[i]=sub(f[i],f[i-2]);
// f[i]=add(f[i],2);
f[i]=sub(mul(f[i-],),f[i-]);
}
print();
return ;
}
果然考场上还是打表找规律比较靠谱...
然而为什么,我的代码跑得好慢...
最新文章
- 启动Tomcat内存溢出解决:java.lang.OutOfMemoryError: PermGen space
- 【leetcode】 Remove Duplicates from Sorted List
- Sprint.Net 笔记
- iOS - UIStoryboard
- linq里的select和selectmany操作
- LCA 笔记
- CSU1612Destroy Tunnels(强连通)
- linux C高手成长过程---书籍推荐
- boosting和bagging
- Javascript Design Patterns - Js Class
- WCF学习系列二_使用IIS发布WCF服务
- asp.net mvc5 设置Area下的为启动页
- spoj 104 Highways (最小生成树计数)
- 转每天一个linux命令(8):cp 命令
- 重温吕鑫MFC教学视频(一)
- bzoj 5297: [Cqoi2018]社交网络
- region server 中的OOM原因
- Python——Flask框架——程序的基本结构
- 我的Git
- 2019 CCPC wannfly winter camp Day 5