洛谷 P4106 / bzoj 3614 [ HEOI 2014 ] 逻辑翻译 —— 思路+递归

题目：https://www.luogu.org/problemnew/show/P4106

https://www.lydsy.com/JudgeOnline/problem.php?id=3614

从很小的情况考虑，看题面上的样例：

x1=+1 x2=+1  0

x1=+1 x2=-1  1

x1=-1 x2=+1  2

x1=-1 x2=-1  3

手算的做法应该是设原式为 a0 + a1x1 + a2x2 + a3x1x2

通过加减，把关于 x2,x1 的项逐渐消去，得到常数项后再回代，逐个得出答案；

然后我们发现加减的两个式子的位置也是有规律的，如果把-1看作0，+1看作1，式子的系数从小到大排序：

x1=-1 x2=-1 3    —— 1

x1=-1 x2=+1 2    —— 2

x1=+1 x2=-1 1    —— 3

x1=+1 x2=+1 0    —— 4

12相加得到 x1 = -1 时无 x2 的值(+a0)，34相加得到 x1 = 1 时无 x2 的值(+a0)；

21相减得到 x2 = 1 时无 x1 的值(-a3)，43相减得到 x2 = 1 时无 x1 的值(+a3)；

两个相加的结果再相加得出a0，相减得出a1，两个相减的结果同理；

可以发现，这样做其实就是把问题的规模变成了一半，可以递归求解；

但空间很小，为了不递归，可以仿照FFT，FWT等，通过位置安排来直接循环做；

实际上，最初的相邻两位置加减，让 xn 的有无与该位置末位的 0/1 对应；

以此类推，如果倍增加减，则每个 xi  的有无就和该位置的第 i 位的 0/1 对应；

“xi 的有无”意思是这个值中关于 xi 的项还是否存在；

所以最后答案的 xi 情况就和其位置上的 0/1 分布情况一样，像FWT一样跳着加减就能得到；

过程中不要除以2，最后输出时把 2^n 放在分母上；

递归输出即可按照字典序，详见代码囧

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef double db;
int const xn=(<<);
int n,lim,f[xn];
char dc[];
int gcd(int a,int b){return b?gcd(b,a%b):a;}
void get()
{
  scanf("%s",dc); int t=; db ff;
  for(int j=;j<n;j++){t<<=; if(dc[j]=='+')t|=;}
  scanf("%lf",&ff);
  if(ff>)f[t]=(int)(ff*+0.5);
  else f[t]=(int)(ff*-0.5);//!
}
void print(int x,int s)
{
  if(x==)return; if(x<)putchar('-'),x=-x;//
  int t=(<<n),g=gcd(x,t);//t:100*(/2)^n
  if(g==t)printf("%d",x/g); else printf("%d/%d",x/g,t/g);
  if(!s){puts(""); return;} putchar(' ');
  for(int i=n-,nw=;i>=;i--,nw++)if(s&(<<i))printf("x%d",nw);
  puts("");
}
void work(int nw,int s)
{
  int t=((s<<|)<<(n-nw));
  print(f[t],t); if(nw==n)return;
  work(nw+,(s<<)|);
  work(nw+,(s<<));
}
int main()
{
  scanf("%d",&n); lim=(<<n);
  for(int i=;i<=lim;i++)get();
  for(int i=;i<n;i++)
    {
      int t=(<<i),d=(t<<);
      for(int j=;j<lim;j+=d)
    for(int k=;k<t;k++)
      {
        int x=f[j+k],y=f[j+t+k];
        f[j+k]=x+y; f[j+t+k]=y-x;
      }
    }
  print(f[],); work(,);
  return ;
}