题意:给定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;
}