Data-2021-01-03_zhr

[CERC2016]Easy Equation

构造一组初始解 \((k,1,0)\)，然后分析那个式子，得到了由一个解转移到另一个解的式子：\((a,b,c)\to(k\cdot (b+c),b,c)\)。

于是就可以做了。

luogu 上的原 SPJ 把我一开始的错误代码判过了，于是我发了个贴，现在管理员已经修掉了。

我当时还以为我写的高精复杂度怒踩 EI，实际上人用 Hash 判重不可能比我这个 log 要慢。

const int M=1e4;
int k,n;
int ps[4]={0,10,100,1000};
struct gyq
{
    int s[33];
    inline void mem0(){memset(s,0,sizeof(s));return;}
    inline void updata(){for(;s[s[0]+1];s[0]++,s[s[0]]%=M)s[s[0]+2]+=s[s[0]+1]/M;return;}
    inline void up(){for(int i=1;i<=s[0];i++)s[i+1]+=s[i]/M,s[i]%=M;updata();return;}
    inline void downdata(){for(;s[0]&&!s[s[0]];s[0]--);return;}
    inline void down(){for(int i=1;i<=s[0];i++)if(s[i]<0)s[i]+=M,s[i+1]--;downdata();return;}
    inline void pri(){if(!s[0])putchar(48);for(int i=s[0];i>=1;i--){if(i!=s[0])for(int j=0;j<4;j++)if(s[i]<ps[j])putchar(48);printf("%d",s[i]);}return;}
};
inline gyq operator +(gyq x,gyq y){gyq z=x;z.s[0]=max(z.s[0],y.s[0]);for(int i=1;i<=y.s[0];i++)z.s[i]+=y.s[i];z.up();return z;}
inline gyq operator -(gyq x,gyq y){gyq z=x;z.s[0]=max(z.s[0],y.s[0]);for(int i=1;i<=y.s[0];i++)z.s[i]-=y.s[i];z.down();return z;}
inline gyq operator *(gyq x,int y){gyq z=x;for(int i=1;i<=z.s[0];i++)z.s[i]*=y;z.up();return z;}
inline bool operator >(gyq x,gyq y){if(x.s[0]!=y.s[0])return x.s[0]>y.s[0];for(int i=x.s[0];i>=1;i--){if(x.s[i]>y.s[i])return true;if(x.s[i]<y.s[i])return false;}return false;}
inline bool operator >=(gyq x,gyq y){if(x.s[0]!=y.s[0])return x.s[0]>y.s[0];for(int i=x.s[0];i>=1;i--){if(x.s[i]>y.s[i])return true;if(x.s[i]<y.s[i])return false;}return true;}
inline bool operator <(gyq x,gyq y){if(x.s[0]!=y.s[0])return x.s[0]<y.s[0];for(int i=x.s[0];i>=1;i--){if(x.s[i]<y.s[i])return true;if(x.s[i]>y.s[i])return false;}return false;}
inline bool operator <=(gyq x,gyq y){if(x.s[0]!=y.s[0])return x.s[0]<y.s[0];for(int i=x.s[0];i>=1;i--){if(x.s[i]<y.s[i])return true;if(x.s[i]>y.s[i])return false;}return true;}
inline bool operator ==(gyq x,gyq y){for(int i=0;i<=x.s[0];i++)if(x.s[i]!=y.s[i])return false;return true;}
inline void jh(gyq x,gyq y){gyq z=x;x=y;y=z;return;}

map<gyq,bool>g;

struct zjj
{
    gyq a,b,c;
    inline void mem0(){a.mem0();b.mem0();c.mem0();return;}
    inline void init(){mem0();a.s[0]=b.s[0]=b.s[1]=1;a.s[1]=k;return;}
    inline void push(gyq x,gyq y,gyq z){a=x;b=y;c=z;if(a<b)jh(a,b);if(a<c)jh(a,c);if(b<c)jh(b,c);return;}
    inline void pri(){a.pri();putchar(' ');b.pri();putchar(' ');c.pri();putchar('\n');return;}
    inline bool check(){if(g[a]||g[b]||g[c]||a==b||b==c||a==c)return false;g[a]=true;g[b]=true;g[c]=true;return true;}
};
queue<zjj>d;
inline zjj start_zjj(){zjj st;st.init();return st;}
inline bool operator >(zjj x,zjj y){return (x.a==y.a)?((x.b==y.b)?(x.c>y.c):(x.b>y.b)):(x.a>y.a);}
inline bool operator <(zjj x,zjj y){return (x.a==y.a)?((x.b==y.b)?(x.c<y.c):(x.b<y.b)):(x.a<y.a);}

inline void add(gyq a,gyq b,gyq c)
{
    gyq z=(b+c)*k;if(a>=z)return;
    zjj ad;ad.push(z-a,b,c);
    d.push(ad);
    return;
}

int cutt;
inline void work()
{
    zjj now=d.front();d.pop();
    if(now.a.s[0]<=25&&now.c.s[0]&&now.check())now.pri(),cutt--;
    add(now.a,now.b,now.c);
    add(now.b,now.a,now.c);
    add(now.c,now.a,now.b);
    return;
}
int main()
{
    gyq Alpha,Beta;
    int i,j;
    k=rin();n=rin();
    d.push(start_zjj());
    for(cutt=n;cutt;)work();
    return 0;
}