【dfs】hdu 1111 Secret Code
题目链接:hdu 1111 Secret Code
题目大意
输入一个复数X=xr+xi以及复数B=br+bi,求数列a[n]使得x=a0+a1*b+a2*b^2+…an*b^n。其中|Xi| <= 1000000, |Bi| <= 16,n<=100,
|b|>ai>=0,|b|>1.
很少有人问津的搜索题目,一开始我以为20^100绝对会超时,数据又可能溢出,结果int直接就过了,时间也还好。
说一下思路
- 暴力搜索,把要搜索的序列{an}想成一个森林,a0在第一层,a1在下一层,以此类推
- 秦九韶算法:x=a0+(a1+(a2+(a3+…)*b)*b)*b,容易想到递归实现
- 为了方便运算,把实部、虚部拆开一起递归,每次递归枚举下一层a(i+1)的所有可能结果,很容易想到{an}必定都为整数序列,每枚举一个a(i+1),判断[X-a(i)]%b能否整除(实部与实部、虚部与虚部),满足继续dfs,不满足继续枚举,直到Xr==0,Xi==0递归结束,当然不能超过100层。最后得到的序列逆序输出即可
参考代码
//*Author:Hacker_vision*/
#include<bits/stdc++.h>
using namespace std;
const int _max=1e3+10;
int Xr,Xi,Br,Bi,limit;
bool ok;
int tot,ans[_max];
void dfs(int Xr,int Xi,int deep){
if(deep>100||ok)return;
if(Xr==0&&Xi==0){
tot=deep;//记录序列{an}的长度
ok=true;//标记ok判断存在性
return;
}
int a,b;//中间变量而已
for(int i = 0; i*i < limit; ++ i) {
ans[deep]=i;
a=(Xr-i)*Br+Xi*Bi;
b=Xi*Br-(Xr-i)*Bi;
if(a%limit==0&&b%limit==0){ //序列{an}一定都为整数,整除是必要的
dfs(a/limit,b/limit,deep+1);
if(ok) return;
}
}
}
int main(){
// freopen("input.txt","r",stdin);
int T;cin>>T;
while(T--){
scanf("%d%d%d%d",&Xr,&Xi,&Br,&Bi);
limit=Br*Br+Bi*Bi;
ok=false;
dfs(Xr,Xi,0);
if(!ok) {puts("The code cannot be decrypted.");continue;}
printf("%d",ans[tot-1]);//逆序输出
for( int i=tot-2;i>=0;i--) printf(",%d",ans[i]);puts("");
}
return 0;
}
- 加粗
Ctrl + B
- 斜体
Ctrl + I
- 引用
Ctrl + Q
- 插入链接
Ctrl + L
- 插入代码
Ctrl + K
- 插入图片
Ctrl + G
- 提升标题
Ctrl + H
- 有序列表
Ctrl + O
- 无序列表
Ctrl + U
- 横线
Ctrl + R
- 撤销
Ctrl + Z
- 重做
Ctrl + Y