BZOJ 3771: Triple(生成函數 FFT)

来源:https://www.cnblogs.com/zwfymqz/archive/2018/06/07/9152951.html
-Advertisement-
Play Games

Description 我們講一個悲傷的故事。 從前有一個貧窮的樵夫在河邊砍柴。 這時候河裡出現了一個水神,奪過了他的斧頭,說: “這把斧頭,是不是你的?” 樵夫一看:“是啊是啊!” 水神把斧頭扔在一邊,又拿起一個東西問: “這把斧頭,是不是你的?” 樵夫看不清楚,但又怕真的是自己的斧頭,只好又答: ...


Time Limit: 20 Sec  Memory Limit: 64 MB
Submit: 911  Solved: 528
[Submit][Status][Discuss]

Description

我們講一個悲傷的故事。 從前有一個貧窮的樵夫在河邊砍柴。 這時候河裡出現了一個水神,奪過了他的斧頭,說: “這把斧頭,是不是你的?” 樵夫一看:“是啊是啊!” 水神把斧頭扔在一邊,又拿起一個東西問: “這把斧頭,是不是你的?” 樵夫看不清楚,但又怕真的是自己的斧頭,只好又答:“是啊是啊!” 水神又把手上的東西扔在一邊,拿起第三個東西問: “這把斧頭,是不是你的?” 樵夫還是看不清楚,但是他覺得再這樣下去他就沒法砍柴了。 於是他又一次答:“是啊是啊!真的是!” 水神看著他,哈哈大笑道: “你看看你現在的樣子,真是醜陋!” 之後就消失了。   樵夫覺得很坑爹,他今天不僅沒有砍到柴,還丟了一把斧頭給那個水神。 於是他準備回家換一把斧頭。 回家之後他才發現真正坑爹的事情才剛開始。 水神拿著的的確是他的斧頭。 但是不一定是他拿出去的那把,還有可能是水神不知道怎麼偷偷從他家裡拿走的。 換句話說,水神可能拿走了他的一把,兩把或者三把斧頭。   樵夫覺得今天真是倒霉透了,但不管怎麼樣日子還得過。 他想統計他的損失。 樵夫的每一把斧頭都有一個價值,不同斧頭的價值不同。總損失就是丟掉的斧頭價值和。 他想對於每個可能的總損失,計算有幾種可能的方案。 註意:如果水神拿走了兩把斧頭a和b,(a,b)和(b,a)視為一種方案。拿走三把斧頭時,(a,b,c),(b,c,a),(c,a,b),(c,b,a),(b,a,c),(a,c,b)視為一種方案。  

 

Input

第一行是整數N,表示有N把斧頭。 接下來n行升序輸入N個數字Ai,表示每把斧頭的價值。  

 

Output

若幹行,按升序對於所有可能的總損失輸出一行x y,x為損失值,y為方案數。  

 

Sample Input

4
4
5
6
7

Sample Output

4 1
5 1
6 1
7 1
9 1
10 1
11 2
12 1
13 1
15 1
16 1
17 1
18 1
樣例解釋
11有兩種方案是4+7和5+6,其他損失值都有唯一方案,例如4=4,5=5,10=4+6,18=5+6+7.

HINT

 

所有數據滿足:Ai<=40000


 

 

Source

應該不難看出是生成函數

我們用$A(x) = a + bx^1 + cx^2 + \dots $表示價值為$1$的方案為$a$,價值為$2$的方案為$b$

那麼很顯然的思路就是:$A(x) + \frac{A(x) * A(x)}{2!} + \frac{A(x) * A(x) * A(x)}{3!}$

但是題目中要求了每種斧子只能拿一次,這樣會多計算重覆拿的一部分

因此我們考慮用容斥的方法將他們減去

定義$B(x) = x ^ i$表示價值為$i$的拿了兩把的方案數

$C(x) = x ^ i$表示價值為$i$的拿了三把的方案數

拿兩把斧子時會計算到$(x, x)$這種情況,所以拿兩把時應該為$\frac{A(x) * A(x) - B(x)}{2!}$

拿三把時有些複雜

我們需要減去$(x, x, x)$和$(x, y, y)$這兩種情況

第一種情況就是$C(x)$

第二種情況可以通過$A(x)* B(x)$計算得到,但此時也計算上了$(x, x, x)$這種情況

$(x, y, y)$有三種組合排列方式,所以要乘$3$,但$(x, x, x)$只有一種排列方式,所以最終統計答案時要加上$2 * C(x)$

最終的答案就是

$A + \frac{A * A - B}{2!} + \frac{A * A * A - 3 * A * B + 2C}{3!}$

多項式乘法可以用NTT,不過模數會炸998244353

看到大佬們都用FFT A了,那我就偷個懶嘍

 

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
//#include<iostream>
const double pi = acos(-1);
using namespace std;
const int MAXN = 150000;
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 * 10 + c - '0', c = getchar();
    return x * f;
}
struct complex {
    double x, y;
    complex(double xx = 0, double yy = 0) {x = xx, y = yy;}
    complex operator + (const complex &rhs) {
        return complex(x + rhs.x, y + rhs.y);
    }
    complex operator - (const complex &rhs) {
        return complex(x - rhs.x, y - rhs.y);
    }
    complex operator * (const complex &rhs) {
        return complex(x * rhs.x - y * rhs.y, x * rhs.y + y * rhs.x);
    }
    complex operator * (const double &rhs) {
        return complex(x * rhs, y * rhs);
    }
    complex operator / (const double &rhs) {
        return complex(x / rhs, y / rhs);
    }
}A[MAXN], B[MAXN], C[MAXN];
int val, n, N, L, len, r[MAXN];
void FFT(complex *A, int type) {
    for(int i = 0; i < N; i++) if(i < r[i]) swap(A[i], A[r[i]]);
    for(int mid = 1; mid < N; mid <<= 1) {
        complex Wn(cos(pi / mid), type * sin(pi / mid));
        for(int j = 0; j < N; j += (mid << 1)) {
            complex w = complex(1, 0);
            for(int i = 0; i < mid; i++, w = w * Wn) {
                complex x = A[j + i], y = w * A[j + i + mid];
                A[j + i] = x + y;
                A[j + i + mid] = x - y;
            }
        }
    }
    if(type == -1) {
        for(int i = 0; i < N; i++)
            A[i].x /= N;
    }
}
void print(complex *a) {
    for(int i = 0; i < N; i++)
        printf("%d %lf %lf\n", i, a[i].x, a[i].y);
}
int main() {
#ifdef WIN32
    freopen("a.in", "r", stdin);
    freopen("b.out", "w", stdout);
#endif
    n = read();
    for(int i = 1; i <= n; i++)  
        val = read(), 
        A[val].x = 1, 
        B[2 * val].x = 1, 
        C[3 * val].x = 1,
        len = max(3 * val, len);
    len = len + 1;//tag
    for(N = 1; N <= len; N <<= 1, L++);
    for(int i = 0; i < N; i++)
        r[i] = (r[i >> 1] >> 1) | ((i & 1) << (L - 1));
        
    FFT(C, 1);
    FFT(A, 1); 
    FFT(B, 1); 
    for(int i = 0; i < N; i++)
        A[i] = A[i] + (A[i] * A[i] - B[i]) / 2.0 + (A[i] * A[i] * A[i] - A[i] * B[i] * 3.0 + C[i] * 2.0) / 6.0;
    FFT(A, -1);
    
    for(int i = 0; i < N; i++) {
        long long cur = (long long )(A[i].x + 0.5);
        if(cur) 
            printf("%d %lld\n", i, cur);
    }
    return 0;
}

 


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 11.1.前臺用戶模型創建 (1)apps/front/models.py 首先安裝:pip install shortuuid (2)manage.py (3)生成表和添加前臺用戶 添加用戶 11.2.註冊界面完成 (1)front/views.py (2)common/images/logo.j ...
  • 因為自己有一些編程經驗,所以在學習Python的過程中過於急功近利,導致自己現在基礎不牢,所以專門又回過頭來複習一遍Python的基礎知識點。算是對自己的學習Python的整理和總結。夯實基礎,打牢基礎,才能走的更遠。 第一部分我打算著重總結一下Python的內置的數據結構。下麵的代碼都是親測可用的 ...
  • 自學記錄: 1、字元串 2、換行/不換行 3、判斷類型 面向對象 ...
  • Java面向對象的三大特性之封裝一:面向對象的三大特性;二:封裝的概念;三:封裝的好處;四:實現Java封裝的步驟; ...
  • 題目描述 給定一個字元串,要求把字元串前面的若幹個字元移動到字元串的尾部,如把字元串“abcdef”前面的2個字元'a'和'b'移動到字元串的尾部,使得原字元串變成字元串“cdefab”。請寫一個函數完成此功能,要求對長度為n的字元串操作的時間複雜度為 O(n),空間複雜度為 O(1)。 分析與解法 ...
  • 1.JDK中的Executor框架是基於生產者-消費者模式的線程池,提交任務的線程是生產者,執行任務的線程是消費者。 Executor線程池可以用於非同步任務執行,而且支持很多不同類型任務執行策略,同時為任務提交和任務執行之間的解耦提供了標準方法。 Executor線程池支持如下三種線程執行策略: ( ...
  • 在電腦科學中,二叉樹是每個結點最多有兩個子樹的樹結構。通常子樹被稱作“左子樹”(left subtree)和“右子樹”(right subtree)。二叉樹常被用於實現二叉查找樹和二叉堆。 ...
  • 需要先把Oracle包倒到項目中ojdbc14-10.2.0.1.0.jar ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發(一):使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說,近來想嘗試開發一款安卓APP,考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP,畢竟是使用Visual Studio開發工具,使用起來也比較的順手,結合微軟官方的教程進行了安卓 ...
  • 前言 QuestPDF 是一個開源 .NET 庫,用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫,它徹底改變了我們生成 PDF 文檔的方 ...
  • 項目地址 項目後端地址: https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址: ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址: http://tree ...
  • 話不多說,直接開乾 一.下載 1.官方鏈接下載: https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe,運行開始下載SQL server; 二. ...
  • 前言 隨著物聯網(IoT)技術的迅猛發展,MQTT(消息隊列遙測傳輸)協議憑藉其輕量級和高效性,已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫,為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog支持多種接收器用於日誌存儲,增強器用於添加屬性,LogContext管理動態屬性,支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項,適用於不同需求。 ...
  • 目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據,爬取信息時需要模擬瀏覽器行為,否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下: 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • 1.前言 什麼是熱更新 游戲或者軟體更新時,無需重新下載客戶端進行安裝,而是在應用程式啟動的情況下,在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR,Xlua,ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles,Addressable, ...
  • 本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺,其實像阿裡雲,騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後,它會送 ...
  • 通過以下方式可以高效,並保證數據同步的可靠性 1.API設計 使用RESTful設計,確保API端點明確,並使用適當的HTTP方法(如POST用於創建,PUT用於更新)。 設計清晰的請求和響應模型,以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證,確保接收到的數據符合預期格 ...