[분치 결정 단순성 최적화 DP] LOJ#6039.'아례 합숙 2017 Day 5'보석

fi,jfi,j를 설정하면 판매가≤i≤i의 보석을 표시하고 jj만 위안의 최대 흡인력을 사용합니다.전이 방정식은 쓰기 쉽다.같은 가격의 보석은 틀림없이 가장 큰 흡인력에서 시작된 것이기 때문에 이동하는 함수는 상철이다.jj를 jmod i jm o d i로 분류하고 각 유형의 이동은 의사결정의 단조성을 만족시킨 다음에 분치로 할 수 있다.시간 복잡도 O(maxciklogk) O(m a x c i k log ⁡k)는 실제 복잡도가 이 값보다 훨씬 작다.

#include
using namespace std;
char nc() {
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
void Read(int& x) {
    char c=nc();
    for(;c<'0'||c>'9';c=nc());
    for(x=0;c>='0'&&c<='9';x=(x<<3)+(x<<1)+c-48,c=nc());
}
typedef long long ll;
const int N=50010;
int k,n,m,mx,x,y,num;
ll f[N],h[N];
vector<int>g[310];
ll s[1000010];
bool b;
bool Cmp(ll x,ll y) {
    return x>y;
}
void Solve(int x,int y,int l,int r,int L,int R) {
    if(l>r) return;
    int Mid=l+r>>1,pos=L;ll mx=0;
    for(int i=max(L,Mid-num);i<=R&&i<=Mid;i++)
        if(f[(i-1)*x+y]+s[Mid-i]>=mx) mx=f[(i-1)*x+y]+s[Mid-i],pos=i;
    h[(Mid-1)*x+y]=mx;
    Solve(x,y,l,Mid-1,L,pos);Solve(x,y,Mid+1,r,pos,R);
}
int main() {
    Read(n);Read(m);
    for(int i=1;i<=n;i++) Read(x),Read(y),g[x].push_back(y),mx=max(mx,x);
    for(int i=1;i<=mx;i++) 
        if(g[i].size()) {
            sort(g[i].begin(),g[i].end(),Cmp);
            num=g[i].size();
            for(int j=0;j1]=s[j]+g[i][j];
            for(int j=0;jint cnt=0;
                for(int k=j;k<=m;k+=i) cnt++;
                Solve(i,j,1,cnt,1,cnt);
                for(int k=j;k<=m;k+=i) f[k]=h[k];
            }
        }
    for(int i=1;i<=m;i++) printf("%lld ",f[i]);
    return 0;
}