[BZOJ]3957: [WF2011]To Add or to Multiply 贪心 思维

Description

工业计算机处理器公司为顾客量身定做了非常快速、用于专门目的的处理单元。a-C-m系列的处理器（比如1-C-2和5-C-3）的指令集只有两种操作：
　　• A 数值加a
　　• M 数值乘m
　　处理器接收一个整数，执行一个A和M的指令序列（即程序）来修改输入，然后输出结果。举个例子，1-C-2处理器执行程序AAAM处理输入2返回输出10（计算过程是2→3→4→5→10）,然而5-C-3处理器用相同的程序和输入返回51（2→7→12→17→51）。
　　你是一个被指定做一个顶级秘密项目的a-C-m程序员。这意味着你不会被告知你的程序执行的精确输入。但会得到四个特别的值p，q，r，s，以及下列条件
　　1、输入保证是一个在p、q之间的数
　　2、输出必须是一个在r，s之间的数。
　　给你一个a-C-m处理器和p，q，r，s这四个数。你的工作是构想最短的a-C-m程序，使得任意任意x保证p≤x≤q，返回一个输出y使得r≤y≤s。如果有多个最短的程序，选择字典序最小的，而一个程序即是由A和M组成的字符串。

Solution

假设一开始的数为$x$，那么经过一堆操作后会变为：$m^{k}x+m^k{b}_{k}a+m^{k-1}{b}_{k-1}a+...m^0{b}_{0}a$，长度就是$k+{\sum }_{i=0}^k{b}_i$。特判掉$m=1$，$k$是$\log$级别的，直接枚举，然后要让$\sum b$尽量小，一定是从大到小，尽量用的最多，一个贪心就解决了。

Code

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const LL inf=4557430888798830399LL;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int a,m,p,q,r,s,Case=0;
int A[35],la,B[35],lb;
struct P
{
    LL x,l,c[35];
}ans;
LL Pow[35];
bool operator < (P a,P b)
{
    if(a.x<b.x)return true;
    if(a.x>b.x)return false;
    for(int i=1;i<=min(a.l,b.l);i++)
    {
        if(a.c[i]>b.c[i])return true;
        if(a.c[i]<b.c[i])return false;
    }return true;
}
vector<LL>Ans;
int main()
{
    while(1)
    {
        a=read(),m=read(),p=read(),q=read(),r=read(),s=read();
        if(a+m+p+q+r+s==0)break;
        printf("Case %d: ",++Case);
        if(q>s){puts("impossible");continue;}
        if(r<=p&&p<=s&&r<=q&&q<=s){puts("empty");continue;}
        if(m==1)
        {
            int L=(r-p+a-1)/a,R=(s-q)/a;
            if(L>R)puts("impossible");
            else printf("%dA\n",L);
            continue;
        }
        int w,tt=q;
        for(int i=1;;i++)
        {
            if(tt>s/m){w=i-1;break;}
            tt*=m;
        }
        Pow[0]=1;for(int i=1;i<=w;i++)Pow[i]=Pow[i-1]*m;
        ans.x=inf;
        for(int k=w;k>=0;k--)
        {
            int tp=Pow[k]*p,tq=Pow[k]*q;
            LL L=r-tp,R=s-tq,tmp;
            if(L>R)continue;
            P t;t.x=k;t.l=0;bool flag=false;
            for(int i=k;i>=0;i--)
            {
                if(L<=0&&0<=R)
                {
                    for(int j=i;j>=0;j--)t.c[++t.l]=0;
                    break;
                }
                LL ll=(L+Pow[i]-1)/Pow[i],rr=R/Pow[i];
                ll=(ll+a-1)/a,rr=rr/a;
                ll=min(ll,rr);
                if(ll*a*Pow[i]>=L){t.x+=ll,t.c[++t.l]=ll;L-=ll*a*Pow[i],R-=ll*a*Pow[i];continue;}
                if(ll==rr)t.x+=ll,t.c[++t.l]=ll,L-=ll*a*Pow[i],R-=ll*a*Pow[i];
                else t.x+=ll+1,t.c[++t.l]=ll+1,L-=(ll+1)*a*Pow[i],R-=(ll+1)*a*Pow[i];
            }
            if(L<=0&&0<=R)flag=true;
            if(!flag)continue;
            if(t<ans)ans=t;
        }
        if(ans.x==inf)puts("impossible");
        else
        {
            Ans.clear();
            for(int i=1;i<=ans.l;i++)
            {
                if(ans.c[i])Ans.push_back(ans.c[i]);
                if(i<ans.l)Ans.push_back(0);
            }
            int now=0;
            while(now<Ans.size())
            {
                LL cnt=0;
                while(!Ans[now]&&now<Ans.size())cnt++,now++;
                if(cnt)printf("%lldM ",cnt);
                else printf("%lldA ",Ans[now++]);
            }
            puts("");
        }
    }
}