题解 AT_nikkei2019ex_e【コラッツ問題】

发布时间 2023-05-30 22:41:21作者: rui_er

啥玩意,诈骗题还能这么诈骗。

\(f(X)\) 就是角谷猜想(冰雹猜想)所需的步数。根据角谷猜想,定义函数 \(g\)

\[g(X)= \begin{cases} \frac{X}{2},&2\mid X\\ 3X+1,&2\nmid X \end{cases} \]

则显然有 \(f(g(X))=f(X)-1\)。样例二已经给出了 \(f(X)=1000\) 的解 \(X=1789997546303\),而 \(P\) 的范围是 \([0,1000]\),根据上面的结论递推即可求出 \(f(X)=P\) 的解。

//By: OIer rui_er
#include <bits/stdc++.h>
#define rep(x,y,z) for(ll x=(y);x<=(z);x++)
#define per(x,y,z) for(ll x=(y);x>=(z);x--)
#define debug(format...) fprintf(stderr, format)
#define fileIO(s) do{freopen(s".in","r",stdin);freopen(s".out","w",stdout);}while(false)
using namespace std;
typedef long long ll;

mt19937 rnd(std::chrono::duration_cast<std::chrono::nanoseconds>(std::chrono::system_clock::now().time_since_epoch()).count());
ll randint(ll L, ll R) {
    uniform_int_distribution<ll> dist(L, R);
    return dist(rnd);
}

template<typename T> void chkmin(T& x, T y) {if(x > y) x = y;}
template<typename T> void chkmax(T& x, T y) {if(x < y) x = y;}

const ll N = 1e3+5;

ll n, ans[N];

int main() {
    ans[1000] = 1789997546303LL;
    per(i, 999, 0) {
        if(ans[i+1] & 1) ans[i] = 3 * ans[i+1] + 1;
        else ans[i] = ans[i+1] / 2;
    }
    scanf("%lld", &n);
    printf("%lld\n", ans[n]);
    return 0;
}