F. 组合数问题 (组合数推公式&逆元) (2021牛客寒假算法基础集训营6)

最新推荐文章于 2022-03-25 20:25:59 发布

S atur

最新推荐文章于 2022-03-25 20:25:59 发布

阅读量200

点赞数

分类专栏：数学妙啊文章标签：算法

本文链接：https://blog.csdn.net/Satur9/article/details/114444075

版权

数学同时被 2 个专栏收录

68 篇文章 0 订阅

订阅专栏

妙啊

16 篇文章 0 订阅

订阅专栏

博客介绍了如何使用数学上的组合数公式来求解模意义下的幂次问题。通过公式推导和代码实现，展示了计算特定指数的快速幂算法，特别是在n为4的倍数或奇数时的情况。最后，给出了C++代码实现，用于计算给定整数n下幂次的模反元素。

摘要由CSDN通过智能技术生成

传送门

思路：学习大佬博客1、大佬博客2

根据数学上的组合数公式：

$(1+i)^{n} = C_{n}^{0} +iC_{n}^{1} - C_{n}^{2} - iC_{n}^{3} + C_{n}^{4}......$ ...................①

$(1-i)^{n} = C_{n}^{0} - iC_{n}^{1} - C_{n}^{2} + iC_{n}^{3} + C_{n}^{4}......$ ...................②

由①和②相加除以二可得到公式③：

$C_{n}^{0} - C_{n}^{2} + C_{n}^{4}......$ ............................................................③

又因为 $(1+i)^{4} = (1-i)^{4} = -4$ ，所以 ③ = $(-4)^{n/4}$ .

又因为公式④：

$C_{n}^{0} + C_{n}^{2} + C_{n}^{4} ...... = 2^{n-1}$ ...............................................④

所以③和④相加除以二便得到了答案。

代码实现：

#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ll long long
#define int long long
#define pii pair<int, int>
#define lowbit(x) (x &(-x))
#define ls(x) x<<1
#define rs(x) (x<<1+1)
#define me(ar) memset(ar, 0, sizeof ar)
#define mem(ar,num) memset(ar, num, sizeof ar)
#define rp(i, n) for(int i = 0, i < n; i ++)
#define rep(i, a, n) for(int i = a; i <= n; i ++)
#define pre(i, n, a) for(int i = n; i >= a; i --)
#define IOS ios::sync_with_stdio(0); cin.tie(0);cout.tie(0);
const int way[4][2] = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}};
using namespace std;
const int  inf = 0x3f3f3f3f;
const double PI = acos(-1.0);
const double eps = 1e-6;
const ll   mod = 998244353;
const int  N = 2e5 + 5;

inline void read(int &x){
    char t=getchar();
    while(!isdigit(t)) t=getchar();
    for(x=t^48,t=getchar();isdigit(t);t=getchar()) x=x*10+(t^48);
}

int qmi(int a, int k){
    int res = 1;
    while(k){
        if(k & 1) res = (ll)res * a % mod;
        k >>= 1;
        a = (ll)a * a % mod;
    }
    return res;
}

signed main()
{
    IOS;

    int n; cin >> n;
    int ans = qmi(4, n/4);
    if((n/4)&1) ans = -ans+mod;
    int res = qmi(2, n-1);
    ans = (ans+res)%mod*qmi(2, mod-2)%mod;
    cout << ans << endl;

    return 0;
}