POJ3150【FFT】

转移很好用矩阵表示.然而矩阵乘法复杂度是O(n^3)的.

很容易发现转移矩阵是【循环矩阵】.而且有一个美妙的性质:【循环矩阵 * 循环矩阵 = 循环矩阵】.

所以我们计算矩阵乘法的时候可以只计算第一行.剩下的可以由第一行递推得出.

一次乘法的复杂度降到了O(n^2).这是可以接受的.

#include <cstdio>  
#include <cstdlib>  
#include <cstring>  
#include <cmath>  
#include <ctime>  
#include <algorithm>  
#include <iostream>  
#include <fstream>  
#include <vector>  
#include <queue>  
#include <deque>  
#include <map>  
#include <set>  
#include <string>  
#define make make_pair  
#define fi first  
#define se second  
  
using namespace std;  
  
typedef long long ll;  
typedef unsigned long long ull;  
typedef pair<int,int> pii;  
  
const int maxn = 1000010;  
const int maxm = 1010;  
const int maxs = 30;  
const int inf = 0x3f3f3f3f;  
const int P = 1000000007;  
const double error = 1e-9;  
const double Pi = 3.141592653589793238462643383279;
  
inline int read()  
{  
	int x = 0, f = 1; char ch = getchar();  
	while (ch <= 47 || ch >= 58)
		f = (ch == 45 ? -1 : 1), ch = getchar();  
	while (ch >= 48 && ch <= 57)  
		x = x * 10 + ch - 48, ch = getchar();  
	return x * f;
}

struct matrix
{
	ll num[501][501];
} f, e;

int n, m, d, 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值