「BZOJ3771」Triple-生成函数+FFT

最新推荐文章于 2020-12-08 09:36:09 发布

DSL_HN_2002

最新推荐文章于 2020-12-08 09:36:09 发布

阅读量320

点赞数

分类专栏：多项式——FFT/NTT 多项式——生成函数文章类型——题解

本文链接：https://blog.csdn.net/DSL_HN_2002/article/details/81453828

版权

文章类型——题解同时被 3 个专栏收录

57 篇文章 0 订阅

订阅专栏

多项式——FFT/NTT

11 篇文章 0 订阅

订阅专栏

多项式——生成函数

4 篇文章 0 订阅

订阅专栏

Decription

给定一个物品集合，每个物品有一个价值。可以在其中使用 $1,2,3$ 个物品，求凑成的物品总价值方案数。

$n,a_i \leqslant 40000$

Solution

设 $a_i$ 表示取一个物品价值和为 $i$ 的方案数, $b_i,c_i$ 分别为两个，三个物品的方案数。

对于每一个出现的价值，就在其对应的生成函数上的位置 $+1$ 。

通过容斥不难得到以下式子

$ans=\sum a_i+\frac{(a_i \times a_i-b_i)}{2}+\frac{a_i \times a_i \times a_i - b_i \times a_i \times 3 + c_i \times 2}{6}$

#include <bits/stdc++.h>
using namespace std;

const int maxn = 40005, maxm = 131075; 
int n, maxa; 

inline int gi()
{
    char c = getchar();
    while (c < '0' || c > '9') c = getchar();
    int sum = 0;
    while ('0' <= c && c <= '9') sum = sum * 10 + c - 48, c = getchar();
    return sum;
}

typedef complex<double> cpx;
const double pi = acos(-1);
cpx a[maxm], b[maxm], c[maxm], ans[maxm];

int m, L, R[maxm];
void fft(cpx *a, int f)
{
    for (int i = 0; i < m; ++i)
        if (i < R[i]) swap(a[i], a[R[i]]);
    for (int i = 1; i < m; i <<= 1) {
        cpx wn(cos(pi / i), sin(f * pi / i)), t;
        for (int j = 0; j < m; j += (i << 1)) {
            cpx w(1, 0);
            for (int k = 0; k < i; ++k, w *= wn) {
                t = w * a[j + i + k];
                a[j + i + k] = a[j + k] - t;
                a[j + k] = a[j + k] + t;
            }
        }
    }
}

int main()
{
    n = gi();
    for (int x, i = 1; i <= n; ++i) {
        x = gi(); maxa = max(maxa, x);
        a[x] = 1; b[x + x] = 1; c[x + x + x] = 1;
    }
    n = maxa * 3;

    for (m = 1; m <= n; m <<= 1) ++L;
    for (int i = 0; i < m; ++i) R[i] = (R[i >> 1] >> 1) | ((i & 1) << (L - 1));
    fft(a, 1); fft(b, 1); fft(c, 1);
    for (int i = 0; i < m; ++i) {
        ans[i] += (a[i] * a[i] * a[i] - a[i] * b[i] * 3. + c[i] * 2.) / 6.;
        ans[i] += (a[i] * a[i] - b[i]) / 2.;
        ans[i] += a[i];
    }
    fft(ans, -1);

    for (int i = 0; i < m; ++i) {
        int x = (int)(ans[i].real() / m + 0.5);
        if (x) printf("%d %d\n", i, x);
    }
    return 0;
}