题目大意:
输入一个复数X=xr+xi以及复数B=br+bi,求数列a[n]使得x=a0+a1b+a2b2+…an*bn
我们对x这个式子进行提公因式的操作,其实就是 秦九韶算法,x=a0+(a1+(a2+(a3+…)*b)*b)*b
对于这个我们一般就使用递归算法, 挨个寻找a[i] 进行试验,n<=100 这个也是一个比较重要的条件了。
dfs时,如果 x-a[i] 可以被 b整除,那么这个 a[i]暂时符合,继续dfs a[i+1] ,如果到最后,n<=100 那么说明前面得到的可取,倒序输出就好了。
如果n>100, 那么前面的不可取,继续退回去dfs。
对于复数的除法
设复数 a + bi , c + di ;
con = c * c + d * d ;
则 (a + bi) / (c + di ) = (ac + bd) / con+ (bc - ad) / con * i ;
//HDU 1111
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <deque>
#define ll long long
const int maxn=2e2+5;
const int inf = 0x3f3f3f3f;
using namespace std;
int xr,xi,br,bi,con;
int flag,t;
int a[maxn];
void dfs(int n){
int x,y,i;
if(n>100) return ;
if(xr==0&& xi==0) {flag=1;t=n;return;}//x都为0了说明到底了。
for( i=0;i*i<con;i++){
/*
设复数 a + bi , c + di ;
t = c * c + d * d ;
则 (a + bi) / (c + di ) = (ac + bd) / t + (bc - ad) / t * i ;
*/
x=(xr-i)*br+xi*bi; //i为所求整数列
y=xi*br-(xr-i)*bi;
a[n]=i;
if(x%con==0 && y%con==0)//因为a[i]都是整数,保证可以整除
{
xr=x/con;xi=y/con;
dfs(n+1);
}
if(flag) return;
}
}
int main(){
int T;
scanf("%d",&T);
while (T--)
{
scanf("%d%d%d%d",&xr,&xi,&br,&bi);
con=br*br+bi*bi;//复数的除法运算需要使用的公式。
flag=0;
dfs(0);
if(!flag) printf("The code cannot be decrypted.\n");
else{
/*for(int i=t-1;i>=0;i--){
if(i!=0) printf("%d,",a[i]);
else printf("%d\n",a[i]);
}*/
cout<<a[t-1];
for(int i=t-2;i>=0;i--)
cout<<','<<a[i];
cout<<endl;
}
}
return 0;
}