Advertisement
radoslav11

FFT with xor

May 24th, 2016
165
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 0.95 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. #define endl '\n'
  3.  
  4. #define int long long
  5.  
  6. using namespace std;
  7. const int MAXN = (1 << 20);
  8.  
  9. typedef vector<int> polynomial;
  10.  
  11. void fft(polynomial &a, int low, int high)
  12. {
  13.     if(low == high - 1)
  14.         return;
  15.  
  16.     int len = (high - low) / 2, mid = low + len;
  17.     fft(a, low, mid);
  18.     fft(a, mid, high);
  19.  
  20.     for(int i = low; i < mid; i++)
  21.     {
  22.         int x1 = a[i];
  23.         int x2 = a[i + len];
  24.  
  25.         a[i] = (x1 - x2);
  26.         a[i + len] = (x1 + x2);
  27.     }
  28. }
  29.  
  30. void inv_fft(polynomial &a, int low, int high)
  31. {
  32.     if(low == high - 1)
  33.         return;
  34.  
  35.     int len = (high - low) / 2, mid = low + len;
  36.  
  37.     for(int i = low; i < mid; i++)
  38.     {
  39.         int y1 = a[i];
  40.         int y2 = a[i + len];
  41.  
  42.         a[i] = (y1 + y2) / 2;
  43.         a[i + len] = (y2 - y1) / 2;
  44.     }
  45.    
  46.     inv_fft(a, low, mid);
  47.     inv_fft(a, mid, high);
  48. }
  49.  
  50. void read()
  51. {
  52.    
  53. }
  54.  
  55. void solve()
  56. {
  57.  
  58. }
  59.  
  60. #undef int
  61. int main()
  62. {
  63.     ios_base::sync_with_stdio(false);
  64.     cin.tie(NULL);
  65.  
  66.     read();
  67.     solve();
  68.     return 0;
  69. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement