[cf1427E]Xum

假设$x$的最高位为$2^{t}$（即$2^{t}\le x<2^{t+1}$），并构造出$y=2^{t}x\oplus x$，不难发现两者仅在第$t$位上均为1，那么根据异或的性质可得$y=(2^{t}+1)x-2^{t+1}$

由于$x$为奇数，即$(x,2^{t+1})=1$，进而也即$(x,y)=1$

通过扩欧求出一组$ax+by=1$的解，并调整使得$0<a\le 2y$且$a\equiv 1(mod\ 2)$，对应的$-2x\le b<0$，根据奇偶性可得$ax$为奇数（且$by$为偶数），那么$ax+by=1$即等价于$ax\oplus (-b)y=1$，由此计算可得

另外，计算过程中需要实现乘法，这借助类似快速幂的做法实现即可

最终操作次数（和时间复杂度）约为$o(\log n)$，数字范围约为$2n^{3}$，可以通过

#include<bits/stdc++.h>
using namespace std;
#define ll long long
struct Data{
    int p;
    ll x,y;
};
vector<Data>ans;
ll n,m,x,y;
void calc(ll n,ll m){
    m--;
    ll s=n,sum=n;
    while (m){
        if (m&1){
            ans.push_back(Data{1,s,sum});
            sum+=s;
        }
        ans.push_back(Data{1,s,s});
        s<<=1,m>>=1;
    }
}
void exgcd(ll a,ll b,ll &x,ll &y){
    if (!b){
        x=1,y=0;
        return;
    }
    exgcd(b,a%b,y,x);
    y-=(a/b)*x;
}
int main(){
    scanf("%lld",&n);
    ll t=2;
    while ((t<<1)<=n)t<<=1;
    calc(n,t);
    m=((n*t)^n),ans.push_back(Data{0,n*t,n});
    exgcd(n,m,x,y);
    x=(x%m+m)%m;
    if (x%2==0)x+=m;
    y=(n*x-1)/m;
    calc(n,x),calc(m,y);
    ans.push_back(Data{0,n*x,m*y});
    printf("%d\n",(int)ans.size());
    for(int i=0;i<ans.size();i++)
        if (ans[i].p)printf("%lld + %lld\n",ans[i].x,ans[i].y);
        else printf("%lld ^ %lld\n",ans[i].x,ans[i].y);
    return 0;
}