前言
并不是所有的高精度题目都可以套模板去解决,很多时候,在一道题目里高精度只是一部分(甚至很多时候只是一小部分)。
题目概述
AC代码
#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
#define maxn 100010
struct Bigint
{
int a[maxn],len;
Bigint(int x=0)
{
memset(a,0,sizeof(a));
for(len=1;x;++len)
a[len]=x%10,x/=10;
--len;
}
int &operator[](int i)
{
return a[i];
}
void flatten(int L)
{
len = L;
for (int i = 1; i <= L; ++i)
{
a[i + 1] += a[i] / 10; a[i] %= 10;
}
for (; !a[len];)
len--;
}
void print()
{
for (int i = max(len, 1); i >= 1; --i)
printf("%d", a[i]);
}
};
Bigint operator+(Bigint a,Bigint b)
{
Bigint c;
int l=max(a.len,b.len);
for(int i=1;i<=l;++i)
c[i]=a[i]+b[i];
c.flatten(l+1);
return c;
}
Bigint operator*(Bigint a,int b)
{
Bigint c;
for(int i=1;i<=a.len;++i)
c[i]=a[i]*b;
c.flatten(a.len+11);
return c;
}
vector<int> fac;
void div(int x)
{
int sum=0;
for(int i=2;sum+i<=x;++i)
sum+=i,fac.push_back(i);
int left=x-sum;
while(left)
{
for(int i=fac.size()-1;i>=0&&left;--i)
fac[i]++,--left;
}
}
int main()
{
int n;
cin>>n;
if(n<5)
{
cout<<"1"<<" "<<n-1<<endl;
cout<<n-1<<endl;
return 0;
}
div(n);
Bigint ans(1);
for(int i=0;i<fac.size();++i)
{
printf("%d ",fac[i]);
ans=ans*fac[i];
}
cout<<endl;
ans.print();
return 0;
}
分析思路
1.把那些该死的高精度先抛掉!先来看看怎么乘积才能最大?小学的时候我们做过一种题目,问相同的周长下长方形和正方形面积哪个大?答案是正方形。,上了初中之后学了二次函数和基本不等式,严格证明了两个数和一定时,当且仅当他们相等时乘积最大。
2.然而这道题目要求不能相等,那么就是他们越接近越好(差1最好),那么多个数字呢?也是越接近越好,而且越多越好。因为每个部分乘积都最大的时候整体乘积也能最大!想到什么,“贪心”二字。
3.那么落实到代码上就很容易,我们从2开始(从1开始只能帮倒忙,当然,3和4是特例,特殊处理就可以了。)比如9,我们把它分解成2+3+4,就很完美。不放心也可以任意举一些其他情况,都是大不过这种情况的。但是,如果不能刚好分解成只相差1的呢?应该从后面加起。比如10,多出来的那个1应该加在最后面,变成2+3+5,如果是11呢?那就是2+4+5。为什么可以这样呢?因为这样可以让保证数字还不相同的情况下让更多的数字增长。
4.然后就是高精度的问题了,这里只会涉及到高精度乘以低精度,所以使用Bigint模板是完全可行的。会贪心策略但是不熟悉这个高精度模板的可以移步下面这题,这道题目没有其他的,就是考察高精度,可以先好好琢磨这个。
(自定义大整数类Bigint)洛谷P1009 阶乘之和**
文末广告
学习算法和数据结构真的是个很累的过程,不会做只能求助于题解。 因为写代码这个东西基本上是千人千面。同时网络上搜到的题解很多要么用到的是自己还没学到的知识,看不懂;要么内核过于简陋,只能糊弄当前题目,不具有普适性。
如果你是一个喜欢做洛谷,ACwing和PTA的题目的同学,欢迎关注我的博客,我主要在这三个平台上做题,认为有价值和有难度的题目我会写题解发布出来。
TreeTraverler的往期文章