結果

問題 No.824 Many Shifts Hard
ユーザー Pulmn
提出日時 2018-09-01 23:55:53
言語 C++11
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,251 bytes
コンパイル時間 489 ms
コンパイル使用メモリ 58,236 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-24 04:58:51
合計ジャッジ時間 1,667 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:12:9: warning: iteration 304 invokes undefined behavior [-Waggressive-loop-optimizations]
   12 |         for(int i=0;i<M;i++) for(int j=0;j<=i;j++){
      |         ^~~
main.cpp:12:22: note: within this loop
   12 |         for(int i=0;i<M;i++) for(int j=0;j<=i;j++){
      |                     ~^~

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
const int M=305;
ll N,K,C[M][M],D[M][M],a[3][M];
int main(){
    cin>>N>>K;
    C[0][0]++;
    for(int i=0;i<M;i++) for(int j=0;j<=i;j++){
        (C[i+1][j]+=C[i][j])%=mod;
        (C[i+1][j+1]+=C[i][j])%=mod;
    }
    for(int i=0;i<3;i++) for(int j=0;j<M;j++){
        if(j) a[i][j]=a[i][j-1]*(N-i)%mod;
        else a[i][j]=1;
    }
    for(int i=1;i<M;i++) for(int j=0;j<M;j++) if(i+j<=K){
        for(int k=i;k<=K-j;k++){
            (D[i][j]+=C[k-1][i-1]*C[K-k][j]%mod*a[2][k-i]%mod*a[1][K-j-k])%=mod;
        }
    }
    cout<<N<<endl;
    for(int i=0;i<3;i++){
        for(int j=0;j<10;j++) cout<<a[i][j]<<' ';
        cout<<endl;
    }
    ll sum=N*(N+1)/2%mod*a[0][K],res=0;
    cout<<sum<<endl;
    for(int i=N;i<=K;i++) (res+=N*(N+1)/2%mod*C[i-1][N-1]*a[1][i-N]%mod*a[0][K-i])%=mod;
    cout<<res<<endl;
    for(int i=1;i<=K;i++) if(i<N) (res+=i*(2*N-i+1)/2%mod*C[K][i]%mod*a[1][K-i])%=mod;
    cout<<res<<endl;
    for(int i=1;i<=K;i++) for(int j=0;j<=K;j++) if(j<i&&i+j<=K&&i<N){
        ll S=(N-i-1)*(i-j)*(N-j+1)/2%mod,t=(C[i+j][i]-C[i+j][i+1]+mod)%mod;
        (res+=C[K][i+j]*a[2][K-i-j]%mod*S%mod*t)%=mod;
        for(int k=0;k<=j;k++){
            S=(i-j)*(i-j+1)/2;
            t=(C[i+k-1][k]-C[i+k-1][i]+mod)%mod;
            (res+=D[i+k][j-k]*S%mod*t)%=mod;
        }
    }
    cout<<(sum-res+mod)%mod<<endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0