結果

問題 No.942 プレゼント配り
ユーザー tarattata1tarattata1
提出日時 2019-12-05 01:22:51
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 4,484 bytes
コンパイル時間 955 ms
コンパイル使用メモリ 86,456 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-21 11:31:09
合計ジャッジ時間 3,423 ms
ジャッジサーバーID
(参考情報)
judge12 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 20 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,384 KB
testcase_04 AC 20 ms
4,376 KB
testcase_05 AC 20 ms
4,380 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,384 KB
testcase_08 AC 15 ms
4,380 KB
testcase_09 AC 17 ms
4,376 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 1 ms
4,376 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 16 ms
4,380 KB
testcase_15 AC 11 ms
4,380 KB
testcase_16 AC 13 ms
4,380 KB
testcase_17 AC 13 ms
4,380 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 1 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

/**********************************************************
Yukicoder No.940 ワープ ε=ε=ε=ε=ε=│;p>д<│
https://yukicoder.me/problems/no/940

<解法>
n=X+Y+Z とおく。

nより十分に大きいNをとり、
T=(1+x+x^2+..+x^(N-1))(1+y+y^2+..+y^(N-1))(1+z+z^2+..+z^(N-1))-1
 =(1-x^N)/(1-x) (1-y^N)/(1-y) (1-z^N)/(1-z) -1
とおくと、
F=Σ(m=1..n)T^m の x^X y^Y z^Z の係数を求めればよい。

※理由:x^i y^j z^k というのが、1回の移動でx方向にi、y方向にj、z方向にk 進むことを表していて、
  m回移動したときに(X,Y,Z)に到達する場合の数が T^m の x^X y^Y z^Z の係数と等しいので。

F=Σ(m=1..n)T^m の x^X y^Y z^Z の係数を求めるときには
FをxでX回、yでY回、zでZ回偏微分してx=y=z=0 を入れた値を X!Y!Z! で割ればよい。
この計算を行うときには、Nを十分大きくとっておけば、
T=(1-x^N)/(1-x) (1-y^N)/(1-y) (1-z^N)/(1-z) -1
のかわりに
T=1/(1-x)(1-y)(1-z) -1
として計算してもよいことがわかる。

S=(1-x)^(-1) (1-y)^(-1) (1-z)^(-1) とおくと、
F=Σ(m=1..n)(S-1)^m
 =(S-1)(1-(S-1)^n)/(1-(S-1))
これは、(S-1)^nを二項定理で展開した後、多項式の掛け算、割り算をすれば、Σ(m=1..∞)Am S^m の形に書ける。
(係数Am を m<=n の範囲のみ計算することにすれば O(n)で求まる)

あとは、各 S^m についてx^X y^Y z^Z の係数を求めるには
S^m=(1-x)^(-m) (1-y)^(-m) (1-z)^(-m)
を xでX回、yでY回、zでZ回偏微分してx=y=z=0 を入れればよいが、
これは m(m+1)..(m+X-1) と m(m+1)..(m+Y-1) と m(m+1)..(m+Z-1) の積になる。

**********************************************************/

#include <stdio.h>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <list>
#include <iterator>
#include <assert.h>
#pragma warning(disable:4996) 
 
typedef long long ll;
#define MIN(a, b) ((a)>(b)? (b): (a))
#define MAX(a, b) ((a)<(b)? (b): (a))
#define LINF 9223300000000000000
#define INF 2140000000
const long long MOD = 1000000007;
//const long long MOD = 998244353;
using namespace std;


int main(int argc, char* argv[])
{
    int n,K;
    scanf("%d%d", &n, &K);

    int m=n/K;
    if(K%2==0 && m%2) {
        printf("No\n"); return 0;
    }
    else if(K%2==0 || m%2==0) {
        printf("Yes\n");
        vector<vector<int> > z(K);
        int i,j;
        int cnt=0;
        for(i=0; i<m; i++) {
            for(j=0; j<K; j++) {
                int j2=(i%2==0? j: K-1-j);
                z[j2].push_back(cnt+1);
                cnt++;
            }
        }
        for(i=0; i<K; i++) {
            for(j=0; j<m; j++) {
                printf("%d ", z[i][j]);
            }
            printf("\n");
        }
    }
    else {
        if(m<=2) {
            printf("No\n"); return 0;
        }
        else {
            printf("Yes\n");
            vector<vector<int> > z(K);
            int i,j;
            int cnt=0;
            for(i=0; i<m; i++) {
                if(i<m-2) {
                    for(j=0; j<K; j++) {
                        int j2=(i%2==0? j: K-1-j);
                        z[j2].push_back(cnt+1);
                        cnt++;
                    }
                }
                else if(i==m-2) {
                    int K0=(K-1)/2;
                    for(j=0; j<K; j++) {
                        int j1;
                        if(j<K0) j1=K0-1-j;
                        else j1=K0*3-j;
                        int j2=(i%2==0? j1: K-1-j1);
                        z[j2].push_back(cnt+1);
                        cnt++;
                    }
                }
                else {
                    vector<pair<int,int> > v;
                    int p;
                    for(p=0; p<K; p++) {
                        int tmp=z[p][m-3]+z[p][m-2];
                        v.push_back(make_pair(tmp,p));
                    }
                    sort(v.rbegin(), v.rend());
                    for(j=0; j<K; j++) {
                        z[v[j].second].push_back(cnt+1);
                        cnt++;
                    }
                }
            }
            for(i=0; i<K; i++) {
                for(j=0; j<m; j++) {
                    printf("%d ", z[i][j]);
                }
                printf("\n");
            }
        }
    }


    return 0;
}
0