Link to this code:
https://cses.fi/paste/f6b687511890d6331a4ea4/
// ░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░
// ░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░
// ░░░░░░░░░░░░░░░░░░░░░░░░░█░░░░░░░░░░░░░░
// ░░░░░░░░░░░░░░░░░░░░░░░███░░░░░░░░░░░░░░
// ░░░░░░░░░░░██░░░░░░░░██░░░░░░░░░░░░░░░░░
// ░░░░░░░░░░██░░░░░░███░░░░░░░░░░░░░░░░░░░
// ░░░░░░░░███░░░░░░███░░░░░░░░░░░░░░░░░░░░
// ░░░░░░░██░░░░░░░██░██░░░░░░░░░░░░░░░░░░░
// ░░░░░░█░░░░░░░░██░░░██░░░░░░░██░░░░░░░░░
// ░░░░░████░░░░░░█░░░░░██░░░░███░░░░░░░░░░
// ░░░░██░░███░░░██░░░░░░██░███░░░░░░░░░░░░
// ░░░██░░░░░█████░░░░░░░░███░░░░░░░░░░░░░░
// ░░██░░░░░░░░██████████████░░░░░░░░███░░░
// ░██░░░░░░░░░█░░░░░░░░░░░░█░░░░░░░██░░░░░
// ░█░░░░░░░░░██░░░░░░░░░░░░░█░░░░░██░░░░░░
// ░░░░░░░░░░█░░░█▀█▌░▀▀█░█░▌██░░███░░░░░░░
// ░░░░░░░░░██░░░█▄▄█░███░██▌░████░░░░░░░░░
// ░░░░░░░░░█░░░░█▄▄█░▌░░░░░▌░░█░░░░░░░░░░░
// ░░░░░░░░░░░░░░░░░░░▀▀▀░░░░░░░░░░░░░░░░░░
#include<bits/stdc++.h>
using namespace std;
#define int long long int
#define ll long long int
#define fi(n) for(int i=0;i<n;i++)
#define fj(n) for(int j=0;j<n;j++)
#define W(x) int x;cin>>x;while(x--)
#define fr(n) for(int i=n-1;i>=0;i--)
#define mp make_pair
#define fo(i,a,n) for(int i=a;i<n;i++)
#define ps(x,y) fixed<<setprecision(y)<<x
#define pb push_back
#define vi vector<int>
#define deb(x) cout<<#x<<"="<<x<<endl;
#define all(a) (a).begin(),(a).end()
#define pii pair<int,int>
#define sortall(x) sort(all(x))
#define mk(arr,n,type) type *arr=new type[n];
#define vpi vector<pair<int,int> >
#define setbits(x) __builtin_popcountll(x)
#define zerbits(x) __builtinctzll(x)
#define mod 1000000007
#define inf 1e18
#define sz(x) (int)x.size()
#define vs vector<string>
#define itr(a) for(auto it=a.begin();it!=a.end();it++)
#define f first
#define s second
#define INT_SIZE 32
int mpow(int exp, int base);
//I don't know what the funciton aeh() does but it speeds up my mess
void aeh()
{
ios_base::sync_with_stdio(false);
cin.tie(NULL); cout.tie(0);
#ifndef ONLINE_JUDGE
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
#endif
}
// int dp[1000][100000];
// int fn(int n,int x,int cpos)
// {
// if(cpos<0)
// return 0;
// if(dp[cpos][x]!=-1)
// return dp[cpos][x];
// if(h[cpos]>x)
// return dp[x][cpos]=fn(n,x,cpos-1);
// return dp[cpos][x]=max(s[cpos]+fn(n,x-h[cpos],cpos-1),fn(n,x,cpos-1));;
// }
int32_t main()
{
aeh();
int n,x;
cin>>n>>x;
vi h(n+1);
vi s(x+1);
fi(n)
cin>>h[i];
fi(n)
cin>>s[i];
int dp[n+1][x+1];
fi(n+1)
{
fj(x+1)
{
if(i==0||j==0)
dp[i][j]=0;
else if(h[i-1]<=j){
dp[i][j]=max(s[i-1]+dp[i-1][j-h[i-1]],dp[i-1][j]);
}
else
dp[i][j]=dp[i-1][j];
}
}
cout<<dp[n][x];
}
int mpow(int base, int exp)
{
base %= mod;
int result = 1;
while (exp > 0)
{
if (exp & 1)result = (result * base) % mod;
base = (base * base) % mod;
exp >>= 1;
}
return result;
}