結果

問題 No.2167 Fibonacci Knapsack
コンテスト
ユーザー C
提出日時 2026-11-05 14:56:41
言語 C++23(gcc16)
(gcc 16.2.0 + boost 1.92.0 + ACL)
コンパイル:
g++-16 -O2 -fplugin=/usr/lib/yukicoder/noavx512-gcc16.so -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 2 ms / 2,000 ms
+ 671µs
コード長 2,086 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 7,158 ms
コンパイル使用メモリ 393,424 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-11-05 14:56:51
合計ジャッジ時間 9,573 ms
ジャッジサーバーID
（参考情報） 		judge4_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#include <atcoder/all>
#include <boost/polygon/point_data.hpp>
#include <boost/polygon/voronoi.hpp>
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
#define all(x) begin(x),end(x)
#define sz(x) (int)(x).size()
using namespace std;
using namespace atcoder;
using ll=long long;
using ull=unsigned long long;
using ld=long double;
using vll=vector<ll>;
using vvll=vector<vll>;
using vi=vector<int>;
using vvi=vector<vi>;
using pll=pair<ll,ll>;
using pi=pair<int,int>;

// using mint=modint998244353;
using mint=modint1000000007;
// using mint2=static_modint<2>;

// using BoostPoint=boost::polygon::point_data<int>;
// using Voronoi=boost::polygon::voronoi_diagram<double>;

// using mint=LazyMontgomeryModInt<998244353>;
using vm=vector<mint>;
// using vvm=vector<vm>;

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
	/**
	 * 要はF_i,F_{i+1},...,F_N任意の部分和a<b<cについて差分が大きいです：
	 * 
	 * b - a >= F_{i-1}  (1)
	 * 
	 * c - a >= F_{i+1}　（2）
	 * 
	 * これで次のように枝刈りできます：
	 * 
	 * 品物i,i+1,...,Nのみ使える時、価値の最大値をB、任意の部分和をvとして
	 * 
	 * 明らかに　v + \sum_{j=1}^{i-1} F_j < B　のような v は捨てれば良い
	 * 
	 * 残された v は必ず　B - F_{i+1} + 2 <= v <= B ので　（ここは\sum_{j=1}^{i-1} F_j　= F_{i+1} - 2を使います）
	 * 
	 * （2）よりこのような v は高々２個しか存在しない
	 */
	int T,N;
	cin>>T;
	ll F[76]={0,1,2},S[76]={};
	rep(i,3,76)F[i]=F[i-1]+F[i-2];
	rep(i,1,76)S[i]=S[i-1]+F[i];
	while(T--){
		ll W;
		cin>>N>>W;
		vll w(N+1);
		rep(i,1,N+1)cin>>w[i];
		vector<pll>dp={{0,0}};
		rrep(i,1,N+1){
			vector<pll>ndp=dp;
			for(auto&[v,tw]:dp){
				if(tw+w[i]<=W)ndp.push_back({v+F[i],tw+w[i]});
			}
			sort(all(ndp));
			ll B=ndp.back().first;
			dp.clear();
			for(auto&[v,tw]:ndp){
				if(v+S[i-1]<B||(!dp.empty()&&v==dp.back().first))continue;
				dp.push_back({v,tw});
			}
		}
		cout<<dp.back().first<<'\n';
	}
}
0