Submission #529919

#	Time	Username	Problem	Language	Result	Execution time	Memory
529919		i_am_noob	Squirrel (RMI18_squirrel)	C++17	60 / 100	4794 ms	19140 KiB

squirrel

#include <bits/stdc++.h>
using namespace std;
#pragma GCC optimize("Ofast,unroll-loops")

#define ll long long
#define int ll
#define ull unsigned ll
#define ld long double
#define rep(a) rep1(i,a)
#define rep1(i,a) rep2(i,0,a)
#define rep2(i,b,a) for(int i=(b); i<((int)(a)); i++)
#define rep3(i,b,a) for(int i=(b); i>=((int)(a)); i--)
#define chkmin(a,b) (a=min(a,b))
#define chkmax(a,b) (a=max(a,b))
#define all(a) a.begin(),a.end()
#define pii pair<int,int>
#define pb push_back
//#define inf 1010000000
#define inf 4000000000000000000
#define eps 1e-9
#define sz(a) ((int)a.size())
#define pow2(x) (1ll<<(x))
#define ceiling(a,b) (((a)+(b)-1)/(b))
#define print0(a) cout << (a) << ' '
#define ykh mt19937 rng(chrono::steady_clock::now().time_since_epoch().count())
#ifdef i_am_noob
#define bug(...) cerr << "#" << __LINE__ << ' ' << #__VA_ARGS__ << "- ", _do(__VA_ARGS__)
template<typename T> void _do(vector<T> x){for(auto i: x) cerr << i << ' ';cerr << "\n";}
template<typename T> void _do(set<T> x){for(auto i: x) cerr << i << ' ';cerr << "\n";}
template<typename T> void _do(unordered_set<T> x){for(auto i: x) cerr << i << ' ';cerr << "\n";}
template<typename T> void _do(T && x) {cerr << x << endl;}
template<typename T, typename ...S> void _do(T && x, S&&...y) {cerr << x << ", "; _do(y...);}
#else
#define bug(...) 777771449
#endif
template<typename T> void print(T && x) {cout << x << "\n";}
template<typename T, typename... S> void print(T && x, S&&... y) {cout << x << ' ';print(y...);}

const int Mod=1000000007,Mod2=998244353;
const int MOD=Mod2;
template <int mod>
struct Modint{
    int val;
    Modint(int _val=0){val=_val%mod;if(val<0) val+=mod;}
    Modint operator +(const Modint& o) const {
        Modint res;
        res.val=val+o.val;
        if(res.val>=mod) res.val-=mod;
        return res;
    }
    Modint operator +(const int& o) const {return Modint(val+o);}
    Modint operator -() const {
        Modint res;
        res.val=-val;
        if(res.val<0) res.val+=mod;
        return res;
    }
    Modint operator -(const Modint& o) const {
        Modint res;
        res.val=val-o.val;
        if(res.val<0) res.val+=mod;
        return res;
    }
    Modint operator -(const int& o) const {return Modint(val-o);}
    Modint operator *(const Modint& o) const {return Modint(val*o.val);}
    Modint operator *(const int& o) const {return Modint(val*(o%mod));}
    Modint operator +=(const Modint& o){*this=(*this)+o;return *this;}
    Modint operator -=(const Modint& o){*this=(*this)-o;return *this;}
    Modint operator *=(const Modint& o){*this=(*this)*o;return *this;}
    Modint Pow(int b) const {
        Modint tmp(val),ret(1);
        while(b){
            if(b&1) ret*=tmp;
            b>>=1;tmp*=tmp;
        }
        return ret;
    }
    Modint Pow(const Modint& a, int b) const {return a.Pow(b);}
    inline Modint inv() const {return (*this).Pow(mod-2);}
    Modint operator /(const Modint& o) const {return *this*o.inv();}
    Modint operator /(const int& o) const {return *this*Modint(o).inv();}
    bool operator ==(const Modint& o) const {return val==o.val;}
};
template<int mod>
ostream& operator << (ostream& o, Modint<mod> x){return o << x.val;}
template<int mod>
Modint<mod> operator +(const int& x, const Modint<mod>& y){return Modint<mod>(x+y.val);}
template<int mod>
Modint<mod> operator -(const int& x, const Modint<mod>& y){return Modint<mod>(x-y.val);}
template<int mod>
Modint<mod> operator *(const int& x, const Modint<mod>& y){return Modint<mod>(x%mod*y.val);}
#define modint Modint<MOD>
vector<modint> inv,fac,invfac;
void init_comb(int N){
    inv.resize(N),fac.resize(N),invfac.resize(N);
    inv[1]=1,fac[0]=1,invfac[0]=1;
    rep2(i,2,N) inv[i]=inv[MOD%i]*(MOD-MOD/i);
    rep2(i,1,N) fac[i]=fac[i-1]*i;
    rep2(i,1,N) invfac[i]=invfac[i-1]*inv[i];
}
inline modint C(int n, int m){return m>n||m<0?modint(0):fac[n]*invfac[m]*invfac[n-m];}
inline modint H(int n, int m){return C(n+m-1,n);}

const int maxn=100005,maxm=5,maxk=7777714;

//i_am_noob
//#define wiwihorz  
int n,m,q,mobius[maxn],res,barrett[maxn];
vector<vector<int>> factors(maxn),pr(maxn);
bool isp[maxn];

int Div(int x, int y){return (x*barrett[y])>>32;}

int count(int k, int l, int r){
    if(k==0){
        if(l<=1&&r>=1) return 1;
        else return 0;
    }
    int tot=0;
    if((r-l)*sz(pr[k])<=sz(factors[k])){
        bool good[r-l+1];
        rep(r-l+1) good[i]=1;
        for(auto p: pr[k]){
            int x=Div(l,p)*p;
            while(x<l) x+=p;
            while(x<=r) good[x-l]=0,x+=p;
        }
        return count(good,good+r-l+1,1);
    }
    if(l==0){
        if(k==1) tot++;
        l++;
    }
    for(auto i: factors[k]) tot+=mobius[i]*(Div(r,i)-Div(l-1,i));
    return tot;
}

void solve(int x, int y, int f, int dir){
    //0: up, 1: down, 2: left, 3: right
    if(dir==0){
        res+=count(y,x-f,x-1);
        if(f==1) return;
        x-=f;
        res+=count(abs(x-y),x-f/2,x-1)+count(x+y,x-f/2,x-1);
        solve(x-f/2,y-f/2,f/2,0);
        solve(x-f/2,y-f/2,f/2,2);
        solve(x-f/2,y+f/2,f/2,0);
        solve(x-f/2,y+f/2,f/2,3);
    }
    else if(dir==1){
        res+=count(y,x+1,x+f);
        if(f==1) return;
        x+=f;
        res+=count(abs(x-y),x+1,x+f/2)+count(x+y,x+1,x+f/2);
        solve(x+f/2,y-f/2,f/2,1);
        solve(x+f/2,y-f/2,f/2,2);
        solve(x+f/2,y+f/2,f/2,1);
        solve(x+f/2,y+f/2,f/2,3);
    }
    else if(dir==2){
        res+=count(x,y-f,y-1);
        if(f==1) return;
        y-=f;
        res+=count(abs(x-y),y-f/2,y-1)+count(x+y,y-f/2,y-1);
        solve(x-f/2,y-f/2,f/2,2);
        solve(x-f/2,y-f/2,f/2,0);
        solve(x+f/2,y-f/2,f/2,2);
        solve(x+f/2,y-f/2,f/2,1);
    }
    else{
        res+=count(x,y+1,y+f);
        if(f==1) return;
        y+=f;
        res+=count(abs(x-y),y+1,y+f/2)+count(x+y,y+1,y+f/2);
        solve(x-f/2,y+f/2,f/2,3);
        solve(x-f/2,y+f/2,f/2,0);
        solve(x+f/2,y+f/2,f/2,3);
        solve(x+f/2,y+f/2,f/2,1);
    }
}

void orzck(){
    /*
    de[0]=1;
    rep2(i,1,11) de[i]=2*pow2(i)+4*de[i-1];
    rep(11) bug(i,de[i]);
    */
    isp[0]=isp[1]=0;
    rep2(i,2,maxn) isp[i]=1;
    rep2(i,2,maxn) if(isp[i]) for(int j=i*2; j<maxn; j+=i) isp[j]=0;
    mobius[1]=1;
    rep2(i,1,maxn) for(int j=i*2; j<maxn; j+=i) mobius[j]-=mobius[i];
    rep2(i,1,maxn) if(mobius[i]) for(int j=i; j<maxn; j+=i) factors[j].pb(i);
    rep2(i,1,maxn) if(isp[i]) for(int j=i; j<maxn; j+=i) pr[j].pb(i);
    rep2(i,1,maxn) barrett[i]=(pow2(32)+i-1)/i;
    cin >> n >> m >> q;
    while(q--){
        int x,y,f;
        cin >> x >> y >> f;
        x--,y--;
        if(__gcd(x,y)==1) res++;
        solve(x,y,f,0);
    }
    print(res);
}

signed main(){
    ios_base::sync_with_stdio(0),cin.tie(0);
    // #ifdef i_am_noob
    // freopen("input1.txt","r",stdin);
    // freopen("output1.txt","w",stdout);
    // freopen("output2.txt","w",stderr);
    // #endif
    cout << fixed << setprecision(15);
    int t;
    #ifdef wiwihorz
    cin >> t;
    #else
    t=1;
    #endif
    ld start=clock();
    while(t--) orzck();
    bug((clock()-start)/CLOCKS_PER_SEC);
    return 0;
}

Compilation message (stderr)

squirrel.cpp: In function 'int main()':
squirrel.cpp:34:18: warning: statement has no effect [-Wunused-value]
   34 | #define bug(...) 777771449
      |                  ^~~~~~~~~
squirrel.cpp:223:5: note: in expansion of macro 'bug'
  223 |     bug((clock()-start)/CLOCKS_PER_SEC);
      |     ^~~
squirrel.cpp:221:8: warning: unused variable 'start' [-Wunused-variable]
  221 |     ld start=clock();
      |        ^~~~~

• Subtask 1: 15 / 15
• Subtask 2: 10 / 10
• Subtask 3: 25 / 25
• Subtask 4: 10 / 50