码迷,mamicode.com
首页 > 其他好文 > 详细

多项式 - 快速沃尔什变换

时间:2019-02-09 22:48:25      阅读:210      评论:0      收藏:0      [点我收藏+]

标签:argc   open   多项式   ast   整数   freopen   line   sum   public   

\(·\)是一种适用于整数域的二元运算,则两多项式关于此运算的方式定义为 \(C_k = \sum_{i·j=k} A_i * B_j\),即 \(C=A·B\)

\(FWT\) 主要解决多项式的常见的三种二元位运算,在三种运算下分别构造出不同的变换方式,个人认为比 \(NTT\) 简单 好背 一些。形式与 \(NTT\) 近似。
 
没有新东西可说,直接放上洛谷模板题的代码好了:

#include <cmath>
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

const int mod = 998244353;
const int maxn = (1 << 17) + 10;
int n, m, limit, f_1[maxn], f_2[maxn], f_3[maxn], f_4[maxn], r[maxn];

class Walsh_transform_and {
public:
  inline void Fast_walsh_transform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1)
      for(int l = 0, length = mid << 1; l < limit; l = l + length)
        for(int k = 0; k < mid; ++k) a[l + k] = (a[l + k] + a[l + mid + k]) % mod;
  }

  inline void Fast_walsh_inverform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1)
      for(int l = 0, length = mid << 1; l < limit; l = l + length)
        for(int k = 0; k < mid; ++k) a[l + k] = (a[l + k] - a[l + mid + k] + mod) % mod;
  }
} calculate_and;

class Walsh_transform_or {
public:
  inline void Fast_walsh_transform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1)
      for(int l = 0, length = mid << 1; l < limit; l = l + length)
        for(int k = 0; k < mid; ++k) a[l + mid + k] = (a[l + mid + k] + a[l + k]) % mod;
  }

  inline void Fast_walsh_inverform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1)
      for(int l = 0, length = mid << 1; l < limit; l = l + length)
        for(int k = 0; k < mid; ++k) a[l + mid + k] = (a[l + mid + k] - a[l + k] + mod) % mod;
  }
} calculate_or;

class Walsh_transform_xor {
public:
  inline void Fast_walsh_transform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1) {
      for(int l = 0, length = mid << 1; l < limit; l = l + length) {
        for(int k = 0; k < mid; ++k) {
          int x = a[l + k], y = a[l + mid + k];
          a[l + k] = (x + y) % mod, a[l + mid + k] = (x - y + mod) % mod;
        }
      }
    }
  }

  inline void Fast_walsh_inverform(int *a) {
    for(int mid = 1; mid < limit; mid = mid << 1) {
      for(int l = 0, length = mid << 1; l < limit; l = l + length) {
        for(int k = 0; k < mid; ++k) {
          int x = a[l + k], y = a[l + mid + k];
          a[l + k] = (x + y) % mod, a[l + mid + k] = (x - y + mod) % mod;
          a[l + k] = (a[l + k] & 1) ? (a[l + k] + mod) >> 1 : a[l + k] >> 1;
          a[l + mid + k] = (a[l + mid + k] & 1) ? (a[l + mid + k] + mod) >> 1 : a[l + mid + k] >> 1;
        }
      }
    }
  }
} calculate_xor;

inline int read() {
  register char ch = 0; register int w = 0, x = 0;
  while( !isdigit(ch) ) w |= (ch == ‘-‘), ch = getchar();
  while( isdigit(ch) ) x = (x * 10) + (ch ^ 48), ch = getchar();
  return w ? -x : x;
}

int main(int argc, char const *argv[])
{
  freopen("..\\nanjolno.in", "r", stdin);
  freopen("..\\nanjolno.out", "w", stdout);

  n = read(), limit = (1 << n);
  for(int i = 0; i < limit; ++i) f_1[i] = read();
  for(int i = 0; i < limit; ++i) f_2[i] = read();

  for(int i = 0; i < limit; ++i) f_3[i] = f_1[i];
  for(int i = 0; i < limit; ++i) f_4[i] = f_2[i];
  calculate_or.Fast_walsh_transform(f_3);
  calculate_or.Fast_walsh_transform(f_4);
  for(int i = 0; i < limit; ++i) r[i] = (1ll * f_3[i] * f_4[i]) % mod;
  calculate_or.Fast_walsh_inverform(r);
  for(int i = 0; i < limit; ++i) printf("%d ", r[i]);
  printf("\n");

  for(int i = 0; i < limit; ++i) f_3[i] = f_1[i];
  for(int i = 0; i < limit; ++i) f_4[i] = f_2[i];
  calculate_and.Fast_walsh_transform(f_3);
  calculate_and.Fast_walsh_transform(f_4);
  for(int i = 0; i < limit; ++i) r[i] = (1ll * f_3[i] * f_4[i]) % mod;
  calculate_and.Fast_walsh_inverform(r);
  for(int i = 0; i < limit; ++i) printf("%d ", r[i]);
  printf("\n");

  for(int i = 0; i < limit; ++i) f_3[i] = f_1[i];
  for(int i = 0; i < limit; ++i) f_4[i] = f_2[i];
  calculate_xor.Fast_walsh_transform(f_3);
  calculate_xor.Fast_walsh_transform(f_4);
  for(int i = 0; i < limit; ++i) r[i] = (1ll * f_3[i] * f_4[i]) % mod;
  calculate_xor.Fast_walsh_inverform(r);
  for(int i = 0; i < limit; ++i) printf("%d ", r[i]);
  printf("\n");

  fclose(stdin), fclose(stdout);
  return 0;
}

 
                  天也醉樱花,云脚乱蹒跚。

多项式 - 快速沃尔什变换

标签:argc   open   多项式   ast   整数   freopen   line   sum   public   

原文地址:https://www.cnblogs.com/nanjoqin/p/10358107.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!