洛谷P1919 【模板】A*B Problem升级版 题解(FFT的第一次实战)

洛谷P1919 【模板】A*B Problem升级版(FFT快速傅里叶)

刚学了FFT,我们来刷一道模板题。

题目描述

给定两个长度为 n 的两个十进制数,求它们的乘积。

n<=100000

如果用 n^2 暴力,肯定会 TLE。

我们把这两个数看成一个多项式。

f(x)=a0+a1*101+a2*102+a3*103+ ...... +an-1*10n-1

然后就可以愉快的FFT求解了!!

 1 #include<iostream>
 2 #include<cmath>
 3 #include<complex>
 4 #include<cstdio>
 5 using namespace std;
 6 
 7 const double PI=acos(-1);
 8 typedef complex<double> cmplx;
 9 int rev[1000005];
10 cmplx a[1000005],b[1000005];
11 int n,bit=2,output[1000005];
12 char s1[60005],s2[60005];
13 
14 void get_rev(){
15     for(int i=0;i<bit;i++)
16     rev[i]=(rev[i>>1]>>1)|(bit>>1)*(i&1);
17 }
18 
19 void FFT(cmplx *a,int dft){
20     for(int i=0;i<bit;i++)
21         if(i<rev[i])swap(a[i],a[rev[i]]);
22     for(int i=1;i<bit;i<<=1){
23         cmplx W=exp(cmplx(0,dft*PI/i));
24         for(int j=0;j<bit;j+=i<<1){
25             cmplx w(1,0);
26             for(int k=j;k<j+i;k++,w*=W){
27                 cmplx x=a[k];
28                 cmplx y=w*a[k+i];
29                 a[k]=x+y;
30                 a[k+i]=x-y;
31             }
32         }
33     }
34     if(!~dft)for(int i=0;i<bit;i++)a[i]/=bit;
35 }
36 
37 int main(){
38     scanf("%d%s%s",&n,s1,s2);
39     for(int i=1;(1<<i)<=2*n;i++)bit<<=1;
40     for(int i=0;i<n;i++){
41         a[i]=s1[n-i-1]-'0';
42         b[i]=s2[n-i-1]-'0';
43     }
44     get_rev();
45     FFT(a,1),FFT(b,1);
46     for(int i=0;i<bit;i++)a[i]=a[i]*b[i];
47     FFT(a,-1);
48     for(int i=0;i<bit;i++){    //确保输出为十进制 
49         output[i]+=a[i].real()+0.5;
50         output[i+1]+=output[i]/10;
51         output[i]%=10;
52     }
53     bool check=false;
54     for(int i=n<<1;i>=0;i--){    //去前导零输出 
55         if(check||output[i]){
56             printf("%d",output[i]);
57             check=true;
58         }
59     }
60     if(!check)printf("0");
61 }

 FFT学习笔记

转载于:https://www.cnblogs.com/ezoiLZH/p/8997960.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述似乎缺失了关键信息,通常我会需要了解“P10780 食物”是什么具体的算法竞赛题目,它来自在线平台洛谷(Luogu),以及该题目的大致背景、条件和目标。洛谷食物(Food)可能是某种数据结构或算法问题,比如贪吃蛇、分配任务等。 然而,我可以给你提供一个通用的模板: **[洛谷 P10780 食物 - 题目解析]** 题目名称:P10780 食物(假设是关于食物分配或者饥饿游戏的问题) 链接:[插入实际题目链接] **背景:** 此题通常涉及动态规划或者搜索策略。场景可能是有n个参与者(选手或角色),每个都有特定的食物需求或者优先级,我们需要在有限的食物资源下合理分配。 **分析:** 1. **输入理解**:首先读入n个参与者的信息,包括每个人的需求量或优先级。 2. **状态定义**:可以定义dp[i][j]表示前i个人分配完成后剩余的食物能满足第j个人的最大程度。 3. **状态转移**:递推式可能涉及到选择当前人分配最多食物的版本,然后更新剩余的食物数。 4. **边界条件**:如果剩余食物不足以满足某人的需求,则考虑无法分配给他;如果没有食物,状态值设为0。 5. **优化策略**:可能需要对状态数组进行滚动更新,以减少空间复杂度。 **代码示例(伪代码或部分关键代码片段):** ```python # 假设函数分配_food(demand, remaining)计算分配给一个人后剩余的食物 def solve(foods): dp = [[0 for _ in range(max_demand + 1)] for _ in range(n)] dp = foods[:] # 从第一个到最后一个参与者处理 for i in range(1, n): for j in range(1, max_demand + 1): if dp[i-1][j] > 0: dp[i][j] = max(dp[i][j], dp[i-1][j] - foods[i]) dp[i][j] = max(dp[i][j], distribute_food_to(i, dp[i-1][j])) return dp[n-1][max_demand] ``` **相关问题--:** 1. 这道题是如何运用动态规划的? 2. 如果有优先级限制,应该如何调整代码? 3. 怎样设计搜索策略来解决类似问题?
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值