Codeforces 70C Lucky Tickets（数学+机智）

最新推荐文章于 2024-02-07 07:00:00 发布

只玩三国的程序猿

最新推荐文章于 2024-02-07 07:00:00 发布

阅读量788

点赞数

分类专栏： 2016寒假集训文章标签： codeforces

本文链接：https://blog.csdn.net/CQUWEL/article/details/50620811

版权

2016寒假集训专栏收录该内容

24 篇文章 0 订阅

订阅专栏

题意：

给你一个maxx和maxy，让你选择一个上限x和y是的对于每一个a<=x，
在b<=y的范围下，满足a*b==rev(a)*rev(b)>=w的个数总和最大，并且
x*y最小（具体的看题目意思吧）。
这个题首先要知道一个性质，那就是x越大，符合条件的越多，y也是
一样的，也就是这个是递增的。首先，我们把等式改变一下，
变成a/rev(a)==rev(b)/b，那么就可以单独的考虑一个数了。之后
我们先固定y使它等于maxy，然后来增加x使得慢慢的可以符合条件。
然后再减少y，然后再增加x是的满足条件，这样就可以得到所有的情况，
取最优的就行。但是对于每一个x和y如何快速的计算呢？我们可以通过
一个pair，把a/rev(a)除以gcd(a,rev(a))，变成最简分式,
这样就利于我们来做统计。具体看代码吧。

代码：

//
//  Created by  CQU_CST_WuErli
//  Copyright (c) 2015 CQU_CST_WuErli. All rights reserved.
//
// #include<bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <sstream>
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define For_UVa if (kase!=1) cout << endl
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
#define BigInteger bign
template <typename T> T gcd(const T& a,const T& b) {return b==0?a:gcd(b,a%b);}
const int MAX_L=2005;// For BigInteger
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

int rev(int x) {
    int ret=0;
    while (x) {
        int tmp=x%10;
        x/=10;
        ret=ret*10+tmp;
    }   
    return ret;
}

map<pair<int,int>,int> mpx,mpy;
map<int,int> mp;
int maxx,maxy,w;

void init() {
    rep(i,1,200000) mp[i]=rev(i);
//  rep(i,1,20) cout << i << ' ' << mp[i]  << endl;
}

pair<int,int> get(int x,int flag) {
    int tmp=__gcd(x,mp[x]);
    pair<int,int> ret;
    ret.first=x/tmp;
    ret.second=mp[x]/tmp;
    if (flag) swap(ret.first,ret.second);
    return ret;
}

int main(){
#ifdef LOCAL
    freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
//  freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
    init();
    while (SIII(maxx,maxy,w)==3) {
        mpx.clear();mpy.clear();
        rep(i,1,maxy) {
            pair<int,int> tmp=get(i,0);
            mpy[tmp]++;
        }
        int tt=0;
        int ok=0,x,ans;
        rep(i,1,maxx) {
            pair<int,int> tmp=get(i,1);
            tt+=mpy[tmp];
            mpx[tmp]++;
            if (tt>=w){
                ok=1;
                x=i;
                break;
            }
        }
        if (!ok) {
            cout << -1 << endl;
            continue;
        }
        int y=maxy;
        ans=y*x;
        int a=x,b=y;
        while (y) {
            mpy[get(y,0)]--;
            tt-=mpx[get(y,0)];
            ok=0;
            while (x<maxx && tt<w) {
                x++;
                pair<int,int> tmp=get(x,1);
                mpx[tmp]++;
                tt+=mpy[tmp];
            }
            y--;
            if (tt>=w && x*y<ans) {
                ans=x*y;
                a=x;b=y;
            }
            if (x>=maxx) break;
        }
        cout << a << ' ' << b << endl;
    }
    return 0;
}