結果

問題 No.1409 Simple Math in yukicoder
ユーザー mjtaimjtai
提出日時 2021-02-27 00:07:13
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 29 ms / 2,000 ms
コード長 8,613 bytes
コンパイル時間 2,898 ms
コンパイル使用メモリ 226,948 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-02 16:35:52
合計ジャッジ時間 27,903 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 4 ms
5,248 KB
testcase_08 AC 12 ms
5,248 KB
testcase_09 AC 4 ms
5,248 KB
testcase_10 AC 10 ms
5,248 KB
testcase_11 AC 11 ms
5,248 KB
testcase_12 AC 9 ms
5,248 KB
testcase_13 AC 5 ms
5,248 KB
testcase_14 AC 4 ms
5,248 KB
testcase_15 AC 5 ms
5,248 KB
testcase_16 AC 5 ms
5,248 KB
testcase_17 AC 14 ms
5,248 KB
testcase_18 AC 8 ms
5,248 KB
testcase_19 AC 9 ms
5,248 KB
testcase_20 AC 7 ms
5,248 KB
testcase_21 AC 6 ms
5,248 KB
testcase_22 AC 11 ms
5,248 KB
testcase_23 AC 4 ms
5,248 KB
testcase_24 AC 3 ms
5,248 KB
testcase_25 AC 10 ms
5,248 KB
testcase_26 AC 11 ms
5,248 KB
testcase_27 AC 16 ms
5,248 KB
testcase_28 AC 17 ms
5,248 KB
testcase_29 AC 17 ms
5,248 KB
testcase_30 AC 17 ms
5,248 KB
testcase_31 AC 16 ms
5,248 KB
testcase_32 AC 17 ms
5,248 KB
testcase_33 AC 17 ms
5,248 KB
testcase_34 AC 16 ms
5,248 KB
testcase_35 AC 16 ms
5,248 KB
testcase_36 AC 17 ms
5,248 KB
testcase_37 AC 27 ms
5,248 KB
testcase_38 AC 28 ms
5,248 KB
testcase_39 AC 29 ms
5,248 KB
testcase_40 AC 27 ms
5,248 KB
testcase_41 AC 28 ms
5,248 KB
testcase_42 AC 28 ms
5,248 KB
testcase_43 AC 28 ms
5,248 KB
testcase_44 AC 27 ms
5,248 KB
testcase_45 AC 28 ms
5,248 KB
testcase_46 AC 27 ms
5,248 KB
testcase_47 AC 27 ms
5,248 KB
testcase_48 AC 27 ms
5,248 KB
testcase_49 AC 28 ms
5,248 KB
testcase_50 AC 27 ms
5,248 KB
testcase_51 AC 27 ms
5,248 KB
testcase_52 AC 27 ms
5,248 KB
testcase_53 AC 27 ms
5,248 KB
testcase_54 AC 27 ms
5,248 KB
testcase_55 AC 28 ms
5,248 KB
testcase_56 AC 28 ms
5,248 KB
testcase_57 AC 28 ms
5,248 KB
testcase_58 AC 28 ms
5,248 KB
testcase_59 AC 28 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

///////////////////////////////////////////////////////////////////////////////

#include <bits/stdc++.h>
#include <unistd.h>
#include <time.h>
#include <sys/time.h>

using namespace std;

///////////////////////////////////////////////////////////////////////////////

#define DEBUG 0

#define pb push_back
#define V vector
#define S static

#define rep(i,n) for(ll i=0LL; i<n; ++i)
#define srep(i,s,n) for(ll i=s; i<n; ++i)
#define rrep(i,n) for(ll i=n-1LL; i>=0LL; --i)
#define ALL(a) (a).begin(),(a).end()
#define RALL(a) (a).rbegin(),(a).rend()
#define CIN(x) do { \
        assert(!cin.eof()); \
        cin >> x; \
        assert(!cin.fail()); \
} while(0);
#if DEBUG
#define debug_print(...) _debug_print(__VA_ARGS__)
#define debug_printf(...) printf(__VA_ARGS__)
#define debug_print_time _debug_print_time
#else // DEBUG
#define debug_print(...)
#define debug_printf(...)
#define debug_print_time
#endif // DEBUG

typedef long long ll;
typedef unsigned long long ull;
typedef __int128_t ll128;
typedef tuple<ll, ll> t2;
typedef tuple<ll, ll, ll> t3;
typedef tuple<ll, ll, ll, ll> t4;
typedef tuple<ll, ll, ll, ll, ll> t5;
template <typename T>
using priority_queue_incr = priority_queue<T, V<T>, greater<T> >;

///////////////////////////////////////////////////////////////////////////////

template <typename TT, typename T>
void incr_m(map<TT, T> &m, TT k)
{
    if (m.find(k) == m.end()) m[k] = 0;
    m[k]++;
}

void llin(ll &a)
{
    CIN(a);
}

void llinl1(V<ll> &v, ll count)
{
    for (ll i = 0LL; i < count; ++i) {
        ll a;
        CIN(a);
        v.push_back(a);
    }
}

void llinl2(V<t2> &v, ll count)
{
    for (ll i = 0LL; i < count; ++i) {
        ll a, b;
        CIN(a >> b);
        v.push_back(t2(a, b));
    }
}

void llinl3(V<t3> &v, ll count)
{
    for (ll i = 0LL; i < count; ++i) {
        ll a, b, c;
        CIN(a >> b >> c);
        v.push_back(t3(a, b, c));
    }
}

void llinl4(V<t4> &v, ll count)
{
    for (ll i = 0LL; i < count; ++i) {
        ll a, b, c, d;
        CIN(a >> b >> c >> d);
        v.push_back(t4(a, b, c, d));
    }
}

void llina(V<ll> &v, ll count)
{
    llinl1(v, count);
}

template <typename T>
void sort(V<T> &v)
{
    sort(v.begin(), v.end());
}

template <typename T>
void sort_reverse(V<T> &v)
{
    sort(v.begin(), v.end(), greater<T>());
}

t2 _ext_gcd(ll a, ll b, ll g)
{
    if (!b) return t2(1, 0);
    ll q = a / b;
    ll r = a % b;
    t2 sol = _ext_gcd(b, r, g);
    ll sx = get<0>(sol);
    ll sy = get<1>(sol);
    ll x = sy;
    ll y = sx - q * sy;
    return t2(x, y);
}

t2 ext_gcd(ll a, ll b)
{
    return _ext_gcd(a, b, gcd(a, b));
}

// x and mod must be coprime
ll mod_inv(ll x, ll mod)
{
    t2 result = ext_gcd(x, mod);
    ll ret = get<0>(result);
    while (ret < 0) ret += mod;
    ret %= mod;
    return ret;
}

ll msec()
{
    struct timeval tv;
    gettimeofday(&tv, NULL);
    ll ret = 0;
    ret += tv.tv_sec * 1000LL;
    ret += tv.tv_usec / 1000LL;
    return ret;
}


template <typename T>
void _debug_print(T x)
{
    cout << x << " ";
}

void _debug_print(const t2 &x)
{
    ll x1 = get<0>(x);
    ll x2 = get<1>(x);
    cout << "(" << x1 << ", " << x2 << ") ";
}

void _debug_print(const t3 &x)
{
    ll x1 = get<0>(x);
    ll x2 = get<1>(x);
    ll x3 = get<2>(x);
    cout << "(" << x1 << ", " << x2 << ", " << x3 << ") ";
}

void _debug_print(const t4 &x)
{
    ll x1 = get<0>(x);
    ll x2 = get<1>(x);
    ll x3 = get<2>(x);
    ll x4 = get<3>(x);
    cout << "(" << x1 << ", " << x2 << ", " << x3 << ", " << x4 << ") ";
}

template <typename T>
void _debug_print(T xarray[], ll n)
{
    rep (i, n) _debug_print(xarray[i]);
    cout << endl;
}

template <typename T>
void _debug_print(const V<T> &xlist)
{
    for (auto x : xlist) _debug_print(x);
    cout << endl;
}

template <typename T>
void _debug_print(const set<T> &xset)
{
    for (auto x : xset) _debug_print(x);
    cout << endl;
}

template <typename TT, typename T>
void _debug_print(const map<TT, T> &xlist)
{
    for (auto x : xlist) {
        TT k = x.first;
        T v = x.second;
        cout << "K=";
        _debug_print(k);
        cout << " V=";
        _debug_print(v);
        cout << endl;
    }
}



ll mod_mlt(ll x, ll y, ll mod)
{
    ll ret = 0LL;
    x %= mod;

    while (y) {
        if (y & 1LL) {
            ret += x;
            ret %= mod;
        }
        y >>= 1;
        x <<= 1;
        x %= mod;
    }

    return ret;
}

// returns t2(solution, mod)
t2 chinese_remainder(ll a1, ll m1, ll a2, ll m2)
{
    assert(a1 >= 0);
    assert(m1 > 0);
    assert(a2 >= 0);
    assert(m2 > 0);

    ll mgcd = gcd(m1, m2);
    if (a1 % mgcd != a2 % mgcd) return t2(0, 0);
    ll mlcm = m1 * m2 / mgcd;
    t2 z = ext_gcd(m1, m2);
    ll z1 = get<0>(z);
    ll z2 = get<1>(z);

    // ll x = a1 + ((a2 - a1) / mgcd) * m1 * z1;

    ll x = z1;
    while (x < 0) x += mlcm;

    x = mod_mlt(x, m1, mlcm);

    ll coef = (a2 - a1) / mgcd;
    while (coef < 0) coef += mlcm;
    x = mod_mlt(x, coef, mlcm);

    x += a1;
    x %= mlcm;

    return t2(x, mlcm);
}

void get_divisors(V<ll> &retlist, ll x)
{
    for (ll i = 1LL; i < sqrt(x) + 3LL; ++i) {
        if (x % i == 0LL) {
            retlist.push_back(i);
            retlist.push_back(x / i);
        }
    }
}

// returns factors and 1
void get_factors(V<ll> &retlist, ll x)
{
    retlist.pb(1LL);
    for (ll i = 2LL; i < (ll)(sqrt(x)) + 3LL; ++i) {
        while (x % i == 0LL) {
            retlist.pb(i);
            x /= i;
        }
    }
    retlist.pb(x);
}

bool is_prime(ll x)
{
    V<ll> factors, factors2;

    get_factors(factors, x);

    for (auto factor : factors) {
        if (factor > 1) factors2.pb(factor);
    }

    return factors2.size() == 1 && x == factors2[0];
}

void eratosthenes(set<ll> &primes, ll n)
{
    bool *is_not_prime = new bool[n+3LL];
    memset(is_not_prime, 0, sizeof(bool) * (n+3LL));

    srep (v, 2LL, (ll)sqrt(n)+10LL) {
        if (is_not_prime[v]) continue;
        for (ll vv = v * 2LL; vv <= n; vv += v) {
            is_not_prime[vv] = true;
        }
    }

    srep (v, 2LL, n+1LL) if (!is_not_prime[v]) primes.insert(v);

    delete [] is_not_prime;
}

ull combination(ll x, ll y)
{
    if (y > x / 2LL) y = x - y;

    ull ret = 1LL;
    for (ll i = 0LL; i < y; ++i) {
        ret *= x--;
        ret /= (i + 1LL);
    }

    return ret;
}




// #define MOD   (1000000007LL)
// #define MOD   (998244353LL)

// O(log(exp))
ll mod_pow(ll base, ll exp, ll mod)
{
    ll ret = 1LL;

    for ( ; exp; ) {
        if (exp & 1LL) {
            ret *= base;
            ret %= mod;
        }
        base = (base * base) % mod;
        exp >>= 1;
    }

    return ret;
}

ull mod_combination(ll x, ll y, ll mod)
{
    if (y > x / 2LL) y = x - y;

    ll ret = 1;

    for (ll i = 0LL; i < y; ++i) {
        ret = (ret * x--) % mod;
        ret = (ret * mod_inv(i + 1LL, mod)) % mod;
    }

    return ret;
}

void make_factorial(ll factorials[], ll inv_factorials[], ll size, ll mod)
{
    factorials[0] = 1LL;
    srep (i, 1, size) {
        factorials[i] = factorials[i-1] * i;
        factorials[i] %= mod;
    }

    inv_factorials[size-1] = mod_inv(factorials[size-1], mod);
    rrep (i, size-1) {
        inv_factorials[i] = inv_factorials[i+1LL] * (i+1LL);
        inv_factorials[i] %= mod;
    }
}

#if 0
S ll factorials[200010];
S ll inv_factorials[200010];
make_factorial(factorials, inv_factorials, 200010, MOD);
#endif




///////////////////////////////////////////////////////////////////////////////

void _main();
int main()
{
    cout << setprecision(12);
#if !DEBUG
    ios::sync_with_stdio(false);
    cin.tie(0);
#endif

    _main();
    return 0;
}

// p must be prime
ll mod_root(ll p)
{
    if (p == 2) return 1;
    if (p == 3) return 2;

    V<ll> flist;
    get_factors(flist, p-1LL);
    set<ll> fs;
    for (auto f : flist) if (f > 1) fs.insert(f);

    srep (a, 1, p) {
        bool ok = true;
        for (auto f : fs) {
            if (mod_pow(a, (p-1LL) / f, p) == 1) {
                ok = false;
                break;
            }
        }
        if (ok) return a;
    }

    assert(false);
}

void solve()
{
    ll v; llin(v);
    ll x; llin(x);
    ll p = v * x + 1LL;

    ll root = mod_root(p);
    debug_printf("---- root=%lld\n", root);
    ll base = mod_pow(root, v, p);
    V<ll> anslist;
    rep (xx, x) anslist.pb(mod_pow(base, xx, p));
    sort(anslist);
    for (auto ans : anslist) cout << ans << " ";
    cout << endl;
}

void _main()
{
    ll t; llin(t);
    rep (i, t) solve();
}

///////////////////////////////////////////////////////////////////////////////
0