UVA-1635(数论)(二项式定理性质+递推+唯一分解定理)

题意:给定n个从1递增的数,将这些数两两相加从而得到一个最终的数,并将该数对m取余,问该余数与哪些数无关?
题解:无关即该数前面的系数是m的倍数即可。
多些几组数据容易发现该数列的系数是一个二项式系数,又因为该数过于大不能直接表示,故可以用唯一分解定理来辨别,另辨别的时候用到了二项式定理的递推式。

hint:可能是因为用到了map的原因,还有就是循环里面写的内容太多了,跑了2000多ms,还是满惊险的。
附上AC代码:

#include <bits/stdc++.h>
#define FOPI freopen("INPUT.TXT", "r", stdin)
#define DOPI freopen("OUTPUT.TXT", "w", stdout)
using namespace std;
typedef long long int ll;
const int ind=0x3f3f3f3f,N=1e5+10;
const ll inlld=0x3f3f3f3f3f3f3f3f,mod=998244353;
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
ll lcm(ll a,ll b){return a*b/gcd(a,b);}
typedef pair<int,int> p;
ll n,m;
int prime[N];
int e[N];
map<int,int> mp;
vector<int>ans;
int main()
{
    ios::sync_with_stdio(false);
    while(cin>>n>>m){
        int num=0;
        ans.clear();
        memset(e,0,sizeof(e));
        memset(prime,0,sizeof(prime));
        n-=1;
        mp.clear();
        int mm=sqrt(m+0.5);
        int j=1;
        for(int i=2;i<=mm;i++){
            if(m%i==0){
                prime[j]=i;
                while(m%i==0){
                    e[j]++;m/=i;
                }
                j++;
            }
            if(m==1)break;
        }
        if(m>1){
            prime[j]=m;e[j]++;j++;
        }
        mp[1]++;
        for(int i=1;i<=n;i++){
            int k=i;int tm=n-k+1;
            for(int t=2;t<=sqrt(k+0.5);t++){
                while(k%t==0){
                    mp[t]--;
                    k/=t;
                }
                if(k==1)break;
            }
            if(k>1)mp[k]--;
            for(int t=2;t<=sqrt(tm+0.5);t++){
                while(tm%t==0){
                    mp[t]++;
                    tm/=t;
                }
                if(tm==1)break;
            }
            if(tm>1)mp[tm]++;int fg=1;
            for(int t=1;t<j;t++){
                if(mp[prime[t]]<e[t]){
                    fg=0;break;
                }
            }
            if(fg){
                ans.push_back(i+1);num++;
            }
        }
        cout<<num<<endl;
        for(int i=0;i<ans.size();i++){
            if(!i)cout<<ans[i];
            else cout<< ' ' << ans[i];
        }cout<<endl;
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值