Submission #992648

#	Time	Username	Problem	Language	Result	Execution time	Memory
992648		PedroBigMan	Fish (IOI08_fish)	C++14	95 / 100	319 ms	65536 KiB

fish

#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
#pragma GCC optimize("Ofast")
#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <string>
#include <map>
#include <unordered_map>
#include <set>
#include <unordered_set>
#include <queue>
#include <deque>
#include <list>
#include <iomanip>
#include <stdlib.h>
#include <time.h>
#include <cstring>
using namespace std;
typedef int ll;
typedef unsigned long long int ull;
typedef long double ld;
#define REP(i,a,b) for(ll i=(ll) a; i<(ll) b; i++)
#define pb push_back
#define mp make_pair
#define pl pair<ll,ll>
#define ff first
#define ss second
#define whole(x) x.begin(),x.end()
#define DEBUG(i) cout<<"Pedro "<<i<<endl
#define INF 1000000000000000000LL
#define EPS ((ld)0.00000000001)
#define pi ((ld)3.141592653589793)
#define VV(vvvv,NNNN,xxxx); REP(iiiii,0,NNNN) {vvvv.pb(xxxx);}
ll mod=1000000007;
 
template<class A=ll> 
void Out(vector<A> a) {REP(i,0,a.size()) {cout<<a[i]<<" ";} cout<<endl;}
 
template<class A=ll>
void In(vector<A> &a, ll N) {A cur; REP(i,0,N) {cin>>cur; a.pb(cur);}} 
 
class ST
{	
    public:
    ll N;
    
    vector<ll> p;
    vector<pl> range;
    ll thisans;
    
    ST() {N=0;}
    ST(vector<ll> arr)
    {
        N = (ll) 1<<(ll) ceil(log2(arr.size()));
        REP(i,0,2*N) {range.pb(mp(0,0));}
        REP(i,0,N) {p.pb(1);}
        REP(i,0,arr.size()) {p.pb(arr[i]); range[i+N]=mp(i,i);}
        REP(i,arr.size(),N) {p.pb(1); range[i+N]=mp(i,i);}
        ll cur = N-1;
        while(cur>0)
        {
            p[cur]=(p[2*cur]*p[2*cur+1])%mod;
            range[cur]=mp(range[2*cur].ff,range[2*cur+1].ss);
            cur--;
        }
    }
    
    void query(ll a,ll b, ll c=1) //range [a,b], current node. initially: query(a,b)
    {
        if(c==1) {thisans=1;}
		if(a>b) {return ;}
        ll x=range[c].ff; ll y=range[c].ss;
        if(y<a || x>b) {return ;}
        if(x>=a && y<=b) {thisans*=p[c]; thisans%=mod;return ;}
        query(a,b,2*c);query(a,b,2*c+1);
    }
    
    void update(ll s, ll a, ll b, ll c=1) //update LV, range [a,b], current node, current range. initially: update(s,a,b)
    {
		if(a>b) {return;}
        ll x=range[c].ff; ll y=range[c].ss;
        if(y<a || x>b) {return ;}
        if(x>=a && y<=b) 
        {
            p[c]+=s; p[c]%=mod;
            return;
        }
        update(s,a,b,2*c); update(s,a,b,2*c+1);
        p[c]=(p[2*c]*p[2*c+1])%mod;
    }
};
 
int main()
{
    ios_base::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
	cout.precision(20);
	ll N, K, M; cin>>N>>K>>M; mod=M;
    vector<pl> p(N,{-1,-1});
    REP(i,0,N) {cin>>p[i].ff>>p[i].ss; p[i].ss--;}
    sort(whole(p));
    vector<pl> to_order(K,{-1,-1}); REP(i,0,K) {to_order[i].ss=i;}
    REP(i,0,N) {to_order[p[i].ss].ff=p[i].ff;}
    sort(whole(to_order));
    vector<ll> limit(K,-1), m(K,-1);
    REP(i,0,K) {m[to_order[i].ss]=i; limit[i]=to_order[i].ff/2;}
    vector<ll> L(N,-1); REP(i,0,N) {L[i]=p[i].ff;} REP(i,0,K) {limit[i]=(ll) (upper_bound(whole(L),limit[i])-L.begin())-1;}
    REP(i,0,N) {p[i].ss=m[p[i].ss];}
    vector<ll> gem(N,-1); REP(i,0,N) {gem[i]=p[i].ss;}
    vector<vector<ll> > pos(K,vector<ll>()); REP(i,0,N) {pos[gem[i]].pb(i);}
    vector<ll> nxt(K,-1); 
    REP(i,0,K) {nxt[i]=*upper_bound(whole(pos[i]),limit[i]);}
    vector<ll> xx(K,1); ST S(xx);
    ll gem_ind = 0; ll ans = 0; ll ind=0; ll thisgem;
    while(gem_ind<K && limit[gem_ind]==-1) {gem_ind++; ans++; ans%=mod;}
    ll G; ll B; vector<ll> oc(K,0); ll val1,val2;
    while(ind<N && gem_ind<K)
    {
        thisgem = gem[ind]; S.update(1,thisgem,thisgem); oc[thisgem]++;
        while(gem_ind<K && limit[gem_ind]==ind) 
        {
            G=gem_ind;
            B = (ll) (lower_bound(whole(limit),nxt[G]) - limit.begin()) -1;
            S.query(0,G-1); val1=S.thisans; S.query(G+1,B); val2=S.thisans;
            ans+=val1*(val2+oc[G]%mod); ans%=mod; 
            gem_ind++;
        }
        ind++;
    }
    if(ans<0LL) {ans+=mod;} cout<<ans<<endl;
    return 0;
}

Compilation message (stderr)

fish.cpp:1: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
    1 | #pragma GCC optimization ("O3")
      | 
fish.cpp:2: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
    2 | #pragma GCC optimization ("unroll-loops")
      |

• Subtask 1: 4 / 4
• Subtask 2: 4 / 4
• Subtask 3: 4 / 4
• Subtask 4: 4 / 4
• Subtask 5: 4 / 4
• Subtask 6: 5 / 5
• Subtask 7: 6 / 6
• Subtask 8: 6 / 6
• Subtask 9: 6 / 6
• Subtask 10: 6 / 6
• Subtask 11: 7 / 7
• Subtask 12: 7 / 7
• Subtask 13: 7 / 7
• Subtask 14: 5 / 5
• Subtask 15: 5 / 5
• Subtask 16: 5 / 5
• Subtask 17: 5 / 5
• Subtask 18: 5 / 5
• Subtask 19: 0 / 5