結果

問題 No.2952 Invision of Multiples
ユーザー t9unkubjt9unkubj
提出日時 2024-09-13 21:32:01
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2,427 ms / 4,000 ms
コード長 11,401 bytes
コンパイル時間 4,040 ms
コンパイル使用メモリ 275,308 KB
実行使用メモリ 192,328 KB
最終ジャッジ日時 2024-09-13 21:33:31
合計ジャッジ時間 77,483 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 782 ms
191,816 KB
testcase_01 AC 753 ms
191,816 KB
testcase_02 AC 1,400 ms
191,816 KB
testcase_03 AC 772 ms
191,744 KB
testcase_04 AC 795 ms
191,816 KB
testcase_05 AC 766 ms
191,692 KB
testcase_06 AC 762 ms
191,816 KB
testcase_07 AC 785 ms
191,820 KB
testcase_08 AC 773 ms
191,820 KB
testcase_09 AC 1,659 ms
192,000 KB
testcase_10 AC 1,723 ms
192,072 KB
testcase_11 AC 1,470 ms
192,068 KB
testcase_12 AC 1,446 ms
192,076 KB
testcase_13 AC 1,441 ms
192,072 KB
testcase_14 AC 1,772 ms
192,128 KB
testcase_15 AC 2,427 ms
192,128 KB
testcase_16 AC 2,198 ms
192,328 KB
testcase_17 AC 2,061 ms
192,196 KB
testcase_18 AC 2,105 ms
192,204 KB
testcase_19 AC 1,696 ms
192,196 KB
testcase_20 AC 1,719 ms
192,204 KB
testcase_21 AC 1,672 ms
192,076 KB
testcase_22 AC 1,675 ms
192,076 KB
testcase_23 AC 1,686 ms
192,256 KB
testcase_24 AC 1,688 ms
192,076 KB
testcase_25 AC 1,706 ms
192,196 KB
testcase_26 AC 1,688 ms
192,256 KB
testcase_27 AC 2,324 ms
192,072 KB
testcase_28 AC 2,279 ms
192,204 KB
testcase_29 AC 2,292 ms
192,256 KB
testcase_30 AC 2,333 ms
192,072 KB
testcase_31 AC 2,287 ms
192,072 KB
testcase_32 AC 2,276 ms
192,200 KB
testcase_33 AC 1,522 ms
192,076 KB
testcase_34 AC 1,568 ms
192,200 KB
testcase_35 AC 1,521 ms
192,200 KB
testcase_36 AC 1,938 ms
192,196 KB
testcase_37 AC 1,928 ms
192,200 KB
testcase_38 AC 1,920 ms
192,200 KB
testcase_39 AC 1,530 ms
192,200 KB
testcase_40 AC 1,549 ms
192,200 KB
testcase_41 AC 1,539 ms
192,200 KB
testcase_42 AC 1,545 ms
192,128 KB
testcase_43 AC 1,519 ms
192,200 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
template<class T,class G,class F>
F modpow(T x,G p,F m){
    T ret=1%m;
    x%=m;
    while(p){
        if(p&1)ret=(1ll*ret*x)%m;
        x=(1ll*x*x)%m;
        p>>=1;
    }
    return ret;
}
template<class T>
T extgcd(T a,T b,T&x,T&y){//ax+by=gcd(a,b)となるようなもの
    if(b==0){
        x=1;
        y=0;
        return a;
    }else{
        T res=extgcd(b,a%b,y,x);
        y-=(a/b)*x;
        return res;
    }
}
template<class T>
pair<T,T> inv(T x,T m){
    T a1,a2;
    T res=extgcd(x,m,a1,a2);
    return {a1,m/res};
}
constexpr int mod=998244353;
struct mint{
    long long val;
    inline long long fast(long long x){
        if(x<mod&&x>=0)return x;
        x%=mod;
        if(x<0)x+=mod;
        return x;
    }
    mint():val(0){}
    mint(long long val):val(fast(val)){}
    mint power(long long m)const {
        mint res(1);
        mint ret(*this);
        while(m){
           if(m&1)res*=ret;
           ret*=ret;
           m>>=1;
        }
        return res;
    }
    mint& operator++() {
        val++;
        if (val == mod) val = 0;
        return *this;
    }
    mint& operator--() {
        if (val == 0)val=mod;
        val--;
        return *this;
    }
    mint operator++(int) {
        mint result = *this;
        ++*this;
        return result;
    }
    mint operator--(int) {
        mint result = *this;
        --*this;
        return result;
    }
    mint operator-() const { 
      return mint(-val);
    }
    friend mint operator +(const mint&a,const mint&b) noexcept{
        return mint(a)+=b;
    }
    friend mint operator -(const mint&a,const mint&b) noexcept{
        return mint(a)-=b;
    }
    friend mint operator *(const mint&a,const mint&b) noexcept{
        return mint(a)*=b;
    }
    friend mint operator /(const mint&a,const mint&b) noexcept{
        return mint(a)/=b;
    }
    mint& operator+=(const mint&a)noexcept{
        val+=a.val;
        if(val>=mod)val-=mod;
        return *this;
    }
    mint& operator-=(const mint&a)noexcept{
        val-=a.val;
        if(val<0)val+=mod;
        return *this;
    }
    mint& operator*=(const mint&a){
        val*=a.val;
        val=fast(val);
        return *this;
    }
    mint& operator/=(const mint&a){
        val*=inv<long long>(a.val,mod).first;
        val=fast(val);
        return *this;
    }
    bool operator == (const mint&x)const noexcept{
        return this->val==x.val;
    }
    bool operator != (const mint&x)const noexcept{
        return this->val!=x.val;
    }
    friend ostream& operator << (ostream &os, const mint &x) noexcept {
        return os << x.val;
    }
    friend istream& operator >> (istream &is,  mint &x) noexcept {
        long long v;
        is >> v;
        x=mint(v);
        return is;
    }
};
vector<mint>fact(1,1),invfact(1,1);
void build(int n){
    if(n<(int)fact.size())return;
    fact=invfact=vector<mint>(n+1);
    fact[0]=1;
    for(int i=1;i<=n;i++)fact[i]=fact[i-1]*i;
    invfact[n]=(1/fact[n]);
    for(int i=n-1;i>=0;i--)invfact[i]=invfact[i+1]*(i+1);
}
mint C(int a,int b){//aCb
    if(a<0||b<0||a-b<0)return mint(0);
    while((int)fact.size()<=a){
        fact.push_back(fact.back()*(fact.size()));
    }
    while((int)invfact.size()<=a){
        invfact.push_back(invfact.back()/invfact.size());
    }
    return fact[a]*invfact[b]*invfact[a-b];
}
mint P(int a,int b){
    if(a<b||b<0)return 0;
    return fact[a]*invfact[a-b];
}
//a個のものからb個を重複を許して選ぶ
mint H(int a,int b){
    return C(a+b-1,b);
}
void print(mint a) { cout << a; }
template<class T>
pair<T,T> mod_solve(T a,T b,T m){//ax=b mod mとなるxを返す
    a%=m,b%=m;if(a<0)a+=m;if(b<0)b+=m;
    T g=gcd(gcd(a,b),m);
    a/=g,b/=g,m/=g;
    if(gcd(a,m)>1)return {-1,-1};
    return {(inv(a,m).first*b)%m,inv(a,m).second};
}
#include<vector>
template <class T = long long >
struct BIT {
    int n;
    std::vector<T>data;
    void add(int i,T x){
        assert(i>=0&&i<n);
        i++;
        for(;i<=n;){
            data[i - 1] += x;
            i += (i & -i);
        }
        return;
    }   
    T internal_sum(int r){
        T ret = 0;
        for(;r;){
            ret += data[r - 1];
            r -= (r & -r);
        }
        return ret;
    }
    T sum(int l,int r){
        return internal_sum(r)-internal_sum(l);
    }
    BIT():n(0), data(0) {}
    BIT(int N):n(N), data(N) {}
    BIT(std::vector<T>& a){
        n = a.size();
        data.assign(n, 0);
        for (int i = 0; i < n; i++)add(i, a[i]);
    }
    int lower_bound(T w) { 
    if (w <= 0) {
        return 0;
    } else {
        int x = 0, r = 1;
        while (r < n) r = r << 1;
        for (int len = r; len > 0; len = len >> 1) { 
            if (x + len < n && data[x + len] < w) {
                w -= data[x + len];
                x += len;
            }
        }
        return x;
    }
}
};
/*
平方分割
d[i]<Bとなるようなものについては前計算をそのまま利用
d[i]>=Bとなるものについてはd[i]<Bの前計算+
BITで点の個数を数えるクエリに変更
*/
constexpr int N=6e4+1;
constexpr int B=200;
mint f1[B][N];
mint f2[B][N];
mint sm[B];
int cnt1[N];
int cnt2[N];
void solve(){
    INT(n,m);
    VC(int,d,n);
	assert(n>=1&&n<=5e4);
	assert(m>=1&&m<=5e4);
    vc<mint>invs(N);
    REP(i,1,N){
        if(m/i==0)continue;
        invs[i]=mint(1)/(m/i);
    }
    
    mint prod=1;
    rep(i,n)prod*=(m/d[i]);
	rep(i,n)assert(1<=d[i]&&d[i]<=m);
    REP(i,1,B){
        REP(j,1,N){
            f1[i][j]=floor_sum<ull>(m/i,j,i,i-1)%mod;
            f2[i][j]=floor_sum<ull>(m/j,i,j,j-1)%mod;
        }
    }
    BIT<mint>bit(N);
    mint ans=0;

    for(auto&x:d)cnt2[x]++;

    for(auto&x:d){
        if(x<B){
            for(int i=1;i<B;i++){
                ans+=f1[i][x]*cnt1[i]*prod*invs[i]*invs[x];
            } 
            cnt1[x]++;
            cnt2[x]--;
            continue;
        }else{
            for(int i=1;i<B;i++){
                ans+=f1[i][x]*cnt1[i]*prod*invs[i]*invs[x];
            }
            for(int i=1;i<B;i++){
                ans+=f2[i][x]*cnt2[i]*prod*invs[i]*invs[x];
            }

            for(int k=x;k<=m;k+=x){
                ans+=bit.sum(k+1,N)*invs[x];
            }
            for(int k=x;k<=m;k+=x){
                bit.add(k,prod*invs[x]);
            }
        }
    }

    PRT(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--)solve();
    pass_time=clock()-pass_time;
    dbg(pass_time/CLOCKS_PER_SEC);
}
0