Codeforces Round #356 (Div. 2)D. Bear and Tower of Cubes

链接:http://codeforces.com/contest/680/problem/D

题意:给定一个m,要求构造一个最优的m'<=m,m'每次分解一个最大的a满足a^3<=m'那么最后m'会被分解成一个唯一的等式m'=a1^3+a2^3+..+ak^3,最优的要求是k最大,k相同时要求a1+a2+..+ak最大。

分析:m是1e15级别。我们设f(m)表示m对应的答案且x为最大的正整数满足x*x*x<=m,我们判断x^3是全部被取完还是剩一个x^3来当x^3-1去凑最优的小的数,那么有f(m)=max(f(m%x^3),f(x^3-1)-1)+m/(x^3)。这样的话就能递归得到最优解了。直接递归的话时间会被爆,我们观察一下这个递归m%x^3<m/2那么f(m%x^3)只有logm层,x^3-1只有m^(1/3)个,那么我们对x^3-1记忆化一下就能处理了。O(m^(1/3)*logm)。

代码:

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=100010;
const int MAX=1000000100;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=998244353;
const int INF=1000000010;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
struct node {
    ll x,y;
}p[100010];
int q[100010];
node getans(ll n) {
    ll i;
    node a,b;
    a.x=a.y=0;
    if (n==0) return a;
    i=(ll)pow(n+1,1/3.0);
    if (n==i*i*i-1&&q[i]) return p[i];
    a=getans(n%(i*i*i));
    a.x+=n/(i*i*i);a.y+=(n/(i*i*i))*i*i*i;
    if (!q[i]) {
        b=getans(i*i*i-1);q[i]=1;p[i]=b;
    } else b=p[i];
    b.x+=n/(i*i*i)-1;b.y+=(n/(i*i*i)-1)*i*i*i;
    if (a.x!=b.x) {
        if (a.x>b.x) return a;
        else return b;
    } else if (a.y>b.y) return a;
        else return b;
}
int main()
{
    ll n;
    node ans;
    scanf("%I64d", &n);
    memset(q,0,sizeof(q));
    ans=getans(n);
    printf("%I64d %I64d\n", ans.x, ans.y);
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值