URAL 1554

题目大意：给出算数函数f(n)（mod 2007）求它的逆函数g(n)（mod 2007）。

Time Limit:2000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u

数据规模：1<=N<=10^4。

理论基础：狄利克雷卷积：参见链接1。

题目分析：是一道阅读题，看完链接你就明白了，时间限制较松，枚举N的所有因子即可。然后用公式法求解。一些基本的同余理论，在本博客置顶文章中已经讲述，这里不再累赘。

代码如下：

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<ctime>
using namespace std;
typedef double db;
#define DBG 1
#define maa (1<<31)
#define mii ((1<<31)-1)
#define sl(c) ((int)(c).size())    //取字符串长度；
#define forl(i, a, b) for(int i = (a); i <  (b); ++i)    //不带边界值循环，正序
#define forle(i, a, b) for(int i = (a); i <= (b); ++i)   //带边界值循环，正序
#define forh(i, a, b) for(int i = (a); i >  (b); --i)     //不带边界值，逆序
#define forhe(i, a, b) for(int i = (a); i >= (b); --i)        //带边界值，逆序
#define forlc(i, a, b) for(int i##_b = (b), i = (a); i <  i##_b; ++i)  //带别名的循环，用于不可改变值
#define forlec(i, a, b) for(int i##_b = (b), i = (a); i <= i##_b; ++i)
#define forgc(i, a, b) for(int i##_b = (b), i = (a); i >  i##_b; --i)
#define forgec(i, a, b) for(int i##_b = (b), i = (a); i >= i##_b; --i)
#define forall(i, v   )  forl(i, 0, sz(v))   //循环所有
#define forallc(i, v   ) forlc(i, 0, sz(v))
#define forlla(i, v   ) forhe(i, sz(v)-1, 0)
#define forls(i, n, a, b) for(int i = a; i != b; i = n[i])   //搜表用
#define rep(n)  for(int               repp = 0; repp <    (n); ++repp)
#define repc(n) for(int repp_b = (n), repp = 0; repp < repp_b; ++repp)
#define rst(a, v) memset(a, v, sizeof a)   //把字符v填充到a  reset 重置
#define cpy(a, b) memcpy(a, b, sizeof a)   //copy b 的sizeof（a）个字符到a
#define rstn(a, v, n) memset(a, v, (n)*sizeof((a)[0]))  //把字符v填充到a[n]之前的字节
#define cpyn(a, b, n) memcpy(a, b, (n)*sizeof((a)[0]))    //copy b 的 n 个字符到a
#define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); }  //调试
#define dout DBG && cout << __LINE__ << ">>| "
#define pr(x) #x"=" << (x) << " | "
#define mk(x) DBG && cout << __LINE__ << "**| "#x << endl
#define pra(arr, a, b)  if(DBG) {\
    dout<<#arr"[] |" <<endl; \
    forlec(i, a, b) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \
    if((b-a+1)%8) puts("");\
}                                                             //数列查看
#define rd(type, x) type x; cin >> x   //读数
inline int     rdi() { int d; scanf("%d", &d); return d; }
inline char    rdc() { scanf(" "); return getchar(); }
inline string  rds() { rd(string, s); return s; }
inline db rddb() { db d; scanf("%lf", &d); return d; }
template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; }
template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; }

typedef long long LL;
typedef long unsigned LU;
int f[10001],g[10001],n,fi;
void breakup(int n,int d[],int &nd)
{
    forlec(i,1,sqrt((float)n))
    {
        if(n%i==0)
        {
            if(n/i!=i&&i!=1)
            {
                d[nd++]=i;
                d[nd++]=n/i;
            }
            else d[nd++]=i;
        }
    }
}
int mod(int a,int b)
{
    if(a>=0)return a%b;
    return (a%b+b)%b;
}
void solve(int i)
{
    if(i==1)
    {
        g[i]=1/f[1];
        return;
    }
    int d[100],nd=0,sum=0;
    breakup(i,d,nd);
    sort(d,d+nd);
    forl(j,0,nd)
    {
        sum=sum+g[d[j]]*f[i/d[j]];
    }
    g[i]=mod(-g[1]*sum,2007);
}
int main()
{
    scanf("%d",&n);
    forle(i,1,n)
    {
        f[i]=rdi();
    }
    forle(i,1,n)
    {
        solve(i);
        if(i==1)printf("%d",g[i]);
        else printf(" %d",g[i]);
    }
    printf("\n");
	return 0;
}

    其中，用mod（）函数防止a%b返回负值的情况。

参考文献：

http://zh.wikipedia.org/wiki/%E7%8B%84%E5%88%A9%E5%85%8B%E9%9B%B7%E5%8D%B7%E7%A9%8D

by:Jsun_moon http://blog.csdn.net/Jsun_moon