Submission #262596

# Submission time Handle Problem Language Result Execution time Memory
262596 2020-08-13T04:58:03 Z ToMmyDong Strange Device (APIO19_strange_device) C++11
100 / 100
1096 ms 69164 KB
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#ifdef tmd
typedef __int128 lll;
#else
typedef __int128 lll;
#endif
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<int, ll> pil;
typedef pair<ll, int> pli;
typedef pair<double, double> pdd;
#define SQ(i) ((i)*(i))
#define MEM(a, b) memset(a, (b), sizeof(a))
#define SZ(i) static_cast<int>(i.size())
#define FOR(i, j, k, in) for (int i=j; i < (k) ; i+=in)
#define RFOR(i, j, k, in) for (int i=j; i >= (k) ; i-=in)
#define REP(i, j) FOR(i, 0, j, 1)
#define REP1(i, j) FOR(i, 1, j+1, 1)
#define RREP(i, j) RFOR(i, j, 0, 1)
#define ALL(_a) _a.begin(), _a.end()
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define X first
#define Y second
template<typename T, typename S>
istream &operator >> (istream &is, pair<T, S> &p) {
    return is >> p.X >> p.Y;
}
#ifdef tmd
#define TIME(i) Timer i(#i)
#define debug(...) do {\
    fprintf(stderr, "%s - %d (%s) = ", __PRETTY_FUNCTION__, __LINE__, #__VA_ARGS__);\
    _do(__VA_ARGS__);\
}while(0)
template<typename T> void _do(T &&_x) {cerr<<_x<<endl;}
template<typename T,typename ...S> void _do(T &&_x, S &&..._t) {cerr <<_x <<" ,"; _do(_t...);}
template<typename _a,typename _b> ostream& operator << (ostream &_s, const pair<_a, _b> &_p) {return _s << "(" << _p.X << "," << _p.Y << ")";}
template<typename It> ostream& _OUTC(ostream &_s, It _ita, It _itb)
{
    _s << "{";
    for (It _it=_ita; _it != _itb; _it++) {
        _s <<(_it == _ita?"":",")<< *_it;
    }
    _s << "}";
    return _s;
}
template<typename _a> ostream &operator << (ostream &_s, const vector<_a> &_c) {return _OUTC(_s,ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const array<_a,2> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const array<_a,3> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const array<_a,4> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const array<_a,5> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const set<_a> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a> ostream &operator << (ostream &_s, const deque<_a> &_c) {return _OUTC(_s, ALL(_c));}
template<typename _a,typename _b> ostream &operator << (ostream &_s, const map<_a,_b> &_c) {return _OUTC(_s,ALL(_c));}
template<typename _t> void pary(_t _a,_t _b){_OUTC(cerr,_a,_b);cerr<<endl;}
#define IOS()
class Timer {
private:
    string scope_name;
    chrono::high_resolution_clock::time_point start_time;
public:
    Timer (string name) : scope_name(name) {
        start_time = chrono::high_resolution_clock::now();
    }
    ~Timer () {
        auto stop_time = chrono::high_resolution_clock::now();
        auto length = chrono::duration_cast<chrono::microseconds>(stop_time - start_time).count();
        double mlength = double(length) * 0.001;
        debug(scope_name, mlength);
    }
};
#else
#define TIME(i)
#define debug(...)
#define pary(...)
#define endl '\n'
#define IOS() ios_base::sync_with_stdio(0);cin.tie(0)
#endif

const ll MOD = 1000000007;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int iNF = 0x3f3f3f3f;
const ll MAXN = 2e5 + 5;

void extGCD(lll A,lll B,lll &x,lll &y) { // A p coprime
    if (B == 0) {
        x = 1;
        y = 0;
        assert(A == 1);
        return;
    }
    lll xx,yy;
    extGCD(B,A%B,xx,yy);
    x = yy;
    y = xx - A/B*yy;
    return;
}

lll ext_inv (lll a, lll p) { // a, p co-prime
    lll x, y;
    extGCD(a,p, x, y);
    x %= p;
    if (x < 0) {
        x += p;
    }
    assert(a * x % p);
    return x;
}

lll di (lll a, lll b, lll md) {
    if (a == 0) {
        return 0;
    }
    assert(b != 0);
    lll cd = __gcd( __gcd(a, b), md);
    a /= cd;
    b /= cd;
    md /= cd;
    
    return ext_inv(b,md) * a % md;
}
signed main () {
    TIME(main);
    IOS();

    debug(RAND_MAX);
    //while (true) {

    ll n, a, b;
#ifdef tmd
    n = rand() % 10000+1;
    //a = ( ( ll(rand())<<32 ) + rand() ) % ll( 1e18 );
    a = rand()%100 + 1;
    b = rand()%3 + 1;
#else
    cin >> n >> a >> b;
#endif
    debug(n, a, b);

    lll sz = b*lll(a/__gcd(a,b+1));
    if (sz == 1) {
        cout << 1 << endl;
        //continue;
        return 0;
    }

    vector<pair<lll,lll>> seg;
    lll g = (b+1) % a;

    lll LS = 0;
    bool gg= false;
    for (int i=0; i<n; i++) {
        ll l, r;
#ifdef tmd
        l = LS + rand();
        r = l + rand() % (sz-1);
        LS = r;
#else        
        cin >> l >> r;
#endif
        if (r-l+1 >= sz) {
            cout << ll(sz) << endl;
            gg = true;
            break;
        } else {
            lll sft = (a+l%a+(l/b)%a-l%b%a)%a;
            lll rw = di(sft, g, a);
            lll st = l%b + b*rw;
            lll ed = st + r-l;
            st %= sz;
            ed %= sz;
            if (ed < st) {
                seg.eb(0, ed);
                seg.eb(st, sz-1);
            } else {
                seg.eb(st, ed);
            }

        }
    }
    if (gg) {
        //continue;
        return 0;
    }
    sort(ALL(seg));

    lll lst = -1;
    lll ans = 0;
    for (auto p : seg) {
        ans += max(lll(0), p.Y - max(lst+1, p.X) + 1);
        lst = max(lst, p.Y);
    }

    cout << ll(ans) << endl;
    //}


    return 0;
}

Compilation message

strange_device.cpp: In function 'int main()':
strange_device.cpp:153:9: warning: unused variable 'LS' [-Wunused-variable]
  153 |     lll LS = 0;
      |         ^~
# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 9 ms 1152 KB Output is correct
3 Correct 8 ms 1024 KB Output is correct
4 Correct 1 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 0 ms 384 KB Output is correct
7 Correct 0 ms 384 KB Output is correct
8 Correct 1 ms 384 KB Output is correct
9 Correct 0 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 1 ms 384 KB Output is correct
12 Correct 1 ms 384 KB Output is correct
13 Correct 0 ms 384 KB Output is correct
14 Correct 0 ms 384 KB Output is correct
15 Correct 0 ms 384 KB Output is correct
16 Correct 8 ms 1024 KB Output is correct
17 Correct 82 ms 4580 KB Output is correct
18 Correct 0 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 384 KB Output is correct
2 Correct 0 ms 384 KB Output is correct
3 Correct 1 ms 384 KB Output is correct
4 Correct 0 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 1 ms 384 KB Output is correct
3 Correct 1 ms 384 KB Output is correct
4 Correct 1 ms 384 KB Output is correct
5 Correct 427 ms 33364 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 384 KB Output is correct
2 Correct 694 ms 33344 KB Output is correct
3 Correct 670 ms 33344 KB Output is correct
4 Correct 682 ms 33560 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 384 KB Output is correct
2 Correct 694 ms 33344 KB Output is correct
3 Correct 670 ms 33344 KB Output is correct
4 Correct 682 ms 33560 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 724 ms 33468 KB Output is correct
7 Correct 705 ms 33344 KB Output is correct
8 Correct 688 ms 33344 KB Output is correct
9 Correct 811 ms 33344 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 384 KB Output is correct
2 Correct 694 ms 33344 KB Output is correct
3 Correct 670 ms 33344 KB Output is correct
4 Correct 682 ms 33560 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 74 ms 4592 KB Output is correct
7 Correct 73 ms 4588 KB Output is correct
8 Correct 69 ms 4584 KB Output is correct
9 Correct 74 ms 4580 KB Output is correct
10 Correct 86 ms 4588 KB Output is correct
11 Correct 72 ms 4584 KB Output is correct
12 Correct 86 ms 4584 KB Output is correct
13 Correct 93 ms 4584 KB Output is correct
14 Correct 65 ms 4588 KB Output is correct
15 Correct 97 ms 4584 KB Output is correct
16 Correct 87 ms 4588 KB Output is correct
17 Correct 102 ms 4556 KB Output is correct
18 Correct 806 ms 33348 KB Output is correct
19 Correct 868 ms 33344 KB Output is correct
20 Correct 986 ms 33344 KB Output is correct
21 Correct 82 ms 4588 KB Output is correct
22 Correct 72 ms 4684 KB Output is correct
23 Correct 175 ms 16980 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 90 ms 4588 KB Output is correct
3 Correct 97 ms 4568 KB Output is correct
4 Correct 1069 ms 33344 KB Output is correct
5 Correct 98 ms 4584 KB Output is correct
6 Correct 96 ms 4588 KB Output is correct
7 Correct 115 ms 4584 KB Output is correct
8 Correct 98 ms 4592 KB Output is correct
9 Correct 105 ms 4640 KB Output is correct
10 Correct 89 ms 4588 KB Output is correct
11 Correct 77 ms 4716 KB Output is correct
12 Correct 64 ms 4588 KB Output is correct
13 Correct 102 ms 4584 KB Output is correct
14 Correct 942 ms 33288 KB Output is correct
15 Correct 81 ms 4588 KB Output is correct
16 Correct 740 ms 33368 KB Output is correct
17 Correct 784 ms 33520 KB Output is correct
18 Correct 0 ms 384 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 9 ms 1152 KB Output is correct
3 Correct 8 ms 1024 KB Output is correct
4 Correct 1 ms 384 KB Output is correct
5 Correct 0 ms 384 KB Output is correct
6 Correct 0 ms 384 KB Output is correct
7 Correct 0 ms 384 KB Output is correct
8 Correct 1 ms 384 KB Output is correct
9 Correct 0 ms 384 KB Output is correct
10 Correct 0 ms 384 KB Output is correct
11 Correct 1 ms 384 KB Output is correct
12 Correct 1 ms 384 KB Output is correct
13 Correct 0 ms 384 KB Output is correct
14 Correct 0 ms 384 KB Output is correct
15 Correct 0 ms 384 KB Output is correct
16 Correct 8 ms 1024 KB Output is correct
17 Correct 82 ms 4580 KB Output is correct
18 Correct 0 ms 384 KB Output is correct
19 Correct 0 ms 384 KB Output is correct
20 Correct 0 ms 384 KB Output is correct
21 Correct 1 ms 384 KB Output is correct
22 Correct 0 ms 384 KB Output is correct
23 Correct 0 ms 384 KB Output is correct
24 Correct 1 ms 384 KB Output is correct
25 Correct 1 ms 384 KB Output is correct
26 Correct 1 ms 384 KB Output is correct
27 Correct 1 ms 384 KB Output is correct
28 Correct 427 ms 33364 KB Output is correct
29 Correct 0 ms 384 KB Output is correct
30 Correct 694 ms 33344 KB Output is correct
31 Correct 670 ms 33344 KB Output is correct
32 Correct 682 ms 33560 KB Output is correct
33 Correct 0 ms 384 KB Output is correct
34 Correct 724 ms 33468 KB Output is correct
35 Correct 705 ms 33344 KB Output is correct
36 Correct 688 ms 33344 KB Output is correct
37 Correct 811 ms 33344 KB Output is correct
38 Correct 0 ms 384 KB Output is correct
39 Correct 74 ms 4592 KB Output is correct
40 Correct 73 ms 4588 KB Output is correct
41 Correct 69 ms 4584 KB Output is correct
42 Correct 74 ms 4580 KB Output is correct
43 Correct 86 ms 4588 KB Output is correct
44 Correct 72 ms 4584 KB Output is correct
45 Correct 86 ms 4584 KB Output is correct
46 Correct 93 ms 4584 KB Output is correct
47 Correct 65 ms 4588 KB Output is correct
48 Correct 97 ms 4584 KB Output is correct
49 Correct 87 ms 4588 KB Output is correct
50 Correct 102 ms 4556 KB Output is correct
51 Correct 806 ms 33348 KB Output is correct
52 Correct 868 ms 33344 KB Output is correct
53 Correct 986 ms 33344 KB Output is correct
54 Correct 82 ms 4588 KB Output is correct
55 Correct 72 ms 4684 KB Output is correct
56 Correct 175 ms 16980 KB Output is correct
57 Correct 1 ms 384 KB Output is correct
58 Correct 90 ms 4588 KB Output is correct
59 Correct 97 ms 4568 KB Output is correct
60 Correct 1069 ms 33344 KB Output is correct
61 Correct 98 ms 4584 KB Output is correct
62 Correct 96 ms 4588 KB Output is correct
63 Correct 115 ms 4584 KB Output is correct
64 Correct 98 ms 4592 KB Output is correct
65 Correct 105 ms 4640 KB Output is correct
66 Correct 89 ms 4588 KB Output is correct
67 Correct 77 ms 4716 KB Output is correct
68 Correct 64 ms 4588 KB Output is correct
69 Correct 102 ms 4584 KB Output is correct
70 Correct 942 ms 33288 KB Output is correct
71 Correct 81 ms 4588 KB Output is correct
72 Correct 740 ms 33368 KB Output is correct
73 Correct 784 ms 33520 KB Output is correct
74 Correct 0 ms 384 KB Output is correct
75 Correct 1 ms 384 KB Output is correct
76 Correct 0 ms 384 KB Output is correct
77 Correct 0 ms 384 KB Output is correct
78 Correct 1 ms 384 KB Output is correct
79 Correct 9 ms 1408 KB Output is correct
80 Correct 1064 ms 69164 KB Output is correct
81 Correct 1075 ms 51904 KB Output is correct
82 Correct 1080 ms 52568 KB Output is correct
83 Correct 1055 ms 52416 KB Output is correct
84 Correct 1033 ms 52160 KB Output is correct
85 Correct 1096 ms 52200 KB Output is correct
86 Correct 197 ms 26708 KB Output is correct
87 Correct 1 ms 384 KB Output is correct