【luogu P6800】Chirp Z-Transform(多项式)(NTT)(bluestein 算法)

Chirp Z-Transform

题目链接:luogu P6800

题目大意

给你一个多项式和 c,m,要你求把 c^0,c^1,…c,^m-1 分别带入多项式得到的值。

思路

考虑把答案也看做是多项式:
a n s i = F ( c i ) = ∑ j = 0 n − 1 a j c i j ans_i=F(c^i)=\sum\limits_{j=0}^{n-1}a_jc^{ij} ansi=F(ci)=j=0n1ajcij

然后又一个东西是: i j = ( i + j 2 ) − ( i 2 ) − ( j 2 ) ij=\binom{i+j}{2}-\binom{i}{2}-\binom{j}{2} ij=(2i+j)(2i)(2j)
简单证明:
( i + j 2 ) − ( i 2 ) − ( j 2 ) = ( i + j ) ( i + j − 1 ) − i ( i − 1 ) − j ( j − 1 ) 2 \binom{i+j}{2}-\binom{i}{2}-\binom{j}{2}=\dfrac{(i+j)(i+j-1)-i(i-1)-j(j-1)}{2} (2i+j)(2i)(2j)=2(i+j)(i+

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1. 数据准备 首先,我们需要准备一些地理坐标数据。这里我们使用一个包含1000个地理坐标的数据集,保存在一个名为“coords.mat”的文件中。 2. 导入数据 我们使用Matlab的load函数导入数据,将数据保存在一个名为“coords”的变量中。代码如下: load('coords.mat'); coords = coords'; 这里需要注意的是,我们需要将数据转置一下,使得每个地理坐标保存在一行中。 3. 初始化聚类中心 接下来,我们需要随机地初始化k个聚类中心。这里我们选择k=5,即将数据集分成5类。初始化聚类中心的代码如下: k = 5; idx = randperm(size(coords, 1)); centroids = coords(idx(1:k), :); 4. 迭代更新聚类中心 接下来,我们需要迭代更新聚类中心,直到聚类中心不再变化为止。在每一轮迭代中,我们需要将每个地理坐标分配到最近的聚类中心所在的类别中,并更新聚类中心的坐标。代码如下: max_iters = 10; for i=1:max_iters % 计算每个地理坐标与各个聚类中心的距离 dists = pdist2(coords, centroids); % 将每个地理坐标分配到最近的聚类中心所在的类别中 [~, idx] = min(dists, [], 2); % 更新聚类中心的坐标 for j=1:k centroids(j, :) = mean(coords(idx==j, :), 1); end end 5. 可视化聚类结果 最后,我们可以将聚类结果可视化出来。代码如下: colors = {'r', 'g', 'b', 'c', 'm'}; figure; hold on; for i=1:k scatter(coords(idx==i, 1), coords(idx==i, 2), 10, colors{i}, 'filled'); end scatter(centroids(:, 1), centroids(:, 2), 50, 'k', 'filled'); hold off; 运行完整的代码后,我们可以得到如下的聚类结果: ![k-means示例图](https://cdn.luogu.com.cn/upload/image_hosting/cv15pj8u.png)

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值