基于有理逼近算法的序列密码的有理分数表示

题目描述

已知一条二元序列a(n)=(a0, a1, ..., an-1). 对1 ≤ k ≤ n, 求有限序列a(k)=(a0, a1, ..., ak-1)的有理分数表示. 序列a(n)请参见附件“sequence.txt”, 其中n = 1966000

设a(n)=(a0, a1, ..., an-1)是一条有限二元序列, 即ai ∈ {0,1},0 ≤ i ≤ n - 1. 若有理分数 p q满足q是正奇数, gcd(p,q) = 1, 并且 p≡ q(a0 +a1*2+⋯+an-1*pow(2, n-1)) mod pow(2,n), 则称 p q是序列a(n)的有理分数表示.  

符号说明

同余符号 ≡: 设n是一个正整数. 对任意的整数a和b, 有a≡ b mod n当且仅 当n整除a-b. 

两个整数a和b的最大公因子记为gcd(a, b). 

对两个整数a和b, 记Φ(a,b) = max {|a|,|b|}. 

题目分析

对于一个二元序列求出有理分数表示,其中关键因素在三点:

1.q需为正奇数。

2.p跟q不能存在公因子

3.[p-q(a0+a1*2+......+an-1*pow(2, n-1))]%pow(2, n)==0  即可被2的n次方整除

设x = a0+a1*2+......+an-1*pow(2, n-1);即 [p-qx]%pow(2,n)==0

有理逼近算法描述

其中1.a为不断累加的功能,也就是上述分析中x的值。

2.f=(f1, f2); g=(g1, g2)

f1, f2, g1, g2的值之间没有任何内在联系,只是这样写可以简化一下过程。

3.g2存在一个坑点,g2即为q的值,虽然算法描述中g2并没有给出取值范围,但题目中明确表示q为正奇数,也即是g2为正奇数

4.d的取值范围,当abs(g1)!=abs(g2)时,d在-(f1+f2)/(g1+g2)和(f1-f2)/(g1-g2)中取

否则,d的范围可设置为一个数据类型的最大值和最小值之间取。

有理逼近算法代码实现

// vs_project1.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include <iostream>
#include <algorithm>
#include <string>
#include <fstream>
#include <cmath>
#include <sstream>
#include <cassert>
#include <cstring>
#include <fstream>
using namespace std;

long long abs_max(long long a, long long b)
{
	return max(abs(a), abs(b));
}


int main()
{
	ios::sync_with_stdio(false);
	ifstream inf;
	inf.open("C://Users//49627//Desktop//sequence.txt");
	if (!inf)
		cout << "error" << endl;
	char c;
	int num = 0; //控制位数
	long double t = 0, f1, f2, g1, g2;
	long long mod;
	bool flag = false;
	while (!inf.eof() && num<45)
	{
		inf >> c;
		int tt = c - '0';
		if (tt == 0 && !flag)
		{
			num++;
			continue;
		}
		else if(!flag)
		{
			t = pow(2, num);
			f1 = 0;
			f2 = 2;
			g1 = t;
			g2 = 1;
			flag = true;
			num++;
			continue;
		}
		t += tt * pow(2, num);
		mod = pow(2, num + 1);
		long long sum = t*g2 - g1;
		long long ss = sum % mod;
		if (ss == 0)
		{
			cout << "the first case: " << num << endl;
			f1 *= 2;
			f2 *= 2;
		}
		else if (abs_max(g1, g2) < abs_max(f1, f2))
		{
			cout << "the second case: " << num << endl;
			long long mixn = 0x3f3f3f3f;
			long long x = 0, y = 0, a = 0, b = 0, d1 = -0x3f3f3f3f, d2 = 0x3f3f3f3f;
			if (abs(g1) != abs(g2))
			{
				d1 = -(f1 + f2) / (g1 + g2);
				d2 = (f1 - f2) / (g1 - g2);
				if (d1 > d2)
				{
					long long d = d1;
					d1 = d2;
					d2 = d;
				}
			}
			if ((d1 - 100) % 2 == 0)
				d1 = d1 - 99;
			else
				d1 = d1 - 100;
			for (long long j = d1; j <= d2 + 100; j+=2)
			{
				long long temp = abs_max(f1 + j*g1, f2 + j*g2);
				if (mixn > temp && (f2 + j*g2) > 0)
				{
					mixn = temp;
					x = f1 + j*g1;
					y = f2 + j*g2;
					a = g1;
					b = g2;
				}
			}
			g1 = x;
			g2 = y;
			f1 = 2 * a;
			f2 = 2 * b;
		}
		else
		{
			cout << "the third case: " << num << endl;
			long long mixn = 0x3f3f3f3f;
			long long x = 0, y = 0, d1 = -9999, d2 = 9999;
			if (abs(g1) != abs(g2))
			{
				d1 = -(f1 + f2) / (g1 + g2);
				d2 = (f1 - f2) / (g1 - g2);
				if (d1 > d2)
				{
					long long d = d1;
					d1 = d2;
					d2 = d;
				}
			}
			if ((d1 - 100) % 2 == 0)
				d1 = d1 - 99;
			else
				d1 = d1 - 100;
			for (long long j = d1; j <= d2 + 10; j+=2)
			{
				long long temp = abs_max(g1 + j*f1, g2 + j*f2);
				if (mixn > temp && (g2 + j*f2)>0)
				{
					mixn = temp;
					x = g1 + j*f1;
					y = g2 + j*f2;
				}
			}
			g1 = x;
			g2 = y;
			f1 *= 2;
			f2 *= 2;
		}
		num++;
	}
	cout << "f1:" << f1 << " " << "f2:" << f2 << " " << "g1:" << g1 << " " << "g2:" << g2 << endl;
	system("pause");
	inf.close();
	return 0;
}

 

**验证算法正确性

#include <iostream>
#include <algorithm>
#include <string>
#include <fstream>
#include <cmath>
#include <sstream>
#include <cstring>
#include <fstream>
#define MAX_L 2005
#define is "0100000000100111111110001000001111111111011100000000111"

using namespace std;

class bign
{
public:
    int len, s[MAX_L];//数的长度,记录数组
//构造函数
    bign();
    bign(const char*);
    bign(int);
    bool sign;//符号 1正数 0负数
    string toStr() const;//转化为字符串,主要是便于输出
    friend istream& operator>>(istream &,bign &);//重载输入流
    friend ostream& operator<<(ostream &,bign &);//重载输出流
//重载复制
    bign operator=(const char*);
    bign operator=(int);
    bign operator=(const string);
//重载各种比较
    bool operator>(const bign &) const;
    bool operator>=(const bign &) const;
    bool operator<(const bign &) const;
    bool operator<=(const bign &) const;
    bool operator==(const bign &) const;
    bool operator!=(const bign &) const;
//重载四则运算
    bign operator+(const bign &) const;
    bign operator++();
    bign operator++(int);
    bign operator+=(const bign&);
    bign operator-(const bign &) const;
    bign operator--();
    bign operator--(int);
    bign operator-=(const bign&);
    bign operator*(const bign &)const;
    bign operator*(const int num)const;
    bign operator*=(const bign&);
    bign operator/(const bign&)const;
    bign operator/=(const bign&);
//四则运算的衍生运算
    bign operator%(const bign&)const;//取模(余数)
    bign factorial()const;//阶乘
    bign Sqrt()const;//整数开根(向下取整)
    bign pow(const bign&)const;//次方
//辅助的函数
    void clean();
    ~bign();
};
#define max(a,b) a>b ? a : b
#define min(a,b) a<b ? a : b

bign::bign()
{
    memset(s, 0, sizeof(s));
    len = 1;
    sign = 1;
}

bign::bign(const char *num)
{
    *this = num;
}

bign::bign(int num)
{
    *this = num;
}

string bign::toStr() const
{
    string res;
    res = "";
    for (int i = 0; i < len; i++)
        res = (char)(s[i] + '0') + res;
    if (res == "")
        res = "0";
    if (!sign&&res != "0")
        res = "-" + res;
    return res;
}

istream &operator>>(istream &in, bign &num)
{
    string str;
    in>>str;
    num=str;
    return in;
}

ostream &operator<<(ostream &out, bign &num)
{
    out<<num.toStr();
    return out;
}

bign bign::operator=(const char *num)
{
    memset(s, 0, sizeof(s));
    char a[MAX_L] = "";
    if (num[0] != '-')
        strcpy(a, num);
    else
        for (int i = 1; i < strlen(num); i++)
            a[i - 1] = num[i];
    sign = !(num[0] == '-');
    len = strlen(a);
    for (int i = 0; i < strlen(a); i++)
        s[i] = a[len - i - 1] - 48;
    return *this;
}

bign bign::operator=(int num)
{
    if (num < 0)
        sign = 0, num = -num;
    else
        sign = 1;
    char temp[MAX_L];
    sprintf(temp, "%d", num);
    *this = temp;
    return *this;
}

bign bign::operator=(const string num)
{
    const char *tmp;
    tmp = num.c_str();
    *this = tmp;
    return *this;
}

bool bign::operator<(const bign &num) const
{
    if (sign^num.sign)
        return num.sign;
    if (len != num.len)
        return len < num.len;
    for (int i = len - 1; i >= 0; i--)
        if (s[i] != num.s[i])
            return sign ? (s[i] < num.s[i]) : (!(s[i] < num.s[i]));
    return !sign;
}

bool bign::operator>(const bign&num)const
{
    return num < *this;
}

bool bign::operator<=(const bign&num)const
{
    return !(*this>num);
}

bool bign::operator>=(const bign&num)const
{
    return !(*this<num);
}

bool bign::operator!=(const bign&num)const
{
    return *this > num || *this < num;
}

bool bign::operator==(const bign&num)const
{
    return !(num != *this);
}

bign bign::operator+(const bign &num) const
{
    if (sign^num.sign)
    {
        bign tmp = sign ? num : *this;
        tmp.sign = 1;
        return sign ? *this - tmp : num - tmp;
    }
    bign result;
    result.len = 0;
    int temp = 0;
    for (int i = 0; temp || i < (max(len, num.len)); i++)
    {
        int t = s[i] + num.s[i] + temp;
        result.s[result.len++] = t % 10;
        temp = t / 10;
    }
    result.sign = sign;
    return result;
}

bign bign::operator++()
{
    *this = *this + 1;
    return *this;
}

bign bign::operator++(int)
{
    bign old = *this;
    ++(*this);
    return old;
}

bign bign::operator+=(const bign &num)
{
    *this = *this + num;
    return *this;
}

bign bign::operator-(const bign &num) const
{
    bign b=num,a=*this;
    if (!num.sign && !sign)
    {
        b.sign=1;
        a.sign=1;
        return b-a;
    }
    if (!b.sign)
    {
        b.sign=1;
        return a+b;
    }
    if (!a.sign)
    {
        a.sign=1;
        b=bign(0)-(a+b);
        return b;
    }
    if (a<b)
    {
        bign c=(b-a);
        c.sign=false;
        return c;
    }
    bign result;
    result.len = 0;
    for (int i = 0, g = 0; i < a.len; i++)
    {
        int x = a.s[i] - g;
        if (i < b.len) x -= b.s[i];
        if (x >= 0) g = 0;
        else
        {
            g = 1;
            x += 10;
        }
        result.s[result.len++] = x;
    }
    result.clean();
    return result;
}

bign bign::operator * (const bign &num)const
{
    bign result;
    result.len = len + num.len;

    for (int i = 0; i < len; i++)
        for (int j = 0; j < num.len; j++)
            result.s[i + j] += s[i] * num.s[j];

    for (int i = 0; i < result.len; i++)
    {
        result.s[i + 1] += result.s[i] / 10;
        result.s[i] %= 10;
    }
    result.clean();
    result.sign = !(sign^num.sign);
    return result;
}

bign bign::operator*(const int num)const
{
    bign x = num;
    bign z = *this;
    return x*z;
}
bign bign::operator*=(const bign&num)
{
    *this = *this * num;
    return *this;
}

bign bign::operator /(const bign&num)const
{
    bign ans;
    ans.len = len - num.len + 1;
    if (ans.len < 0)
    {
        ans.len = 1;
        return ans;
    }

    bign divisor = *this, divid = num;
    divisor.sign = divid.sign = 1;
    int k = ans.len - 1;
    int j = len - 1;
    while (k >= 0)
    {
        while (divisor.s[j] == 0) j--;
        if (k > j) k = j;
        char z[MAX_L];
        memset(z, 0, sizeof(z));
        for (int i = j; i >= k; i--)
            z[j - i] = divisor.s[i] + '0';
        bign dividend = z;
        if (dividend < divid) { k--; continue; }
        int key = 0;
        while (divid*key <= dividend) key++;
        key--;
        ans.s[k] = key;
        bign temp = divid*key;
        for (int i = 0; i < k; i++)
            temp = temp * 10;
        divisor = divisor - temp;
        k--;
    }
    ans.clean();
    ans.sign = !(sign^num.sign);
    return ans;
}

bign bign::operator/=(const bign&num)
{
    *this = *this / num;
    return *this;
}

bign bign::operator%(const bign& num)const
{
    bign a = *this, b = num;
    a.sign = b.sign = 1;
    bign result, temp = a / b*b;
    result = a - temp;
    result.sign = sign;
    return result;
}

bign bign::pow(const bign& num)const
{
    bign result = 1;
    for (bign i = 0; i < num; i++)
        result = result*(*this);
    return result;
}

bign bign::factorial()const
{
    bign result = 1;
    for (bign i = 1; i <= *this; i++)
        result *= i;
    return result;
}

void bign::clean()
{
    if (len == 0) len++;
    while (len > 1 && s[len - 1] == '\0')
        len--;
}

bign bign::Sqrt()const
{
    if(*this<0)return -1;
    if(*this<=1)return *this;
    bign l=0,r=*this,mid;
    while(r-l>1)
    {
        mid=(l+r)/2;
        if(mid*mid>*this)
            r=mid;
        else
            l=mid;
    }
    return l;
}

bign::~bign()
{
}


///求出序列长度
/*
int main()
{
    ios::sync_with_stdio(false);
    ifstream myfile("C:\\Users\\49627\\Desktop\\全国高校密码竞赛\\赛题一\\sequence.txt");
    string temp;
    if (!myfile.is_open())
    {
        cout << "open file error!" << endl;
    }
    int sum = 0;
    while(getline(myfile,temp))
    {
        int len = temp.length();
        sum += len;
    }
    cout<<sum<<endl;
    myfile.close();
    return 0;
}

///1966000

*/
bign mixn(bign a, bign b)
{
    if(a<0)
        a *= -1;
    if(a>b)
        return a;
    else
        return b;
}

bool judge(bign a, bign b)
{
    bign c=a<b?a:b;
    for(bign i=2; i<=c; ++i)
        if(a%i==0 && b%i==0)
            return false;
    return true;
}


int main()
{
    ios::sync_with_stdio(false);
    int n=0;
    ofstream outf;
    outf.open("output.txt");
    //outf<<"n\tp\tq\tx\tpow(2,n)"<<endl;
    while(n<200)
    {
        bign a=0, b;
        for(int i=0; i<n; ++i)
        {
            bool tt = is[i] - '0';
            a += tt * pow(2, i);
        }
        b = pow(2, n);
        bign p, q, abm=0x3f3f3f3f;
        bign i=0, j=0;
        for(p=-9999999; p<=9999999; ++p)
        {
            cout<<p<<endl;
            for(q=3;q<=9999999;++q)
            {
                if((p-q*a)%b==0 && judge(q, p))
                {
                    bign x = mixn(p, q);
                    if(abm > x)
                    {
                        abm = x;
                        i=p;
                        j=q;
                        cout<<abm<<" "<<p<<" "<<q<<endl;
                    }
                }
            }
        }
        outf<<n<<"\t"<<i<<"\t"<<j<<"\t"<<a<<"\t"<<b<<endl;
        cout<<n<<endl;
        n++;
    }
    outf.close();
    return 0;
}

 

版权声明:本文为AcSuccess原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/AcSuccess/article/details/80587514

智能推荐

ActiveMQ学习4-ActiveMQ的安全机制和集群模式

ActiveMQ的安全机制和集群模式 20 ActiveMQ安全机制 20.1 Web 控制台安全 20.2 消息服务器Broker安全 21 ActiveMQ主从集群 21.1 使用集群的重要性 20.2 主从集群的方式 20.2.1 shared filesystem Master-Slave方式主从集群 20.2.2 shared database Master-Slave方式主从集群 20...

说说 Python Django 应用的基础目录结构

通过以下 django-admin 指令创建应用之后,就会生成应用的基础目录结构。 比如,我们建立了一个叫 ‘first’ 的应用,它的目录结构是这样的: 目录或文件 说明 最外层的 first/ 这是新应用的根目录,所有与该应用相关的内容都放在这里。 manage.py 用于管理 Django 项目的命令行工具。 里面一层的 first/ 目录 是一个...

Springboot整合rabbitMQ

依赖: 配置文件application.yml RabbitConfig 消息生产者RabbitProducer 消息消费者RabbitCustomer 通过Controller进行调用 启动项目后调用接口: 结果:...

Thread.join()方法的使用

如果一个线程A执行了thread.join()语句,代表当前线程A等待thread线程终止后才从thread.join()方法返回 并且这个方法具有超时特性,可以添加参数设置 输出结果: jdk中Thread.join()方法的源码(进行了部门调整)   每个线程终止的条件是前驱线程的终止,每个线程等待前驱线程终止后,才从join()方法返回,  当线程终止时,会调用自身的no...

linux服务器部署jenkins笔记

安装jenkins参考文档:https://blog.csdn.net/tomatocc/article/details/83930714 1. 打开jenkins官网:https://jenkins.io/download/ 将war包下载到本地 **ps:**这里要注意的是要下载左边下方的war包,不要下载右边下面的war包。左边是稳定版本,右边是最新版本,建议大家使用稳定版本(我刚开始下载的...

猜你喜欢

k8s部署elasticsearch集群

百度营销大学     环境准备 我们使用的k8s和ceph环境见: https://blog.51cto.com/leejia/2495558 https://blog.51cto.com/leejia/2499684 ECK简介 Elastic Cloud on Kubernetes,这是一款基于 Kubernetes Operator 模式的新型编排产品,用户可使用该产品在...

saas-export项目-AdminLTE介绍与入门

AdminLTE介绍 (1)AdminLTE是什么? AdminLTE是一款建立在bootstrap和jquery之上的开源的模板主题工具 (2)AdminLTE有什么特点? 提供一系列响应的、可重复使用的组件, 并内置了多个模板页面 自适应多种屏幕分辨率,兼容PC和移动端 快速的创建一个响应式的Html5网站 AdminLTE 不但美观, 而且可以免去写很大CSS与JS的工作量 AdminLTE...

MyBatis中ResultMap结果集映射

用于解决属性名和字段名不一致的情况: resultMap 元素是 MyBatis 中最重要最强大的元素。...

编写一个shell

编写shell的过程: 1.从标准输入中读入一个字符串。 2.解析字符串 3.创建一个子进程的执行程序。 4.子进程程序替换。 5.父进程等待子进程退出。...

WEB自动化测试中Xpath定位方法

前言: Xpath是在XML文档中查找信息的一种语言,使用路径表达式来选取XML文档中的节点或节点集,由于XML与HTML结构类似(前者用于传输数据,后者用于显示数据),所以Xpath也常用于查找HTML文档中的节点或节点集。 一  路径表达式: 路径以“/”开始     表示找到满足该绝对路径的元素; 路径以//”开始  ...