数和分解,乘积最大-vijos整数分解(版本2)

https://vijos.org/p/1033
把一个数和分解,让其乘积最大;
我们尽量分解3,然后分解2;
换句话说,分解好了,只有3,2,,且3最多;
为什么我们分解3,2;
为什么我们不分解5呢?
因为5可以分解成3,2
3*2=6>5;
另外的数同理;
为什么我们让3最多呢?
因为3分解2,1
2*1=2<3;
2*2*2<3*3;
所以我口胡好了;
然后直接高精模版(然后我发现我的ksm模版有bug)

#include<iostream>
#include<cstdio>
#include<cstring>
#define Ll long long
using namespace std;
const Ll NN=1e8;
const int N=8;
struct H{
    Ll a[1250],len;
    H(){memset(a,0,sizeof a);len=1;}
    void clear(){memset(a,0,sizeof a);len=1;}
};
int n,x,y;
void outit(H a){
    printf("%lld",a.a[a.len]);
    for(int i=a.len-1;i;i--){
        for(int k=NN/10;a.a[i]<k;k/=10)printf("0");
        if(a.a[i])printf("%lld",a.a[i]); 
    }printf("\n");  
}
H in(int x){
    H a;
    if(!x)return a;
    a.len=0;
    while(x){a.a[++a.len]=x%NN;x=x/NN;}
    return a;
}
H chen(H a,H b){
    H z;
    z.len=a.len+b.len+2;
    for(int i=1;i<=a.len;i++)
    for(int j=1;j<=b.len;j++)z.a[i+j-1]+=(a.a[i]*b.a[j]);
    for(int i=1;i<=z.len;i++)z.a[i+1]+=z.a[i]/NN,z.a[i]%=NN;    
    while(z.len>1&&!z.a[z.len])z.len--;
    return z;
}
H ksm(H a,int n){
    if(n==0)return in(1);
    if(n==1)return a;
    H c=ksm(a,n>>1);
    c=chen(c,c);
    if(n&1)c=chen(c,a);
    return c;
}
int main()
{
    scanf("%d",&n);
    if(n==1){
        cout<<1;return 0;
    }
    if(n==0){
        cout<<0;return 0;
    }
    x=n/3;
    y=n%3;
    if(y==1){
        y=2;x-=1;
    }else
    if(y==2)y=1;
    H ans;
    H a;
    a=in(3);
    ans=ksm(a,x);
    a=in(2);
    ans=chen(ans,ksm(a,y));
    outit(ans);
}
posted @ 2017-03-20 10:38  largecube233  阅读(172)  评论(0编辑  收藏  举报