結果

問題 No.886 Direct
ユーザー kakao745kakao745
提出日時 2024-10-16 22:43:29
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 278 ms / 4,000 ms
コード長 5,372 bytes
コンパイル時間 6,121 ms
コンパイル使用メモリ 337,616 KB
実行使用メモリ 41,600 KB
最終ジャッジ日時 2024-10-16 22:43:40
合計ジャッジ時間 9,996 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 1 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 2 ms
5,248 KB
testcase_11 AC 2 ms
5,248 KB
testcase_12 AC 1 ms
5,248 KB
testcase_13 AC 1 ms
5,248 KB
testcase_14 AC 2 ms
5,248 KB
testcase_15 AC 2 ms
5,248 KB
testcase_16 AC 2 ms
5,248 KB
testcase_17 AC 2 ms
5,248 KB
testcase_18 AC 3 ms
5,248 KB
testcase_19 AC 2 ms
5,248 KB
testcase_20 AC 3 ms
5,248 KB
testcase_21 AC 3 ms
5,248 KB
testcase_22 AC 2 ms
5,248 KB
testcase_23 AC 94 ms
21,288 KB
testcase_24 AC 106 ms
23,040 KB
testcase_25 AC 59 ms
15,492 KB
testcase_26 AC 83 ms
19,200 KB
testcase_27 AC 247 ms
38,588 KB
testcase_28 AC 210 ms
37,452 KB
testcase_29 AC 244 ms
41,344 KB
testcase_30 AC 243 ms
41,344 KB
testcase_31 AC 249 ms
41,344 KB
testcase_32 AC 278 ms
41,472 KB
testcase_33 AC 247 ms
41,344 KB
testcase_34 AC 239 ms
41,600 KB
testcase_35 AC 241 ms
41,344 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define itn int
#define rep(i,n) for(long long i=0;i<(long long)n;i++)
#define reps(i,n) for(long long i=1;i<=(long long)n;i++)
#define repi(i,n) for(int i=0;i<(int)n;i++)
#define loop(i,l,r) for(long long i=(long long)l;i<=(long long)r;i++)
#define loopi(i,l,r) for(int i=(int)l;i<=(int)r;i++)
#define drep(i,n) for(long long i=(long long)n-1;i>=0;i--)
#define drepi(i,n) for(int i=(int)n-1;i>=0;i--)
#define dreps(i,n) for(int i=(int)n;i>=1;i--)
#define dloop(i,l,r) for(long long i=(long long)l;i>=(long long)r;i--)
#define dloopi(i,l,r) for(int i=(int)l;i>=(int)r;i--)
#define all(v) v.begin(), v.end()
#define rall(v) v.rbegin(), v.rend()
#define yn(x) cout << (x? "Yes":"No") << endl;
#define cou(x) cout << x << endl;
#define emp emplace_back
#pragma GCC target ("avx,avx2")//四則演算
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")//ループ
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")//浮動小数点
const long long mod=998244353LL;
const long long mods=1000000007LL;
const double inf=numeric_limits<double>::infinity();
const int kaz=1000000005;
const long long yab=2500000000000000000LL;
const long long aho =-yab;
const long double eps=1.0e-14L;
const long double pi=acosl(-1.0L);
using ll=long long;
using st=string;
using P=pair<ll,ll>;
using tup=tuple<ll,ll,ll>;
using vi=vector<ll>;
using vin=vector<int>;
using vc=vector<char>;
using vb=vector<bool>;
using vd=vector<double>;
using vs=vector<string>;
using vp=vector<P>;
using sp=set<P>;
using si=set<ll>;
using vvi=vector<vector<ll>>;
using vvin=vector<vin>;
using vvc=vector<vc>;
using vvb=vector<vb>;
using vvvi=vector<vvi>;
using vvvin=vector<vvin>;
const int dx[4]={0,1,0,-1};
const int dy[4]={1,0,-1,0};
const vector<int> ex = {-1, -1, -1, 0, 0, 1, 1, 1};
const vector<int> ey = {-1, 0, 1, -1, 1, -1, 0, 1};
template<typename T>istream&operator>>(istream&is,vector<T>&v){for(T&in:v)is>>in;return is;}
template<typename T>ostream&operator<<(ostream&os,vector<T>v){rep(i,v.size())os<<v[i]<<(i+1!=v.size()?" ":"\n");return os;}
template<typename T1,typename T2>
void co(bool x,T1 y,T2 z){
  if(x)cout << y << endl;
  else cout << z << endl;
}
template<typename T>
bool chmax(T &a, T b){
	if(a<b){
		a=b;
		return true;
	}
	return false;
}
template <typename T, typename U>
T ceil(T x, U y) {
  return (x > 0 ? (x + y - 1) / y : x / y);
}

template <typename T, typename U>
T floor(T x, U y) {
  return (x > 0 ? x / y : (x - y + 1) / y);
}
template<typename T>
bool chmin(T &a, T b){
	if(a>b){
		a=b;
		return true;
	}
	return false;
}
template<typename T>
void print(vector<T> &a){
	for(int i=0;i<a.size();i++){
		cout << a[i];
		if(i==(long long)a.size()-1)cout << endl;
		else cout << " ";
	}
}
template<typename T>
void her(vector<T> &a){
  for(auto &g:a)g--;
}
vector<long long> fac,finv,invs;
// テーブルを作る前処理
void COMinit(int MAX,int MOD) {
	fac.resize(MAX);
	finv.resize(MAX);
	invs.resize(MAX);
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    invs[1] = 1;
    for (int i = 2; i < MAX; i++){
        fac[i] = fac[i - 1] * i % MOD;
        invs[i] = MOD - invs[MOD%i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * invs[i] % MOD;
    }
}

// 二項係数計算
long long COM(int n, int k,int MOD){
    if (n < k) return 0;
    if (n < 0 || k < 0) return 0;
    return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}
ll mypow(ll x,ll y){
	ll ret=1;
	while(y>0){
		if(y&1)ret=ret*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return ret;
}
ll nopow(ll x,ll y){
	ll ret=1;
	while(y>0){
		if(y&1)ret*=x;
		x*=x;
		y>>=1;
	}
	return ret;
}
#include <atcoder/all>
using namespace atcoder;
using minta=modint998244353;
using mintb=modint1000000007;
using mint=mintb;
ostream&operator<<(ostream&os,vector<mintb> v){rep(i,v.size())os<<v[i].val()<<(i+1!=v.size()?" ":"\n");return os;}
//高速ゼータ約数変換
template<typename T>
vector<T> zeta(const vector<T> &x,const bool inv=false,const bool subset=false){
    const size_t n=x.size();
    vector<uint_fast8_t> sieve(x.size(),true);
    vector<T> a=x;
    /*if(inv){
        for(size_t i=1;i<n;i++){
            a[i]-=a[0];
        }
    }*/
    for(size_t i=2;i<n;i++){
        if(!sieve[i])continue;
        if(!subset){
            if(inv){
                for(size_t j=1;i*j<n;j++){
                    sieve[i*j]=false;
                    a[j]-=a[i*j];
                }
            } else{
                for(size_t j=(n-1)/i;j > 0;j--){
                    sieve[i*j]=false;
                    a[j]+=a[i*j];
                }
            }
        } else{
            if(inv){
                for(size_t j=(n-1)/i;j > 0;j--){
                    sieve[i*j]=false;
                    a[i*j]-=a[j];
                }
            } else{
                for(size_t j=1;i*j<n;j++){
                    sieve[i*j]=false;
                    a[i*j]+=a[j];
                }
            }
        }
    }
    /*if(!inv){
        for(size_t i=1;i<n;i++){
            a[i]+=a[0];
        }
    }*/
    return a;
}
int main(){
	cin.tie(nullptr);
	int h,w;
	cin >> h >> w;
	mint mh=h;
	mint mw=w;
	mint res=mh*mw*2-mh-mw;
	int n=max(h,w);
	vector<mint> f(n+1);
	vector<mint> g(n+1);
	reps(i,h)f[i]=h-i;
	reps(j,w)g[j]=w-j;
	f=zeta(f);
	g=zeta(g);
	rep(i,n)f[i]*=g[i];
	f=zeta(f,true);
	res+=f[1]*2;
	cout << res.val() << endl;
}
0