2023牛客暑期多校训练营10 C Multiplication（复杂度证明）

拿了个一血，队友说复杂度可能有问题，我就来证一证
做法就不说了，题解里有，贴个代码

def gcd(a,b):
    if (b==0):return a
    return gcd(b,a%b)
n=int(input())
k=int(input())
ans=0
for x in range(0,102):
    for y in range(0,102-x):
        xx=10**x
        yy=10**y
        a=yy-k
        b=k*xx-1
        if (a<=0 and b<=0):
            a=-a
            b=-b
        if (a<0 or b<0): continue
        g=gcd(a,b)
        if (g==0):continue
        a//=g
        b//=g
        t=1e100
        if (a>0):t=min(t,(yy-1)//a)
        if (b>0):t=min(t,(xx-1)//b)
        if (a*xx+b>0):t=min(t,n//(a*xx+b))
        if (t==1e100):continue
        a*=t
        b*=t
        if (a<(yy//10) or b<(xx/10)):continue
        tmp=a*xx+b
        ans=max(ans,tmp)
print(ans)

其实题解里的$O(L^{4}logL)$有点假，
枚举$L^2$，$gcd$理论次数$log(10^L)=LlogL$，还有除法复杂度设为$O(P)$
总复杂度是$O(L^{3}logL\cdot P)$，$P暴力是L^2$，应该可以做到$LlogL$，所以题解里说的复杂度应该是$O(L^{4}lo{g}^{2}L)$
如果用$stein$算法求$gcd$的话是比较严格的$O(L^{4}logL)$，因为只有除二和减法，都是$O(L)$的，但会跑满，肯定很慢

回到正文，先贴式子$(10^y-k,k\cdot 10^b-1)$
如果$min(x,y)<18$做一次高精度除法以后就相当于普通$gcd$了，复杂度$O(P+log(10^{18}))$
否则假设$x>y$，该式$=(10^y-k,k^2\cdot 10^{x-y}-1)$，如果把$x+y$作为整体的话，每次都会减至少$18$，否则会跳转上一种情况
总次数$200/18<12$，抛开枚举的复杂度$O(12P)$，实际上比如$100,100$这种情况是不可能出现的，我算他是$(100+100)/18$，实际只有一次
如果打个表，这个数不会超过$7$（我没打，毛估估）
$x<=y$其实也一样，因为哪怕$k^{12}$也是个无关紧要的数，远小于$10^{18}$
故总复杂度$O(12\cdot L^{2}P)$，其实还远跑不满，这只是我用纯数学能证到的地步

回头一看发现傻逼了，根本不需要这么麻烦的证明
不妨令$x>y$，若$x>2y$，取模后$x$规模一定减半
若$x<=2y$，根据$(10^y-k,k\cdot 10^b-1)=(10^y-k,k^2\cdot 10^{x-y}-1)$，$x-y$后规模一定减半
也不用按$18$来分块了，每次高精除的复杂度和长度是相关的（应该可以视作$P=LlogL$），
故总复杂度$O(L^2\cdot (L+L/2+L/4+...)logL)=O(L^{3}logL)$