第九集,结束亦是开始

题意:

大致意思就是给你n个3进制的数字,让你计算有多少对数字的哈夫曼距离等于i(0<=i<=2^m)

思路:

这个是一个防ak题,做法是要手推公式的fwt

大概就这个意思

把n个数字标记到大小为3^m的数组里

然后一个简单的方法就是,假设a是标记数组

for i=0 i<3^m i++

for j=0 j<3^m j++

ans[dis(a[i],a[j])]+=a[i]*a[j]

可能i==j的时候被算重复了,大概特判减去一下n就行了

我们发现,如果dis(a[i],a[j])是位运算就好了,这样就直接拍个fwt就行了

然后我们注意到dis(a[i],a[j])其实可以看成是一个三进制的位运算,这时候我们要推一下公式

我们定义一下位运算法则假设符号为*,则i*j=abs(i-j),且i,j∈[0,2]

给出我当时推导的过程:

a=0+2

b=0-2

c=1

d=0+1+2

C0=(a*a+b*b)/2+c*c

C1=d*d-C0-C2

C2=(a*a-b*b)/2

大概就是把一个长度为n的多项式乘法变成了4个长度为n/3的多项式乘法

0,1,2是卷积前,最高位三进制为0,1,2的三个长度为n/3的多项式

然后我们构造a,b,c,d四个长度为n/3的多项式

C0,C1,C2是卷积后的多项式

这时候我们只需要计算a*a,b*b,c*c,d*d就行了!

至于怎么构造出来a,b,c,d的,可能就看感觉了吧。。。。。。

然后分析一下复杂度

然后就看出来这是一个公比为4/3的等比数列,求和一下就知道T(m)=4^m-3^m

时间复杂度为O(m)=4^m

代码实现

#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = ;
ll pool[/+],*p=pool;
ll ans[];
ll A[maxn];
/*
a=0+2
b=0-2
c=1
d=0+1+2
C0=(a*a+b*b)/2+c*c
C1=d*d-C0-C2
C2=(a*a-b*b)/2
*/
void fwt(ll *A,int len){
if(len==){
A[]=A[]*A[];
return ;
}
int nlen=len/;
ll *a=p;
p+=nlen;
ll *b=A;
ll *c=A+nlen;
ll *d=A+nlen*;
for(int i=;i<nlen;i++){
ll x=A[i],y=A[i+nlen],z=A[i+*nlen];
a[i]=x+z;
b[i]=x-z;
c[i]=y;
d[i]=x+y+z;
}
fwt(a,nlen);fwt(b,nlen);
fwt(c,nlen);fwt(d,nlen);
for(int i=;i<nlen;i++){
ll x=a[i],y=b[i],z=c[i],w=d[i];
A[i]=(x+y)/+z;
A[i+nlen*]=(x-y)/;
A[i+nlen]=w-A[i+nlen*]-A[i];
}
} int cal(int x){
int res=;
while(x){
res+=x%;
x/=;
}
return res;
} int main(){
int n,m;
scanf("%d%d",&n,&m);
for(int i=,x;i<=n;i++){
scanf("%d",&x);
A[x]++;
}
int len=;
for(int i=;i<=m;i++)
len*=;
fwt(A,len);
for(int i=;i<len;i++)
ans[cal(i)]+=A[i];
ans[]-=n;
for(int i=;i<=m*;i++){
if(i)putchar(' ');
printf("%lld",ans[i]);
}puts("");
return ;
}

福建工程学院第十四届ACM校赛M题题解 fwt进阶,手推三进制fwt的更多相关文章

  1. 福建工程学院第十四届ACM校赛J题题解

    第六集,想不到你这个浓眉大眼的都叛变革命了 题意: 给你两个只包含01的字符串S和T,问你在允许一次错误的情况下,T是否能成为S的子串 思路: 这个问题的解法挺多,我是用fft匹配的,也比较简单,针对 ...

  2. 福建工程学院第十四届ACM校赛G题题解

    外传:编剧说了不玩游戏不行 题意: 有n个石堆,我每次只能从某一堆中取偶数个石子,你取奇数个,我先手,先不能操作的人输.问最后谁能赢. 思路: 这个题仔细想想,就发现,取奇数的人有巨大的优势,因为假设 ...

  3. 福建工程学院第十四届ACM校赛B题题解

    第二集,未来的我发量这么捉急的吗 题意: 有n个数,请问有多少对数字(i,j)(1<=i<j<=n),满足(a[i]^a[j])+((a[i]&a[j])<<1) ...

  4. 福建工程学院第十四届ACM程序设计大赛 - E - 外传:小晋逃生记

    http://www.fjutacm.com/Contest.jsp?cid=705#P4 其实想清楚了就很简单,之前想了很多种方法,以为是二分什么的,看起来就像是一个单峰函数.但是发现直接暴力一波就 ...

  5. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  6. 湖南大学第十四届ACM程序设计新生杯(重现赛)G a+b+c+d=? (16进制与LL范围)

    链接:https://ac.nowcoder.com/acm/contest/338/G来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 32768K,其他语言65536K6 ...

  7. 湖南大学第十四届ACM程序设计新生杯(重现赛)I:II play with GG(博弈论||DP)

    链接:https://ac.nowcoder.com/acm/contest/338/I 来源:牛客网 题目描述 IG won the S championship and many people a ...

  8. 湖南大学第十四届ACM程序设计新生杯(重现赛)

    RANK  0 题数 0 期末复习没有参加,补几道喜欢的题. A: AFei Loves Magic  签到 思路 :不需考虑 碰撞 直接计算最终状态即可. #include<bits/stdc ...

  9. 湖南大学第十四届ACM程序设计新生杯 E.Easy Problem

    E.Easy Problem Description: Zghh likes number, but he doesn't like writing problem description. So h ...

随机推荐

  1. UEFI和GPT

    好就没用linux了,这几天在win8笔记本上用虚拟机装了下,也准备装到硬盘上和win8双系统使用,发现一些概念已经跟不上时代了. 一个是在虚拟机中装的时候,分配了虚拟硬盘分区时,提示选择分区表类型, ...

  2. 按下home键,重新打开,应用重启

    其实不是重启,只是重新打开了luncher的那个activity.只要通过判断把它finish,就会显示按下home键前的页面. 解决方法: 在重启的页面中加入一下代码,注意加在setContentV ...

  3. LC 963. Minimum Area Rectangle II

    Given a set of points in the xy-plane, determine the minimum area of any rectangle formed from these ...

  4. [Err] ORA-00942: table or view does not exist

    [Err] ORA-00942: table or view does not exist 当前用户加表明 例如:SCOTT."replyInfo"

  5. openstack核心组件--keystone(1)

    一.keystone介绍 keystone 是OpenStack的组件之一,用于为OpenStack家族中的其它组件成员提供统一的认证服务,包括身份验证.令牌的发放和校验.服务列表.用户权限的定义等等 ...

  6. zabbix(2)使用指南

    一.邮件报警(一个客户端安装server,agent) 管理->报警媒介类型->email 管理->用户->Admin->报警媒介 配置->动作->Repor ...

  7. RxJava2实战---第五章 变换操作符和过滤操作符

    RxJava2实战---第五章 变换操作符和过滤操作符 RxJava的变换操作符主要包括以下几种: map():对序列的每一项都用一个函数来变换Observable发射的数据序列. flatMap() ...

  8. windows客户端如果通过cmd窗口连接到远程linux服务器,可以使用telnet;

     linux系统打开telnet端口的方法 2016-03-11 16:02:25 标签:linux telnet 原创作品,允许转载,转载时请务必以超链接形式标明文章 原始出处 .作者信息和本声明. ...

  9. MSSQL数据库日志满的快速解决办法

    先提供一种复杂的方法压缩日志及数据库文件如下:   1.清空日志 DUMP TRANSACTION 库名 WITH NO_LOG  2.截断事务日志: BACKUP LOG 数据库名 WITH NO_ ...

  10. 【.NET】ASP.Net IE10+ SCRIPT:XXX_doPostBack 未定义

    问题描述 GridView中分页控件,点击分页无反应,Linkbutton点击无反应,打开Web控制台,发现如下错误:SCRIPTXXX:_doPostBack 未定义:查询后得知,是由于.NET F ...