#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define X first #define Y second #define pb push_back #define rep(X,Y) for (int (X) = 0;(X) < (int)(Y);++(X)) #define reps(X,S,Y) for (int (X) = (int)(S);(X) < (int)(Y);++(X)) #define rrep(X,Y) for (int (X) = (int)(Y)-1;(X) >=0;--(X)) #define rreps(X,S,Y) for (int (X) = (int)(Y)-1;(X) >= (int)(S);--(X)) #define repe(X,Y) for ((X) = 0;(X) < (Y);++(X)) #define peat(X,Y) for (;(X) < (Y);++(X)) #define all(X) (X).begin(),(X).end() #define rall(X) (X).rbegin(),(X).rend() #define eb emplace_back #define UNIQUE(X) (X).erase(unique(all(X)),(X).end()) #define Endl endl #define NL <<"\n" #define cauto const auto using namespace std; using ll=long long; using pii=pair; using pll=pair; template using vv=vector>; template inline bool MX(T &l,const T &r){return l inline bool MN(T &l,const T &r){return l>r?l=r,1:0;} //#undef NUIP #ifdef NUIP #include "benri.h" #else #define out(args...) #endif #ifdef __cpp_init_captures templatevector table(int n, T v){ return vector(n, v);} template auto table(int n, Args... args){auto val = table(args...); return vector(n, move(val));} #endif template pair operator+(const pair &p,const pair &q){ return {p.X+q.X,p.Y+q.Y};} template pair& operator+=(pair &p,const pair &q){ p.X+=q.X; p.Y+=q.Y; return p;} template pair operator-(const pair &p,const pair &q){ return {p.X-q.X,p.Y-q.Y};} template pair& operator-=(pair &p,const pair &q){ p.X-=q.X; p.Y-=q.Y; return p;} template istream& operator>>(istream &is, pair &p){ is>>p.X>>p.Y; return is;} template T read(){ T re; cin>>re; return move(re);} template T read(const T &dec){ T re; cin>>re; return re-dec;} template vector readV(const int sz){ vector re(sz); for(auto &x:re) x=read(); return move(re);} template vector readV(const int sz, const T &dec){ vector re(sz); for(auto &x:re) x=read(dec); return move(re);} vv readG(const int &n,const int &m){ vv g(n); rep(_,m){ cauto a=read(1),b=read(1); g[a].pb(b); g[b].pb(a);} return move(g);} vv readG(const int &n){ return readG(n,n-1);} const ll MOD=1e9+7; //998244353 vector shu(vector a){ vector re; const int n=a.size()/2; rep(i,n){ re.eb(a[i]); re.eb(a[i+n]); } return re; } int bf(int n){ vector orig(2*n); iota(all(orig),0); int re=1; auto x=shu(orig); while(x!=orig){ ++re; x=shu(x); } return re; } ll extgcd(ll a,ll b,ll &x,ll &y){ ll g=a; x=1,y=0; if(b){ g=extgcd(b,a%b,y,x); y-=(a/b)*x; } return g; } ll inv(ll n,ll mod){ ll x,y; extgcd(n,mod,x,y); x%=mod; if(x<0) x+=mod; return x; } #include using namespace __gnu_pbds; const int SQ=33434; int solve(int n){ const ll m=2*n+1; vector mp; ll x=1; { rep(i,SQ){ if(i && x==1) return i; mp.eb(x,i); // mp[x]=i; x+=x; if(x>=m) x-=m; } } sort(all(mp)); mp.eb(m,0); int re=0; ll pw=x; reps(i,1,MOD){ cauto in=inv(pw,m); cauto it=lower_bound(all(mp),pii(in,0)); if(it->X==in) return i*SQ+it->Y; (pw*=x)%=m; } } int main(){ ios_base::sync_with_stdio(false); cin.tie(0); cout<