标签：des   style   blog   http   io   ar   color   os   sp   

Ignatius‘s puzzle

11
100
9999

22
no
43

题意：f(x)=5*x^13+13*x^5+k*a*x,找出最小的a，使得对于给定的k，任意x下f（x）能被65整除。

思路：打表找规律。

你会发现，（5*x^13+13*x^5）%65每65个x成循环。要使结果成立（k*a*x）%65，也得每65个循环，且a的值在0-64之间，所以就有如下代码-> _ ->

代码：

/************************************************
* Author: Ac_sorry
* File:
* Create Date:
* Motto: One heart One life
* CSDN: http://blog.csdn.net/code_or_code
*************************************************/

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<vector>
#include<string>
#include<map>
#include<utility>
#include<queue>
#include<stack>
#define INF 0x3f3f3f3f
#define MOD 1000000007
#define seed_ 131
#define eps 1e-8
#define mem(a,b) memset(a,b,sizeof a)
#define w(i) T[i].w
#define ls(i) T[i].ls
#define rs(i) T[i].rs
using namespace std;
typedef long long LL;

const int N=100010;

int a[100];

int main()
{
    //int T;scanf("%d",&T);
    //int AC=0;
    for(int x=1;x<=65;x++)
        a[x]=(5*x%65*x*x%65*x*x%65*x*x%65*x*x%65*x*x%65*x*x%65+13*x%65*x*x%65*x*x%65)%65;
    int k;
    while(scanf("%d",&k)==1)
    {
        int ok,ans;
        for(int aa=0;aa<65;aa++)
        {
            ok=1;
            for(int i=1;i<=65;i++)
            {
                if((a[i]+k*i%65*aa%65)%65!=0)
                {
                    ok=0;
                    break;
                }
            }
            if(ok)
            {
                ans=aa;
                break;
            }
        }
        if(ok)
            printf("%d\n",ans);
        else
            printf("no\n");
    }
    return 0;
}

HDU_1098 Ignatius's puzzle[规律题]

标签：des   style   blog   http   io   ar   color   os   sp   

原文地址：http://blog.csdn.net/code_or_code/article/details/41870841