Lưu ý: Các code mẫu dưới đây chỉ mang tính tham khảo và có thể không AC được bài tập này

Code mẫu của skyvn97

#include<bits/stdc++.h>
#define MAXK   555
#define MAXM   111
#define FOR(i,a,b) for (int i=(a),_b=(b);i<=_b;i=i+1)
#define REP(i,n) for (int i=0,_n=(n);i<_n;i=i+1)
using namespace std;
typedef long long ll;
ll f[MAXM][MAXK];
int cnt[MAXM];
int n,m,q;
template<class T>
    inline void minimize(T &x,const T &y) {
        if (x>y) x=y;
    }
void init(void) {
    scanf("%d%d%d",&n,&m,&q);
    REP(zz,n) {
        int t;
        scanf("%d",&t);
        cnt[t]++;
    }
}
ll sum(int n) {
    return (1LL*n*(n+1)/2);
}
ll divide(int n,int k) {
    if (n<=k) return (n);
    return (sum(n/k)*(k-n%k)+sum(n/k+1)*(n%k));
}
void optimize(void) {
    memset(f,0x3f,sizeof f);
    REP(j,q+1) f[0][j]=0;
    FOR(i,1,m) FOR(j,0,q) FOR(k,0,j)
        minimize(f[i][j],f[i-1][j-k]+divide(cnt[i],k+1));
    cout<<f[m][q];
}
int main(void) {
    init();
    optimize();
    return 0;
}