matu 534 约瑟夫游戏

题目:

用循环链表实现:N个乘客同乘一艘船,因为严重超载,加上风高浪大,危险万分,因此船长告诉乘客,只有将部分乘客投入海中,其余人才能幸免于难。

无奈,大家只得同意这种办法。于是N个人围成一圈(从1,2,3...N分别编号)。由编号为1的人开始,依次报数,数到第M人,便把他投入大海中,

然后再从他的下一个人数起,数到第M人,再将他扔到大海中,如此循环地进行,直到剩下K个乘客为止。按顺序依次输出被扔下大海的乘客的编号。

提交格式:

实现int * solve(int N,int M,int K)函数。

函数参数为乘客人数N、间隔人数M和剩余乘客人数K,1<=N<=1000,1<=M<=500000,0<=K<N。

函数返回值为按顺序被扔下大海乘客编号的数组。

请不要printf输出任何内容。

输入样例1:

9 3 2

输出样例1:

3 6 9 4 8 5 2

输入样例2:

12 5 6

输出样例2:

5 10 3 9 4 12


#include <iostream>
#include <stdlib.h>
using namespace std;
typedef struct Node {
    int data;
    Node* before;
    Node* next;
}Node;
int * solve(int N,int M,int K){
    int *res=(int *)malloc((N-K)*sizeof(int));
    Node* head=new Node();
    head->next=head;
    head->data=1;
    Node *tail = head;
    for(int i=0;i<N-1;i++){
        Node *new_node=new Node();
        new_node->data=i+2;
        new_node->next=head;
        new_node->before=tail;
        tail->next=new_node;
        tail=new_node;
    }
    head->before=tail;
    Node *current=head;
    int count;int r=0;
    int n=N,p=1;
    while(N!=K){
        count=1;
        if(count==M) {
        res[r++]=current->data;
        current->next->before=current->before;
        current->before->next=current->next;
        current=current->next;
        N--;
        continue;
        }
        while(count<M){
            current=current->next;
            count++;
        }
        res[r++]=current->data;
        current->next->before=current->before;
        current->before->next=current->next;
        N--;
        current=current->next;
    }
    return res;
}

使用双向循环链表解决,注意M==1的情况要单拎出来

怎么将代码转换成自己熟悉的语言Eigen::Isometry3d Get3DR_TransMatrix(const std::vector<Eigen::Vector3d>& srcPoints, const std::vector<Eigen::Vector3d>& dstPoints) { Eigen::Isometry3d T = Eigen::Isometry3d::Identity(); double srcSumX = 0.0f; double srcSumY = 0.0f; double srcSumZ = 0.0f; double dstSumX = 0.0f; double dstSumY = 0.0f; double dstSumZ = 0.0f; //至少三组点 if (srcPoints.size() != dstPoints.size() || srcPoints.size() < 3) { return T; } int pointsNum = srcPoints.size(); for (int i = 0; i < pointsNum; ++i) { srcSumX += srcPoints[i].x; srcSumY += srcPoints[i].y; srcSumZ += srcPoints[i].z; dstSumX += dstPoints[i].x; dstSumY += dstPoints[i].y; dstSumZ += dstPoints[i].z; } Eigen::Vector3d centerSrc, centerDst; centerSrc.x = double(srcSumX / pointsNum); centerSrc.y = double(srcSumY / pointsNum); centerSrc.z = double(srcSumZ / pointsNum); centerDst.x = double(dstSumX / pointsNum); centerDst.y = double(dstSumY / pointsNum); centerDst.z = double(dstSumZ / pointsNum); Eigen::MatrixX3d srcMat; Eigen::MatrixX3d dstMat; for (int i = 0; i < pointsNum; ++i)//N组点 { //三行 srcMat(0, i) = srcPoints[i].x - centerSrc.x; srcMat(1, i) = srcPoints[i].y - centerSrc.y; srcMat(2, i) = srcPoints[i].z - centerSrc.z; dstMat(0, i) = dstPoints[i].x - centerDst.x; dstMat(1, i) = dstPoints[i].y - centerDst.y; dstMat(2, i) = dstPoints[i].z - centerDst.z; } Eigen::Matrix3d matS = srcMat * dstMat.transpose(); Eigen::JacobiSVD<Eigen::MatrixXf> svd(matS, Eigen::ComputeThinU | Eigen::ComputeThinV); Eigen::Matrix3d matV = svd.matrixV(), matU = svd.matrixU(); Eigen::Matrix3d matM; matM << 1, 0, 0, 0, 1, 0, 0, 0, matU*matV.determinant(); Eigen::Matrix3d matR = matV.transpose() * matM * matU.transpose(); double delta_X = centerDst.x - (centerSrc.x * matR(0, 0) + centerSrc.y * matR(0,1) + centerSrc.z * matR(0,2)); double delta_Y = centerDst.y - (centerSrc.x * matR(1, 0) + centerSrc.y * matR(1,1) + centerSrc.z * matR(1,2)); double delta_Z = centerDst.z - (centerSrc.x * matR(2, 0) + centerSrc.y * matR(2,1) + centerSrc.z * matR(2,2)); Eigen::Vector3d t(delta_X, delta_Y, delta_Z); T.rotate(matR); T.translate(t); return T; } 请写出这个源代码是什么意思
03-08
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值