結果
| 問題 |
No.2866 yuusaan's Knapsack
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-10-12 12:34:10 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 3,431 bytes |
| コンパイル時間 | 1,364 ms |
| コンパイル使用メモリ | 125,920 KB |
| 最終ジャッジ日時 | 2025-02-24 17:28:00 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 4 RE * 22 |
ソースコード
#include<iostream>
#include<string>
#include<queue>
#include<vector>
#include<cassert>
#include<random>
#include<set>
#include<map>
#include<cassert>
#include<unordered_map>
#include<bitset>
#include<numeric>
#include<algorithm>
using namespace std;
typedef long long ll;
const int inf=1<<30;
const ll INF=1LL<<62;
using P=pair<ll,ll>;
typedef pair<int,P> PP;
const ll MOD=998244353;
const int mod = 998244353;
struct mint {
long long x; // typedef long long long long;
mint(long long x=0):x((x%mod+mod)%mod){}
mint operator-() const { return mint(-x);}
mint& operator+=(const mint a) {
if ((x += a.x) >= mod) x -= mod;
return *this;
}
mint& operator-=(const mint a) {
if ((x += mod-a.x) >= mod) x -= mod;
return *this;
}
mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;}
mint operator+(const mint a) const { return mint(*this) += a;}//後ろのconstはメンバxを変更しないことを示す
mint operator-(const mint a) const { return mint(*this) -= a;}
mint operator*(const mint a) const { return mint(*this) *= a;}
bool operator==(const mint a) const {return a.x==x;}
mint pow(unsigned long long int t) const {
assert(t>=0);
if (!t) return 1;
//aがtで割り切れない場合には t%=(mod-1)をして良い
mint a = pow(t>>1);
a *= a;
if (t&1) a *= *this;
return a;
}
// for prime mod
mint inv() const { return pow(mod-2);}
mint& operator/=(const mint a) { return *this *= a.inv();}
mint operator/(const mint a) const { return mint(*this) /= a;}
};
std::istream& operator>>(std::istream& is, mint& a) { return is >> a.x;}
std::ostream& operator<<(std::ostream& os, const mint& a) { return os << a.x;}
struct combination {
std::vector<mint> fact, ifact;
combination(int n):fact(n+1),ifact(n+1) {
assert(n < mod);
fact[0] = 1;
for (int i = 1; i <= n; ++i) fact[i] = fact[i-1]*i;
ifact[n] = fact[n].inv();
for (int i = n; i >= 1; --i) ifact[i-1] = ifact[i]*i;
}
mint operator()(int n, int k) {
if (k < 0 || k > n) return 0;
return fact[n]*ifact[k]*ifact[n-k];
}
};
int main(){
int N,W;
cin>>N>>W;
vector<ll> v(N),w(N);
vector<P> vw(N);
for(int i=0;i<N;i++){
cin>>vw[i].first>>vw[i].second;
}
sort(vw.begin(),vw.end(),
[&](P& a,P& b){
//重さの昇順
return a.second<b.second;
});
using X=pair<ll,mint>;
int OFFSET=10000;
const int MAXW = 2*10005 + OFFSET;//負の重さ
vector<vector<X>> dp(N+3,vector<X>(MAXW+5,X{-INF,0}));
dp[0][OFFSET+0]=X{0,1};
auto f=[&](X x,X y){
if(x.first!=y.first){
if(x.first<y.first){
swap(x,y);
}
//x.first>y.first
return x;
}
return X{x.first,x.second+y.second};
};
for(int i=0;i<N;i++){
auto [v,w]=vw[i];
for(int j=0;j<=MAXW;j++){
dp[i+1][j]=dp[i][j];
}
for(int j=0;j<=MAXW;j++){
//{価値の合計のmax,その価値の組み合わせの数}
if(j-w>=0 && j-w<=2*MAXW && dp[i][j-w].first!=-INF) {
dp[i+1][j]=f(dp[i+1][j],X{dp[i][j-w].first+v,dp[i][j-w].second});
}
}
}
X ans{-INF,0};
for(int j=0;j<=W+OFFSET;j++){
ans=f(ans,dp[N][j]);
}
cout<<ans.first<<' '<<ans.second<<endl;
}