最大MOD值

最大MOD值

题意:有一个长度为n的数列a,求数列中 a i % a j ( a i &gt; a j ) a_i\%a_j(a_i&gt;a_j) ai%aj(ai>aj)的最大值。n<=2e5 a i &lt; = 1 e 6 a_i&lt;=1e6 ai<=1e6

例:
n=3
a: 3 4 5
m a x ( a i % a j ) = 2 max(a_i\%a_j)=2 max(ai%aj)=2

原题

思路:
对于取模运算 b % a = c b\%a=c b%a=c ,我们可以把 b 看 成 b = a k + c b看成b=ak+c bb=ak+c 现在问题则是要使c最大。
若一个数 x = a k + d , x % a &gt; b % a x=ak+d,x\%a&gt;b\%a x=ak+dx%a>b%a 则x比b更加接近 a ( k + 1 ) a(k+1) a(k+1),换句话说:则是对于模数a
可将所有数字分在区间 [ 0 , a − 1 ] , [ a , 2 a − 1 ] , [ 2 a , 3 a − 1 ] . . . [ k a , ( k + 1 ) a − 1 ] [0,a-1],[a,2a-1],[2a,3a-1]...[ka,(k+1)a-1] [0,a1],[a,2a1],[2a,3a1]...[ka,(k+1)a1]的区间里。再对于每一个
区间,我们都采用二分找到属于该区间的最大的数x,那么这一段区间的数中模a最大的就是x。
如此对每一个数这样计算一遍,可以发现复杂度是可过的。(调和级数的复杂度)

/*******************************
Author:galaxy yr
LANG:C++
Created Time:2019年06月10日 星期一 22时09分07秒
*******************************/
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=2e5+10;
int n,a[maxn],ans;
int main()
{
        cin>>n;
        for(int i=1;i<=n;i++)
            cin>>a[i];
        sort(a+1,a+n+1);
        for(int i=1;i<=n;i++)
            for(int l=a[i],r=2*l-1;l<=a[n];l=r+1,r+=a[i])
                ans=max(ans,*(lower_bound(a+i,a+n+1,r+1)-1)%a[i]);
        cout<<ans<<endl;
        return 0;
}

推荐一波个人博客

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
最大matlab代码clear clc total=1;%总的进化代数 Num=0; MaN1=0; MaN2=0; flag=0; N1=initialize(); N2=initialize(); N3=initialize(); N4=initialize(); %复制-reproduction while (1) while (1) DecN1=BinToDec(N1);%计算适配度 DecN2=BinToDec(N2); DecN3=BinToDec(N3); DecN4=BinToDec(N4); f1=DecN1^2; f2=DecN2^2; f3=DecN3^2; f4=DecN4^2; ttt=[f1,f2,f3,f4]; if flag==0 MaN1=max(ttt); else MaN2=max(ttt); end if MaN1==MaN2 Num=Num+1; else Num=0; end sum=f1+f2+f3+f4; aver=sum/4; n1=round(f1/aver); %实际得到的复制个数 n2=round(f2/aver); n3=round(f3/aver); n4=round(f4/aver); N=n1+n2+n3+n4; %总的复制个数 if N~=4 %引入N:保持原种群大小不变 N1=initialize(); N2=initialize(); N3=initialize(); N4=initialize(); total=1; Num=0; MaN1=0; MaN2=0; flag=0; continue; else break; end end temp1=cal(N1,n1); temp2=cal(N2,n2); temp3=cal(N3,n3); temp4=cal(N4,n4); M=[temp1;temp2;temp3;temp4]; N1=M(1,:); N2=M(2,:); N3=M(3,:); N4=M(4,:); %交叉-crossover while 1 p=floor(1000*rand); if p>1 break; else continue; end end k=mod(p,3); switch k case 0 [N1,N2,N3,N4,ps]=crossover(N1,N2,N3,N4); case 1 [N1,N3,N2,N4,ps]=crossover(N1,N3,N2,N4); case 2 [N1,N4,N2,N3,ps]=crossover(N1,N4,N2,N3); end %变异-mutation U=[N1,N2,N3,N4]; pos=mod(floor(1000*rand),20)+1; if U(1,pos)==0 U(1,pos)=1; else U(1,pos)=0; end N1=U(1,1:5); N2=U(1,6:10); N3=U(1,11:15); N4=U(1,16:20); %遗传算法结束条件:连续10代最大均保持一致 if Num==10 disp('进化代数:') total-10 disp('现在的种群:') N1 N2 N3 N4 disp('最大:') MaN1 break; end flag=mod(flag+1,2); total=total+1; end

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值