#1008359 (C++23) No.2880 Max Sigma Mod

提出ソース

結果

問題	No.2880 Max Sigma Mod
ユーザー	t9unkubj
提出日時	2024-09-08 13:01:29
言語	C++23 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	1,917 ms / 3,000 ms
コード長	6,027 bytes
コンパイル時間	3,732 ms
コンパイル使用メモリ	262,276 KB
実行使用メモリ	6,948 KB
最終ジャッジ日時	2024-09-08 13:01:59
合計ジャッジ時間	29,842 ms
ジャッジサーバーID （参考情報）	judge1 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,812 KB
testcase_01	AC	2 ms 6,940 KB
testcase_02	AC	3 ms 6,940 KB
testcase_03	AC	2 ms 6,944 KB
testcase_04	AC	2 ms 6,944 KB
testcase_05	AC	1,507 ms 6,944 KB
testcase_06	AC	1,166 ms 6,948 KB
testcase_07	AC	705 ms 6,944 KB
testcase_08	AC	529 ms 6,940 KB
testcase_09	AC	248 ms 6,944 KB
testcase_10	AC	495 ms 6,944 KB
testcase_11	AC	817 ms 6,940 KB
testcase_12	AC	231 ms 6,940 KB
testcase_13	AC	58 ms 6,944 KB
testcase_14	AC	718 ms 6,944 KB
testcase_15	AC	1,914 ms 6,940 KB
testcase_16	AC	1,917 ms 6,944 KB
testcase_17	AC	1,913 ms 6,940 KB
testcase_18	AC	1,916 ms 6,940 KB
testcase_19	AC	1,901 ms 6,940 KB
testcase_20	AC	12 ms 6,940 KB
testcase_21	AC	333 ms 6,940 KB
testcase_22	AC	801 ms 6,940 KB
testcase_23	AC	11 ms 6,940 KB
testcase_24	AC	438 ms 6,944 KB
testcase_25	AC	1,461 ms 6,944 KB
testcase_26	AC	1,561 ms 6,940 KB
testcase_27	AC	1,529 ms 6,940 KB
testcase_28	AC	1,766 ms 6,944 KB
testcase_29	AC	564 ms 6,944 KB
testcase_30	AC	2 ms 6,944 KB
testcase_31	AC	2 ms 6,940 KB
testcase_32	AC	2 ms 6,940 KB
testcase_33	AC	2 ms 6,944 KB
testcase_34	AC	2 ms 6,944 KB
testcase_35	AC	2 ms 6,940 KB
testcase_36	AC	2 ms 6,940 KB
testcase_37	AC	2 ms 6,940 KB
testcase_38	AC	2 ms 6,940 KB
testcase_39	AC	2 ms 6,940 KB
testcase_40	AC	2 ms 6,944 KB
testcase_41	AC	2 ms 6,944 KB
testcase_42	AC	2 ms 6,940 KB
testcase_43	AC	2 ms 6,940 KB
testcase_44	AC	2 ms 6,940 KB
testcase_45	AC	2 ms 6,940 KB
testcase_46	AC	2 ms 6,940 KB
testcase_47	AC	2 ms 6,944 KB
testcase_48	AC	2 ms 6,944 KB
testcase_49	AC	2 ms 6,944 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#ifdef t9unkubj
#include"template.h"
//#include"template_no_debug.h"
#else
#undef _GLIBCXX_DEBUG
#pragma GCC optimize("O3")
#define dbg(...) 199958
using namespace std;
#include<bits/stdc++.h>
using uint=unsigned;
using ll=long long;
using ull=unsigned long long;
using ld=long double;
using pii=pair<int,int>;
using pll=pair<ll,ll>;
template<class T>using vc=vector<T>;
template<class T>using vvc=vc<vc<T>>;
template<class T>using vvvc=vvc<vc<T>>;
using vi=vc<int>;
using vvi=vc<vi>;
using vvvi=vc<vvi>;
using vl=vc<ll>;
using vvl=vc<vl>;
using vvvl=vc<vvl>;
template<class T>using smpq=priority_queue<T,vector<T>,greater<T>>;
template<class T>using bipq=priority_queue<T>;
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
#define REP(i,j,n) for(ll i=(j);i<(ll)(n);i++)
#define DREP(i,n,m) for(ll i=(n);i>=(m);i--)
#define drep(i,n) for(ll i=((n)-1);i>=0;i--)
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define fi first
#define se second
#define is insert
#define bg begin()
#define ed end()
void scan(int&a) { cin >> a; }
void scan(ll&a) { cin >> a; }
void scan(string&a) { cin >> a; }
void scan(char&a) { cin >> a; }
void scan(uint&a) { cin >> a; }
void scan(ull&a) { cin >> a; }
void scan(bool&a) { cin >> a; }
void scan(ld&a){ cin>> a;}
template<class T> void scan(vector<T>&a) { for(auto&x:a) scan(x); }
void read() {}
template<class Head, class... Tail> void read(Head&head, Tail&... tail) { scan(head); read(tail...); }
#define INT(...) int __VA_ARGS__; read(__VA_ARGS__);
#define LL(...) ll __VA_ARGS__; read(__VA_ARGS__);
#define ULL(...) ull __VA_ARGS__; read(__VA_ARGS__);
#define STR(...) string __VA_ARGS__; read(__VA_ARGS__);
#define CHR(...) char __VA_ARGS__; read(__VA_ARGS__);
#define DBL(...) double __VA_ARGS__; read(__VA_ARGS__);
#define LD(...) ld __VA_ARGS__; read(__VA_ARGS__);
#define VC(type, name, ...) vector<type> name(__VA_ARGS__); read(name);
#define VVC(type, name, size, ...) vector<vector<type>> name(size, vector<type>(__VA_ARGS__)); read(name);
void print(int a) { cout << a; }
void print(ll a) { cout << a; }
void print(string a) { cout << a; }
void print(char a) { cout << a; }
void print(uint a) { cout << a; }
void print(bool a) { cout << a; }
void print(ull a) { cout << a; }
void print(double a) { cout << a; }
void print(ld a){ cout<< a; }
template<class T> void print(vector<T>a) { for(int i=0;i<(int)a.size();i++){if(i)cout<<" ";print(a[i]);}cout<<endl;}
void PRT() { cout <<endl; return ; }
template<class T> void PRT(T a) { print(a); cout <<endl; return; }
template<class Head, class... Tail> void PRT(Head head, Tail ... tail) { print(head); cout << " "; PRT(tail...); return; }
template<class T,class F>
bool chmin(T &x, F y){
    if(x>y){
        x=y;
        return true;
    }
    return false;
}
template<class T, class F>
bool chmax(T &x, F y){
    if(x<y){
        x=y;
        return true;
    }
    return false;
}
void YesNo(bool b){
    cout<<(b?"Yes":"No")<<endl;
}
void Yes(){
    cout<<"Yes"<<endl;
}
void No(){
    cout<<"No"<<endl;
}
template<class T>
int popcount(T n){
    return __builtin_popcountll(n);
}
template<class T>
T sum(vc<T>&a){
    return accumulate(all(a),T(0));
}
template<class T>
T max(vc<T>&a){
    return *max_element(all(a));
}
template<class T>
T min(vc<T>&a){
    return *min_element(all(a));
}
template<class T>
void unique(vc<T>&a){
    a.erase(unique(all(a)),a.end());
}
vvi readgraph(int n,int m,int off = -1){
    vvi g(n);
    rep(i, m){
        int u,v;
        cin>>u>>v;
        u+=off,v+=off;
        g[u].push_back(v);
        g[v].push_back(u);
    }
    return g;
}
vvi readtree(int n,int off=-1){
    return readgraph(n,n-1,off);
}
template<class T>
vc<T> presum(vc<T> &a){
    vc<T> ret(a.size()+1);
    rep(i,a.size())ret[i+1]=ret[i]+a[i];
    return ret;
}
template<class T, class F>
vc<T> &operator+=(vc<T> &a,F b){
    for (auto&v:a)v += b;
    return a;
}
template<class T, class F>
vc<T> &operator-=(vc<T>&a,F b){
    for (auto&v:a)v-=b;
    return a;
}
template<class T, class F>
vc<T> &operator*=(vc<T>&a,F b){
    for (auto&v:a)v*=b;
    return a;
}
#endif
double pass_time=0;

//Σ{0<=i<=n-1}floor((a*i+b)/m)
//基本的なアイデアは寄与を考えて小さいケースへ互除法的に帰着するというもの
template<class T=long long>
//a*b+m,n^2がオーバーフローしないように!
//つまりa=b=1e9などでintにしたらダメ
T floor_sum(T n,T m,T a,T b){
    if(a==0)return b/m*n;
    T ans=0;
    if(a>=m){
        ans+=(n*(n-1))/2*(a/m);
        a%=m;
    }
    if(b>=m){
        ans+=n*(b/m);
        b%=m;
    }
    T y=(a*n+b)/m;
    T z=(a*n+b)%m;
    
    if(y==0)return ans;

    return floor_sum(y,a,m,z)+ans;
}
#define int ll
ll sm(ll l){
    return (l*(l+1))/2;
}
ll solve(ll n,ll m){
    ll ans=0;
    for(int i=1;i<=n;i++){
        ll now=0;
        ll last=0;
        for(int j=1;j*j<=i&&j<=m;j++){
            now+=ll(i/j)*j;
            last=i/j;
        }
        for(int j=1;j<last;j++){
            int l=i/(j+1),r=i/j;
            //(l,r]
            chmin(r,m);
            chmin(r,i);
            if(l>=r)continue;
            chmax(l,0);
            ll w=r-l;
            dbg(l,r);
            now+=(sm(r)-sm(l))*j;
        }
        //[M+1,n]
        chmax(ans,ll(i)*m-now);
    }
    return ans;
}
ll brute(ll n,ll m){
    ll ans=0;
    REP(i,1,n+1){
        ll res=0;
        REP(j,1,m+1){
            res+=i%j;
        }
        chmax(ans,res);
    }
    return ans;
}
signed main(){
    #ifdef t9unkubj
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
    #endif
    cin.tie(0)->sync_with_stdio(0);
    pass_time=clock();
    int t=1;
    //cin>>t;
    /*while(t++){
        int n=(rand()%int(200))+1,m=rand()%int(200)+1;
        if(solve(n,m)!=brute(n,m)){
            dbg(n,m,solve(n,m),brute(n,m));
            return 0;
        }
    }*/
   INT(n,m);
   PRT(solve(n,m));
    pass_time=clock()-pass_time;
    dbg(pass_time/CLOCKS_PER_SEC);
}

yukicoder

結果

テストケース

ソースコード