[hdu 1111] Secret Code(dfs + 秦九韶算法)

在这里插入图片描述

题目大意:
输入一个复数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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值