poj3708:函数式化简+高精度进制转换+同余方程组
题目大意
给定一个函数
找出满足条件
等于 k
的最小的x
m,k,d已知 其中 m,k 很大需要使用高精度存储
思路:
对 函数f(m)进行化简 ,令t=ceil( log(d,m) )
可以得到 f(m)=d ^ t * ( a [ m / (d^t) ] ) + d ^ (t-1) * ( b[ m/( d^(t-1) ) ] )......+b [m%d] ;
我们一看,每一项都是 跟 d 的次方有关,所以考虑使用 d 进制进行计算
设 m=a1b1b2b3b4(d进制)
那么 f(m)=a[a1]b[b2]b[b3]b[b4](d进制)
以此类推
fx(m)=k,即 m通过 x次上述变换后得到了k (这里的 m 和 k 均为 d 进制)
于是 我们可以找到 m 每一位第一次等于k的情况以及这一位的循环节 分别存在数组中
那么问题就转化为了解同余方程
这个题还有很多无解的情况,需要注意特判
代码:
#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
using namespace std;
#define MAXN 100000
char m[];
char k[];
long long p[];
long long q[];
long long start[];
long long ans[];
long long res[];
long long mm[];
long long kk[];
long long a[];
long long b[];
long long d;
long long vi[];
long long getnum(char c)
{
return c-'';
}
void trans(char* str,long long base0,long long base1)
{
memset(res,,sizeof(res));
long long y,i,j;
start[]=strlen(str);
for(i=;i<=start[];i++)
{
start[i]=getnum(str[i-]);
}
while(start[]>=)
{
y=; //余数
ans[]=start[];
for(i=;i<=start[];i++)
{
y=y*base0+start[i];
ans[i]=y/base1;
y%=base1;
}
res[++res[]]=y; //这一轮的余数
i=;
while(i<=ans[]&&ans[i]==)
i++;
memset(start,,sizeof(start));
for(j=i;j<=ans[];j++)
start[++start[]]=ans[j];
memset(ans,,sizeof(ans));
}
return;
}
long long save(long long *dis)
{
for(long long i=res[];i;i--)
{
dis[res[]-i]=res[i];
}
return res[];
}
long long exgcd(long long a,long long b,long long &x,long long &y)
{
if(!b)
{
x=;
y=;
return a;
}
long long tt=exgcd(b,a%b,x,y);
long long t;
t=x;
x=y;
y=(t-a/b*y);
return tt;
}
long long solve(long long r)
{
long long a1,a2,b1,b2,x,y,A,B,C,d,t;
a1=a[];
b1=b[];
for(long long i=;i<r;i++)
{
a2=a[i];
b2=b[i];
A=a1;
B=a2;
C=b2-b1;
d=exgcd(A,B,x,y);
if(C%d)
{
return -;
}
t=B/d;
x=(x*(C/d)%t+t)%t;
b1=a1*x+b1;
a1=a1/d*a2;
}
return b1;
}
int main()
{
while(scanf("%lld",&d),d!=-)
{
for(long long i=;i<d;i++)
scanf("%lld",p+i);
for(long long i=;i<d;i++)
scanf("%lld",q+i);
scanf("%s%s",m,k);
trans(m,,d);
long long nm=save(mm);
trans(k,,d);
long long nk=save(kk);
if(nm!=nk)
{
puts("NO");
continue;
}
memset(a,-,sizeof(a));
memset(b,-,sizeof(b));
bool have=;
for(long long i=;i<nm;i++)
{
bool ok=;
long long tmp=;
memset(vi,,sizeof(vi));
for(long long j=mm[i];;)
{
if(j==kk[i])
{
if(ok==)
{
b[i]=tmp;
ok=;
}
else
{
a[i]=tmp-b[i];
break;
}
}
tmp++;
if(i==)
j=p[j];
else
j=q[j];
if(tmp>=)
{
ok=;
break;
}
}
if(!ok)
{
have=;
break;
}
}
if(!have)
{
puts("NO");
continue;
}
long long ans=solve(nm);
if(ans==-)
{
puts("NO");
continue;
}
printf("%lld\n",ans);
}
return ;
}
poj3708:函数式化简+高精度进制转换+同余方程组的更多相关文章
- poj1220:高精度进制转换模板题
今天撸3708 一直奇怪的re 就先放下了,写这个题的过程中学习了一个高精度进制转换,用这个模板写了1220 记录一下: #include <iostream> #include < ...
- 高精度进制转换(poj1220)
常规短除法原理 高精度进制转换是对于特别大的数字来说的,当数字特别大时,难以进行除法和取余的操作,此时通过字符串模拟的办法可以解决. #include <iostream> #includ ...
- java高精度进制转换
POJ1131 由于本题只有小数部分(整数部分均为0),故在进制转换的之后只能自己手写转换方法了. 8进制转换10进制的方法为,以0.75为例,应是7*8^-1 + 5*8^-2.所以呢,可以 ...
- POJ 1220 高精度/进制转换
n进制转m进制,虽然知道短除法但是还是不太理解,看了代码理解一些了: 记住这个就好了: for(int k=0;l; ){ for(int i=l ; i>=1 ; i--){ num[i - ...
- 高精度的N进制转换模板(转K神)
/* 高精度进制转换 把oldBase 进制的数转化为newBase 进制的数输出. 调用方法,输入str, oldBase newBase. change(); solve(); output(); ...
- poj1220 (高精度任意进制转换)
http://poj.org/problem?id=1220 高精度任意进制转换 代码是从discuss里找到的,据说是maigo神牛写的. 超精简!! 我自己第一写的时候,还把n进制先转成10进制, ...
- zznu 1068: 进制转换
进制应该属于程序员的看家本事了,也是大家水平告别菜鸟的一个转折,所以进制转换题目是很有意义的, 这个题目是最简单的把二进制数化简成十进制,因为输入有可能有31位,所以无法使用int或者long lon ...
- python 实现进制转换(二进制转十进制)
摘自https://baike.baidu.com/item/%E5%8D%81%E8%BF%9B%E5%88%B6%E8%BD%AC%E4%BA%8C%E8%BF%9B%E5%88%B6 pytho ...
- 华为笔试——C++进制转换
题目:2-62进制转换 题目介绍:输入一个n1 进制的整数(包括负数),将其转换成n2 进制,其中n1 .n2 的范围是 [ 2,62 ] .每个数字的范围是0-9.a-z.A-Z.不用考虑非法输入. ...
随机推荐
- MyBatis的事务处理
先来假设这样一个问题:如果数据库里面有一个用户表和一个作家表,那么当要添加一条数据到作家表中时,作家表的id必须是用户表中的其中一个id,因为作家一定也要是一个用户.这时就涉及到事务处理. 在上一篇博 ...
- 【winform程序】自定义webrowser控件调用IE的版本
修改注册表: bit: HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Internet Explorer\MAIN\FeatureControl\FEATURE_BROW ...
- 【性能优化】优化笔记之一:图像RGB与YUV转换优化
本文主要介绍如何优化您自己的CODE,实现软件的加速.我们一个图象模式识别的项目,需要将RGB格式的彩色图像先转换成黑白图像.图像转换的公式如下: Y = 0.299 * R + 0.587 * G ...
- [React] Extracting Private React Components
we leverage private components to break our render function into more manageable pieces without leak ...
- HDU 1015 Safecracker 解决问题的方法
Problem Description === Op tech briefing, 2002/11/02 06:42 CST === "The item is locked in a Kl ...
- 在VMware中为Linux系统安装VM-Tools的详解教程
在安装Linux的虚拟机中,单击“虚拟机”菜单下的“安装Vmware-Tools”. 先介绍一下下面安装该工具时要用到的几个目录: /mnt 挂载目录,用来临时挂载别的文件系统,硬件设备 /tmp临时 ...
- Xshell 4的上传与下载
http://blog.csdn.net/lioncode/article/details/7921525
- SQLServer2008/2012 删除所有表视图存储过程
SQLServer2008/2012 删除所有表视图存储过程 -------------------删除所有的表-------------------use xuwenbin111--/第1步**** ...
- C#程序调用cmd执行命令-MySql备份还原
1.简单实例 //备份还原mysql public static void TestOne() { Process p = new Process(); p.StartInfo.FileName = ...
- 包装 request Demo
//包装request,增强getParameter方法 class MyReq extends HttpServletRequestWrapper{ private HttpServletReque ...