ARC106-D.数学,二项式定理

题目大意:

很简洁,直接放题面吧.
在这里插入图片描述

题目思路:

很容易看出要使用二项式定理。但是需要一点推公式的能力。

回顾一下二项式定理: ( a + b ) x = ∑ i = 0 x C x i ∗ a i ∗ b x − i (a+b)^x=\sum_{i=0}^{x}C_x^i*a^i*b^{x-i} (a+b)x=i=0xCxiaibxi

先将式子调整一下:
原式= 1 2 [ ∑ i = 1 N ∑ j = 1 N ( A i + A j ) x − ∑ i = 1 N ( 2 A i ) x ] \frac{1}{2}[\sum_{i=1}^{N}\sum_{j=1}^{N}(A_i+A_j)^x-\sum_{i=1}^{N}(2A_i)^x] 21[i=1Nj=1N(Ai+Aj)xi=1N(2Ai)x]

后面部分可以O(n)直接算。

关注前半部分: 原式 = ∑ i = 1 N ∑ j = 1 N ∑ k = 0 x C x k ∗ ( A i ) k ∗ ( A j ) x − k \sum_{i=1}^{N}\sum_{j=1}^{N}\sum_{k=0}^{x}C_x^k*(A_i)^k*(A_j)^{x-k} i=1Nj=1Nk=0xCxk(Ai)k(Aj)xk

= ∑ k = 0 x C x k [ ∑ i = 1 N ( A i ) k ] [ ∑ j = 1 N ( A j ) x − k ] =\sum_{k=0}^{x}C_x^k[\sum_{i=1}^{N}(A_i)^k][\sum_{j=1}^{N}(A_j)^{x-k}] =k=0xCxk[i=1N(Ai)k][j=1N(Aj)xk]

显然,上述式子可以独立算了。

所以 O ( n k ) O(nk) O(nk)预处理 v a l [ j ] = ∑ i = 1 N ( A i ) j val[j]=\sum_{i=1}^{N}(A_i)^j val[j]=i=1N(Ai)j,再 O ( k 2 ) O(k^2) O(k2)预处理组合数.

然后枚举 k k k, O ( k ) O(k) O(k)的计算 ∑ k = 0 x C x k \sum_{k=0}^{x}C_x^k k=0xCxk.然后三者相乘即可.最后记得乘2的逆元.

总时间复杂度: O ( n k ) O(nk) O(nk)

AC代码:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
#define pb push_back
#define mp make_pair
const int maxn = 2e5 + 5;
const int mod = 998244353;
ll a[maxn] , b[maxn] , val[305];
ll c[305][305];
ll ksm (ll a , ll b)
{
    ll ans = 1 , base = a;
    while (b){
        if (b & 1) ans = ans * base % mod;
        b >>= 1;
        base = base * base % mod;
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    for (int i = 1 ; i <= 300 ; i++){
        c[i][0] = c[i][i] = 1;
        for (int j = 1 ; j < i ; j++)
            c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
    }
    ll n , k; cin >> n >> k;
    for (int i = 1 ; i <= n ; i++)
        cin >> a[i] , b[i] = 1;
    val[0] = n;
    for (int j = 1 ; j <= k ; j++){
        for (int i = 1 ; i <= n ; i++){
            b[i] = (b[i] * a[i]) % mod;
            val[j] = (val[j] + b[i]) % mod;
        }
    }
    ll base = 1 , inv = ksm(2 , mod - 2);
    for (int i = 1 ; i <= k ; i++){
        ll res =  0;
        for (int j = 0 ; j <= i ; j++)
            res = (res + c[i][j] * val[j] % mod * val[i - j] % mod) % mod;
        base = base * 2 % mod;
        res = (res - base * val[i] % mod + mod) % mod;
        res = res * inv % mod;
        cout << res << endl;
    }
    return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
碰撞检测是指在游戏或者模拟中检测两个物体是否发生了碰撞。在网页开发中,我们可以使用JavaScript实现简单的碰撞检测。在01-object-hit-test.html中,实现了一个基于Canvas的简单碰撞检测示例,通过判断两个圆形物体之间的距离是否小于它们的半径之和来判断它们是否发生了碰撞。 具体实现步骤如下: 1. 定义两个圆形物体(对象),包括其位置、半径、颜色等属性。 ``` var circle1 = { x: 100, y: 100, r: 50, color: "#FF0000" }; var circle2 = { x: 200, y: 200, r: 80, color: "#0000FF" }; ``` 2. 在Canvas中绘制两个圆形物体。 ``` ctx.beginPath(); ctx.arc(circle1.x, circle1.y, circle1.r, 0, 2 * Math.PI); ctx.fillStyle = circle1.color; ctx.fill(); ctx.beginPath(); ctx.arc(circle2.x, circle2.y, circle2.r, 0, 2 * Math.PI); ctx.fillStyle = circle2.color; ctx.fill(); ``` 3. 计算两个圆形物体之间的距离。 ``` var distance = Math.sqrt(Math.pow(circle1.x - circle2.x, 2) + Math.pow(circle1.y - circle2.y, 2)); ``` 4. 判断两个圆形物体是否发生碰撞。 ``` if (distance < circle1.r + circle2.r) { console.log("Collision detected!"); } ``` 完整代码如下: ``` <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title>Object Hit Test</title> </head> <body> <canvas id="canvas" width="400" height="400"></canvas> <script> var canvas = document.getElementById("canvas"); var ctx = canvas.getContext("2d"); var circle1 = { x: 100, y: 100, r: 50, color: "#FF0000" }; var circle2 = { x: 200, y: 200, r: 80, color: "#0000FF" }; ctx.beginPath(); ctx.arc(circle1.x, circle1.y, circle1.r, 0, 2 * Math.PI); ctx.fillStyle = circle1.color; ctx.fill(); ctx.beginPath(); ctx.arc(circle2.x, circle2.y, circle2.r, 0, 2 * Math.PI); ctx.fillStyle = circle2.color; ctx.fill(); var distance = Math.sqrt(Math.pow(circle1.x - circle2.x, 2) + Math.pow(circle1.y - circle2.y, 2)); if (distance < circle1.r + circle2.r) { console.log("Collision detected!"); } </script> </body> </html> ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值