描述 Description
给定一个整数n，求一个整数m，满足m<=n，并且m/phi(m)的值最大。
注：phi(m)代表m的欧拉函数，即不大于m且与m互质的数的个数。

题解：
m/phi（m) 很容易化成 连积（p/（p-1））  p|m
所以就很简单了，将最小的质数乘起来，直到>n，输出前一个。
因为保证最小所以只乘一次，因为p/（p-1）单调减，所以从小的开始选。
高精度写错搞了好久，然后有卡了几次时才过了
代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 100000

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 10

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }

 int n,m,tot,p[maxn];

 bool v[maxn];

 char ch[maxn];

 class bigg{

 public:

     int num[maxn],len;

     bigg()

     {

         memset(num,,sizeof(num));

         len=;

     }

     inline bigg operator =(const bigg &b)

     {

         memset(num,,sizeof(num));

         len=b.len;

         for1(i,len)num[i]=b.num[i];

         return(*this);

     }

     inline bigg operator =(int b)

     {

         memset(num,,sizeof(num));

         len=;

         while(b){num[++len]=b%mod;b/=mod;}

         return(*this);

     }

     inline bigg operator *(int b)

     {

         for1(i,len)num[i]*=b;

         for1(i,len)

         {

             num[i+]+=num[i]/mod;

             num[i]%=mod;

             if(num[len+])len++;

         }

         return(*this);

     }

     inline bool operator <(const bigg&b)

    {

       if(len!=b.len)return len<b.len;

       for3(i,len,)if(num[i]!=b.num[i])return num[i]<b.num[i];

       return ;

    }

     inline void print()

     {

         printf("%d",num[len]);

         for3(i,len-,)printf("%d",num[i]);printf("\n");

     }

 };

 bigg a[],b,c[];

 int rk[];

 bool cmp(int x,int y){return a[x]<a[y];}

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     for2(i,,maxn)

     {

         if(!v[i])p[++tot]=i;

         for1(j,tot)

         {

             int k=p[j]*i;

             if(k>maxn)break;

             v[k]=;

             if(i%p[j]==)break;

         }

     }

     n=read();

     for1(i,n)

     {

         scanf("%s",ch+);

         a[i].len=strlen(ch+);

         for1(j,a[i].len)a[i].num[j]=ch[a[i].len+-j]-'';

         rk[i]=i;

     }

     sort(rk+,rk+n+,cmp);

     int j=;

     b=j;

     //for1(i,100)b=b*p[i],b.print();

     for1(i,n)

     {

         c[rk[i]]=c[rk[i-]];

         while(b<a[rk[i]])c[rk[i]]=b,b=b*p[j++];

         //b.print();

     }

     //for1(i,n)a[rk[i]].print(),c[rk[i]].print();

     for1(i,n)c[i].print();

     return ;

 }

不知道每组询问暴力求会不会T，我为了保险拍了个序233

UPD：这个程序交到bz也T了。。。TAT

无奈看了lyd的程序，居然预处理压了8位！

orzz

憋了1h+终于写出来了，二分+预处理。。。

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 100000

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 100000000

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }

 int n,m,tot,p[maxn];

 bool v[maxn];

 char s[maxn];

 class bigg{

 public:

     int num[],len;

     bigg()

     {

         memset(num,,sizeof(num));

         len=;

     }

     inline bigg operator =(const bigg &b)

     {

         memset(num,,sizeof(num));

         len=b.len;

         for1(i,len)num[i]=b.num[i];

         return(*this);

     }

     inline bigg operator =(int b)

     {

         memset(num,,sizeof(num));

         len=;

         while(b){num[++len]=b%mod;b/=mod;}

         return(*this);

     }

     inline bigg operator *(int b)

     {

        ll x=;

         for1(i,len)

         {

             x+=(ll)num[i]*b;

             num[i]=x%mod;

             x/=mod;

         }

         if(x)num[++len]=x;

         return(*this);

     }

     inline bool operator <(const bigg&b)

    {

       if(len!=b.len)return len<b.len;

       for3(i,len,)if(num[i]!=b.num[i])return num[i]<b.num[i];

       return ;

    }

     inline void print()

     {

         printf("%d",num[len]);

         for3(i,len-,)printf("%08d",num[i]);printf("\n");

     }

 };

 bigg a,b[];

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     for2(i,,maxn)

     {

         if(!v[i])p[++tot]=i;

         for1(j,tot)

         {

             int k=p[j]*i;

             if(k>maxn)break;

             v[k]=;

             if(i%p[j]==)break;

         }

     }

     b[]=;

     for1(i,)b[i]=b[i-],b[i]=b[i]*p[i];

     int cs=read();

     while(cs--)

     {

         memset(s,,sizeof(s));

         scanf("%s",s);

         n=strlen(s);

         reverse(s,s+n);

         for0(i,n-)s[i]-='';

         a.len=(n+)/;

         for0(i,a.len-)

         a.num[i+]=s[i*]+*s[i*+]+*s[i*+]+*s[i*+]+*(s[i*+]+*s[i*+]+*s[i*+]+*s[i*+]);

         int l=,r=,mid;

         while(l<=r)

         {

             mid=(l+r)>>;

             if(b[mid]<a)l=mid+;else r=mid-;

         }

         b[r].print();

     }

     return ;

 }

UPD：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 100000

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 100000000

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }

 int n,m,tot,p[maxn];

 bool v[maxn];

 char s[maxn];

 class bigg{

 public:

     int num[],len;

     bigg()

     {

         memset(num,,sizeof(num));

         len=;

     }

     inline bigg operator =(const bigg &b)

     {

         memset(num,,sizeof(num));

         len=b.len;

         for1(i,len)num[i]=b.num[i];

         return(*this);

     }

     inline bigg operator =(int b)

     {

         memset(num,,sizeof(num));

         len=;

         while(b){num[++len]=b%mod;b/=mod;}

         return(*this);

     }

     inline bigg operator *(int b)

     {

        ll x=;

         for1(i,len)

         {

             x+=(ll)num[i]*b;

             num[i]=x%mod;

             x/=mod;

         }

         if(x)num[++len]=x;

         return(*this);

     }

     inline bool operator <(const bigg&b)

    {

       if(len!=b.len)return len<b.len;

       for3(i,len,)if(num[i]!=b.num[i])return num[i]<b.num[i];

       return ;

    }

     inline void print()

     {

         printf("%d",num[len]);

         for3(i,len-,)printf("%08d",num[i]);printf("\n");

     }

 };

 bigg a[],b,c[];

 int rk[];

 bool cmp(int x,int y){return a[x]<a[y];}

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     for2(i,,maxn)

     {

         if(!v[i])p[++tot]=i;

         for1(j,tot)

         {

             int k=p[j]*i;

             if(k>maxn)break;

             v[k]=;

             if(i%p[j]==)break;

         }

     }

     n=read();

     for1(j,n)

     {

         memset(s,,sizeof(s));

         scanf("%s",s);

         m=strlen(s);

         reverse(s,s+m);

         for0(i,m-)s[i]-='';

         a[j].len=(m+)/;

         for0(i,a[j].len-)

         a[j].num[i+]=s[i*]+*s[i*+]+*s[i*+]+*s[i*+]+*(s[i*+]+*s[i*+]+*s[i*+]+*s[i*+]);

         rk[j]=j;

     }

     sort(rk+,rk+n+,cmp);

     int j=;

     b=;

     for1(i,n)

     {

         c[rk[i]]=c[rk[i-]];

         while(b<a[rk[i]])c[rk[i]]=b,b=b*p[j++];

     }

     for1(i,n)c[i].print();

     return ;

 }

哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈

用我的方法+压8位 怒排rank3

代码：

「Poetize3」Heaven Cow与God Bull的更多相关文章

  1. bzoj3034: Heaven Cow与God Bull

    Description __int64 ago,there's a heaven cow called sjy...A god bull named wzc fell in love with her ...

  2. 【Heaven Cow与God Bull】题解

    题目 Description __int64 ago,there's a heaven cow called sjy... A god bull named wzc fell in love with ...

  3. LG5200 「USACO2019JAN」Sleepy Cow Sorting 树状数组

    \(\mathrm{Sleepy Cow Sorting}\) 问题描述 LG5200 题解 树状数组. 设\(c[i]\)代表\([1,i]\)中归位数. 显然最终的目的是将整个序列排序为一个上升序 ...

  4. 「Poetize3」导弹防御塔

    描述 Description Freda控制着N座可以发射导弹的防御塔.每座塔都有足够数量的导弹,但是每座塔每次只能发射一枚.在发射导弹时,导弹需要T1秒才能从防御塔中射出,而在发射导弹后,发射这枚导 ...

  5. 「Poetize3」绿豆蛙的归宿

    描述 Description 给出一个有向无环图,起点为1终点为N,每条边都有一个长度,并且从起点出发能够到达所有的点,所有的点也都能够到达终点.绿豆蛙从起点出发,走向终点.到达每一个顶点时,如果有K ...

  6. joyoi1935 「Poetize3」导弹防御塔

    #include <iostream> #include <cstring> #include <cstdio> #include <queue> #i ...

  7. 「POJ3613」Cow Relays

    「POJ3613」Cow Relays 传送门 就一个思想:\(N\) 遍 \(\text{Floyd}\) 求出经过 \(N\) 个点的最短路 看一眼数据范围,想到离散化+矩阵快速幂 代码: #in ...

  8. 「笔记」折半搜索（Meet in the Middle）

    思想 先搜索前一半的状态,再搜索后一半的状态,再记录两边状态相结合的答案. 暴力搜索的时间复杂度通常是 \(O(2^{n})\) 级别的.但折半搜索可以将时间复杂度降到 \(O(2 \times 2^ ...

  9. 「译」JUnit 5 系列：条件测试

    原文地址:http://blog.codefx.org/libraries/junit-5-conditions/ 原文日期:08, May, 2016 译文首发:Linesh 的博客:「译」JUni ...

随机推荐

  1. heap堆内存不足

    配置jvm内存,加个启动参数 -Xms512m -Xmx512m -XX:PermSize=128m  -XX:MaxPermSize=256m   start 下面是这个: -Dcatalina.h ...

  2. asp.net 事件模型

    asp.net的原始设计构想,就是要让开发人员能够像 VB 开发工具那样,可以使用事件驱动式程序开发模式 (Event-Driven Programming Model) 的方法来开发网页与应用程序, ...

  3. hadoop集群环境搭建准备工作

    一定要注意hadoop和linux系统的位数一定要相同,就是说如果hadoop是32位的,linux系统也一定要安装32位的. 准备工作: 1 首先在VMware中建立6台虚拟机(配置默认即可).这是 ...

  4. gulp的常用api

    gulp是什么? http://gulpjs.com/ 相信你会明白的! 与著名的构建工具grunt相比,有什么优势呢? 易于使用,代码优于配置 高效,不会产生过多的中间文件,减少I/O压力 易于学习 ...

  5. SQL:42601

    以前遇到SQL中触发器的问题, 添加一个字段后,发现以前的程序的SQL报错 看了下表定义,有这样的一行代码 CREATE TRIGGER chintai_keiyaku_audit AFTER INS ...

  6. java构造函数也可以用private开头

    private 构造函数一般用于Singleton模式,指的是整个应用只有本类的一个对象,一般这种类都有一个类似getInstance()的方法!下面是一个Singleton的例子:public cl ...

  7. VS2015升级Update2之后Cordova程序提示:此应用程序无法在此电脑上运行

    VS2015在升级到Update2之后,有可能出现如下异常,在运行Cordova项目时提示: 查看输出面板会有乱码错误信息: 出现此问题的原因是在于npm程序损坏了.vs调用的npm程序并不是在nod ...

  8. jquery,javascript -设置某一ul下的li下的 a的属性

    //javascriptvar ul = document.getElementById('ul); var as = ul.getElementsByTagName('a'); for(var i ...

  9. 文字排版--字体(font-family)

    我们可以使用css样式为网页中的文字设置字体.字号.颜色等样式属性.下面我们来看一个例子,下面代码实现:为网页中的文字设置字体为宋体. body{font-family:"宋体"; ...

  10. EA创建用例图步骤详解

    EA创建用例图步骤详解 1 创建一个项目 2 选择需要的模型 3 新建模型包 4 新建图表 5 新建模型包 6 创建用户角色Actor 7 新建用例 8 关联用户和用例 9 最后整个项目浏览器目录结构 ...