HDU A simple stone game K倍动态规划减法游戏的一些理解。

两人取一堆n个石子 先手不能全部取完 之后每人取的个数不能超过另一个人上轮取的数*K，给n，K判断先手必胜并求第一步走法。
先看k=2

斐波拉契数列
0 1 1 2 3 5 8 13，
有关这个数列首先明确这几个事实
1.任意一个整数可以写成斐波那契数列中的不相邻的项的和，
2.任何一个整数n都可以写成两项斐波那契数的和
3.
$f(i)=f(i-1)+f(i-2)$
$f(i)=f(i-2)\ast 2+s$; $s=f(i-1)-f(i-2)=f(i-3)$
$f(i)=f(i-2)\ast 2+f(i-3)$ $<3\ast f(i-2)$
这样 一个数分解为两个斐波拉切数，1，b，c，
先取走1，那么，下一个人无法完全取走b，但是第一个人可以根据第二个人的情况凑出一个较小的斐波拉切数取走，那么，下一个人又不能完全取走下一个斐波拉切数，但前提是这个是数本身不是斐波拉切数。
1,3，8，
先取1，然后能凑出3，然后能凑出8，都是后手能凑出。
然后观察波拉切序列的性质，有这样一条，
$f(i)=f(i-2)\ast 2+s$; $s=f(i-1)-f(i-2)=f(i-3)$
$f(i)=f(i-2)\ast 2+f(i-3)$ $<3\ast f(i-2)$
也就是2倍不满，3倍能超，这个序能够展现k=2所体现的一种临界状态。任意一个整数可以写成斐波那契数列中的不相邻的项的和，为什么要搞成不相邻的两项的组合，就是为了这样的状态展现。
在看k=1, 有序列
1 ,2,4,8，32,（二进制序列）
能看出 一倍不满下一个，两倍刚好下一个，完美展现k=1的一种临界，
那么，我们要找一个普适的办法完美展现k=k的临界的序列。
将n写成数列中一些项的和，使得这些被取到的项的相邻两个倍数差距>k，k倍不满，（k+1）倍能凑，设这个数列已经被构造了i 项，第 i 项为a[ i ]，前 i 项可以完美对1…b[ i ] 编码使得每个编码的任意两项倍数>K 那么有a[ i+1 ] = b[ i ] + 1;这是显然的 因为b[ i ] + 1没法构造出来，只能新建一项表示，然后关于b[i]的计算，是a[0…i]构造最大数字，那么 a[i]是一定要选用的，那么要选用的下一项只能递减寻找，直到找到 a[t] 满足 a[t] * K < a[i] ，而b[t]就是a[0…t]所能构造的最大数字，再加上a[i], 即为a[0…i]能构造的最大数字，于是b[i] = b[t] + a[i]。
然后，就是记板子了。

然后有个详细的讲解
神牛曹钦翔的论文：从“k倍动态减法游戏”出发探究一类组合游戏问题
https://wenku.baidu.com/view/79b2b2360b4c2e3f572763a8.html

#include <iostream>
#include <cstdio>
#include <cmath>
#include <bits/stdc++.h>
#define int long long
using namespace std;
//#define int long long
// 这里的read与下面的相比更加灵活。
// Complex 写成complex有可能导致命名冲突。
inline int read()
{
    char c = getchar();
    int x = 0, f = 1;
    while (c < '0' || c > '9')
    {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9')
    {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return x * f;
}
/*template<class cl>void read(cl& x)
{
    x=0;int f=0;char ch;ch=getchar();
    while(!isdigit(ch)){f=f|(ch=='-'),ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
    x=f?-x:x;
    return ;
}*/
template <class cl>
void put(cl x)
{
    if (x < 0)
        putchar('-'), x = -x;
    if (x > 9)
        put(x / 10);
    putchar(x % 10 + '0');
    return;
}
const int N=20000005;
int a[N], b[N];
int init(int n, int k)
{
    int i = 0, j = 0;
    while(n > a[i]){
            i ++;
            a[i] = b[i - 1] + 1;
            while(a[j + 1] * k < a[i])
                j ++;
            if(k * a[j] < a[i])
                b[i] = b[j] + a[i];
            else
				b[i] = a[i];
        }
    return i;
}
int n,k;
int ans;
int cnt=0;
void to()
{
   n=read();k=read();
   int x=init(n,k);
   printf("Case %d: " ,++cnt);
   if(a[x]==n)
   {
      puts("lose");
   }
   else 
   {
      while(n)
      {
          if(n>=a[x])
          {
              n=n-a[x];
              ans=a[x];
          }
          x--;
      }
      put(ans);
      puts("");
   }
}
signed main()
{
   int t;
   t=read();
   while(t--)
   {
       to();
   }
   return 0;
}