結果

問題 No.1825 Except One
ユーザー srjywrdnprkt
提出日時 2025-03-29 02:42:21
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 80 ms / 3,000 ms
コード長 1,283 bytes
コンパイル時間 3,437 ms
コンパイル使用メモリ 283,024 KB
実行使用メモリ 7,328 KB
最終ジャッジ日時 2025-03-29 02:42:27
合計ジャッジ時間 5,643 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
//#include <atcoder/modint>

using namespace std;
//using namespace atcoder;
using ll = long long;
//using mint = modint998244353;

int main(){
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);

    /*
       S回全体からインクリメントしてS回1を足す。
       S >= max(A)でないといけない。
       SN - sum(A) = S
       S(N-1) = sum(A) and S >= max(A)
       後者の条件はAをソートすれば良い。
       dp(i, j, k) = i番目まででj個選んで、総和がkであるような状態数
       O(N^3A)
    */
    ll N, M, ans=0;
    cin >> N;
    M = N*100;
    vector dp(N+1, vector<ll>(M+1));
    vector<ll> A(N);
    for (int i=0; i<N; i++) cin >> A[i];
    sort(A.begin(), A.end());

    dp[0][0] = 1;
    for (auto a : A){
        vector pd(N+1, vector<ll>(M+1));
        for (int j=0; j<=N; j++){
            for (int k=0; k<=M; k++){
                pd[j][k] += dp[j][k];
                if (j && k-a >= 0){
                    pd[j][k] += dp[j-1][k-a];
                    if (j >= 2 && k % (j-1) == 0 && k / (j-1) >= a){
                        ans += dp[j-1][k-a];
                    }
                }
            }
        }
        swap(dp, pd);
    }
    cout << ans << endl;

    return 0;
}
0